Search

We present a high-order boundary integral equation (BIE) method for the frequency-domain acoustic scattering of a point source by a singly-periodic, infinite, corrugated boundary. We apply it to the accurate numerical study of acoustic radiation in the neighborhood of a sound-hard two-dimensional staircase modeled after the El Castillo pyramid. Such staircases support trapped waves which travel along the surface and decay exponentially away from it. We use the array scanning method (Floquet--Bloch transform) to recover the scattered field as an integral over the family of quasiperiodic solutions parameterized by their on-surface wavenumber. Each such BIE solution requires the quasiperiodic Green's function, which we evaluate using an efficient integral representation of lattice sum coefficients. We avoid the singularities and branch cuts present in the array scanning integral by complex contour deformation. For each frequency, this enables a solution accurate to around 10 digits in a couple of seconds. We propose a residue method to extract the limiting powers carried by trapped modes far from the source. Finally, by computing the trapped mode dispersion relation, we use a simple ray model to explain an observed acoustic "raindrop" effect (chirp-like time-domain response)., Comment: 25 pages, 11 figures. Submitted to JCP

Ion implantation has been a key technology for the controlled surface modification of materials in microelectronics and generally, for tribology, biocompatibility, corrosion resistance and many more. To form shallow junctions in Ge is a challenging task. In this work the formation and accumulation of shallow damage profiles was studied by in-situ spectroscopic ellipsometry (SE) for the accurate tracking and evaluation of void and damage fractions in crystalline Ge during implantation of 200-keV Sb ions with a total fluence up to 1E16 cm-2 and an ion flux of 2.1E12 cm-2 s-1. The consecutive stages of damage accumulation were identified using optical multi-layer models with quantitative parameters of the thickness of modified layers as well as the volume fractions of amorphized material and voids. The effective size of damaged zones formed from ion tracks initiated by individual bombarding ions can be estimated by numerical simulation compared with the dynamics of damage profiles measured by ion beam analysis and ellipsometry. According to our observations, the formation of initial partial disorder was followed by complete amorphization and void formation occurring at the fluence of about 1E15 cm-2, leading to a high volume fraction of voids and a modified layer thickness of approx. 200 nm by the end of the irradiation process. This agrees with the results of numerical simulations and complementary scanning electron microscopy (SEM) measurements. In addition, we found a quasi-periodic time dependent behavior of amorphization and void formation represented by alternating accelerations and decelerations of different reorganization processes, respectively., Comment: 23 pages, 8 figures

We introduce an efficient numerical method for second order linear ODEs whose solution may vary between highly oscillatory and slowly changing over the solution interval. In oscillatory regions the solution is generated via a nonoscillatory phase function that obeys the nonlinear Riccati equation. We propose a defect-correction iteration that gives an asymptotic series for such a phase function; this is numerically approximated on a Chebyshev grid with a small number of nodes. For analytic coefficients we prove that each iteration, up to a certain maximum number, reduces the residual by a factor of order of the local frequency. The algorithm adapts both the step size and the choice of method, switching to a conventional spectral collocation method away from oscillatory regions. In numerical experiments we find that our proposal outperforms other state-of-the-art oscillatory solvers, most significantly at low-to-intermediate frequencies and at low tolerances, where it may use up to $10^6$ times fewer function evaluations. Even in high frequency regimes, our implementation is on average 10 times faster than other specialized solvers., Comment: 25 pages, 5 figures, 2 tables. Submitted to SINUM. Code available at https://github.com/fruzsinaagocs/riccati

We discuss the challenges of motivating, constructing, and quantizing a canonically normalized inflationary perturbation in spatially curved universes. We show that this has historically proved challenging due to the interaction of nonadiabaticity with spatial curvature. We construct a novel curvature perturbation that is canonically normalized in the sense of its equation of motion and is unique up to a single scalar parameter. With this construction it becomes possible to set initial conditions invariant under canonical transformations, overcoming known ambiguities in the literature. This corrected quantization has potentially observational consequences via modifications to the primordial power spectrum at large angular scales, as well as theoretical implications for quantization procedures in curved cosmologies filled with a scalar field., Comment: 15 pages, 3 figures; v3: As published in PRD

We investigate the effects of non-zero spatial curvature on cosmic inflation in the light of cosmic microwave background (CMB) anisotropy measurements from the Planck 2018 legacy release and from the 2015 observing season of BICEP2 and the Keck Array. Even a small percentage of non-zero curvature today would significantly limit the total number of e-folds of the scale factor during inflation, rendering just-enough inflation scenarios with a kinetically dominated or fast-roll stage prior to slow-roll inflation more likely. Finite inflation leads to oscillations and a cutoff towards large scales in the primordial power spectrum and curvature pushes them into the CMB observable window. Using nested sampling, we carry out Bayesian parameter estimations and model comparisons taking into account constraints from reheating and horizon considerations. We confirm the preference of CMB data for closed universes with Bayesian odds of over $100:1$ and with a posterior on the curvature density parameter of $\Omega_{K,0}=-0.051\pm0.017$ for a curvature extension of LCDM and $\Omega_{K,0}=-0.031\pm0.014$ for Starobinsky inflation. Model comparisons of various inflation models give similar results as for flat universes with the Starobinsky model outperforming most other models., Comment: 46 pages, 32 figures, 3 tables, submitted to PRD, from chapter 5 of L. T. Hergt's doctoral thesis at https://www.repository.cam.ac.uk/handle/1810/323129, data available at https://zenodo.org/record/6547872

Context:Optical interferometry is at a key development stage. ESO's VLTI has established a stable, robust infrastructure for long-baseline interferometry for general astronomical observers. The present second-generation instruments offer a wide wavelength coverage and improved performance. Their sensitivity and measurement accuracy lead to data and images of high reliability. Aims:We have developed MATISSE, the Multi AperTure mid-Infrared SpectroScopic Experiment, to access high resolution imaging in a wide spectral domain and explore topics such: stellar activity and mass loss; planet formation and evolution in the gas and dust disks around young stars; accretion processes around super massive black holes in AGN. Methods:The instrument is a spectro-interferometric imager covering three atmospheric bands (L,M,N) from 2.8 to 13.0 mu, combining four optical beams from the VLTI's telscopes. Its concept, related observing procedure, data reduction and calibration approach are the product of 30 years of instrumental research. The instrument utilizes a multi-axial beam combination that delivers spectrally dispersed fringes. The signal provides the following quantities at several spectral resolutions: photometric flux, coherent fluxes, visibilities, closure phases, wavelength differential visibilities and phases, and aperture-synthesis imaging. Results:We provide an overview of the physical principle of the instrument and its functionalities, the characteristics of the delivered signal, a description of the observing modes and of their performance limits. An ensemble of data and reconstructed images are illustrating the first acquired key observations. Conclusion:The instrument has been in operation at Cerro Paranal, ESO, Chile since 2018, and has been open for science use by the international community since April 2019. The first scientific results are being published now., Comment: 24 pages, 26 figures, submitted to Astronomy & Astrophysics

The nature of circumstellar envelopes (CSE) around Cepheids is still a matter of debate. The physical origin of their infrared (IR) excess could be either a shell of ionized gas, or a dust envelope, or both. This study aims at constraining the geometry and the IR excess of the environment of the long-period Cepheid $\ell$ Car (P=35.5 days) at mid-IR wavelengths to understand its physical nature. We first use photometric observations in various bands and Spitzer Space Telescope spectroscopy to constrain the IR excess of $\ell$ Car. Then, we analyze the VLTI/MATISSE measurements at a specific phase of observation, in order to determine the flux contribution, the size and shape of the environment of the star in the L band. We finally test the hypothesis of a shell of ionized gas in order to model the IR excess. We report the first detection in the L band of a centro-symmetric extended emission around l Car, of about 1.7$R_\star$ in FWHM, producing an excess of about 7.0\% in this band. In the N band, there is no clear evidence for dust emission from VLTI/MATISSE correlated flux and Spitzer data. On the other side, the modeled shell of ionized gas implies a more compact CSE ($1.13\pm0.02\,R_\star$) and fainter (IR excess of 1\% in the L band). We provide new evidences for a compact CSE of $\ell$ Car and we demonstrate the capabilities of VLTI/MATISSE for determining common properties of CSEs. While the compact CSE of $\ell$ Car is probably of gaseous nature, the tested model of a shell of ionized gas is not able to simultaneously reproduce the IR excess and the interferometric observations. Further Galactic Cepheids observations with VLTI/MATISSE are necessary for determining the properties of CSEs, which may also depend on both the pulsation period and the evolutionary state of the stars., Comment: 13 pages, 8 figures, accepted in Astronomy and Astrophysics

Ground-based thermal-infrared observations have a unique scientific potential, but are also extremely challenging due to the need to accurately subtract the high thermal background. Since the established techniques of chopping and nodding need to be modified for observations with the future mid-infrared ELT imager and spectrograph (METIS), we investigate the sources of thermal background subtraction residuals. Our aim is to either remove or at least minimise the need for nodding in order to increase the observing efficiency for METIS. To this end we need to improve our knowledge about the origin of chop residuals and devise observing methods to remove them most efficiently, i.e. with the slowest possible nodding frequency. Thanks to dedicated observations with VLT/VISIR and GranTeCan/CanariCam, we have successfully traced the origin of three kinds of chopping residuals to (1) the entrance window, (2) the spiders and (3) other warm emitters in the pupil, in particular the VLT M3 mirror cell in its parking position. We conclude that, in order to keep chopping residuals stable over a long time (and therefore allow for slower nodding cycles), the pupil illumination needs to be kept constant, i.e. (imaging) observations should be performed in pupil-stabilised, rather than field-stabilised mode, with image de-rotation in the post-processing pipeline. This is now foreseen as the default observing concept for all METIS imaging modes., Comment: 17 pages, 9 figures, submitted to the SPIE proceedings as paper no. 11447-403

Physical theories that depend on many parameters or are tested against data from many different experiments pose unique challenges to statistical inference. Many models in particle physics, astrophysics and cosmology fall into one or both of these categories. These issues are often sidestepped with statistically unsound ad hoc methods, involving intersection of parameter intervals estimated by multiple experiments, and random or grid sampling of model parameters. Whilst these methods are easy to apply, they exhibit pathologies even in low-dimensional parameter spaces, and quickly become problematic to use and interpret in higher dimensions. In this article we give clear guidance for going beyond these procedures, suggesting where possible simple methods for performing statistically sound inference, and recommendations of readily-available software tools and standards that can assist in doing so. Our aim is to provide any physicists lacking comprehensive statistical training with recommendations for reaching correct scientific conclusions, with only a modest increase in analysis burden. Our examples can be reproduced with the code publicly available at https://doi.org/10.5281/zenodo.4322283., Comment: 15 pages, 4 figures. extended discussions. closely matches version accepted for publication

Context. The inner few au region of planet-forming disks is a complex environment. High angular resolution observations have a key role in understanding the disk structure and the dynamical processes at work. Aims. In this study we aim to characterize the mid-infrared brightness distribution of the inner disk of the young intermediate-mass star HD 163296, from VLTI/MATISSE observations. Methods. We use geometric models to fit the data. Our models include a smoothed ring, a flat disk with inner cavity, and a 2D Gaussian. The models can account for disk inclination and for azimuthal asymmetries as well. We also perform numerical hydro-dynamical simulations of the inner edge of the disk. Results. Our modeling reveals a significant brightness asymmetry in the L-band disk emission. The brightness maximum of the asymmetry is located at the NW part of the disk image, nearly at the position angle of the semimajor axis. The surface brightness ratio in the azimuthal variation is $3.5 \pm 0.2$. Comparing our result on the location of the asymmetry with other interferometric measurements, we confirm that the morphology of the $r<0.3$ au disk region is time-variable. We propose that this asymmetric structure, located in or near the inner rim of the dusty disk, orbits the star. For the physical origin of the asymmetry, we tested a hypothesis where a vortex is created by Rossby wave instability, and we find that a unique large scale vortex may be compatible with our data. The half-light radius of the L-band emitting region is $0.33\pm 0.01$ au, the inclination is ${52^\circ}^{+5^\circ}_{-7^\circ}$, and the position angle is $143^\circ \pm 3^\circ$. Our models predict that a non-negligible fraction of the L-band disk emission originates inside the dust sublimation radius for $\mu$m-sized grains. Refractory grains or large ($\gtrsim 10\ \mu$m-sized) grains could be the origin for this emission., Comment: accepted for publication in A&A

We present a method to construct a continuous extension (otherwise known as dense output) for a numerical routine in the special case of the numerical solution being a scalar-valued function exhibiting rapid oscillations. Such cases call for numerical routines that make use of the known global behaviour of the solution, one example being methods using asymptotic expansions to forecast the solution at each step of the independent variable. An example is oscode, numerical routine which uses the Wentzel-Kramers-Brillouin (WKB) approximation when the solution oscillates rapidly and otherwise behaves as a Runge-Kutta (RK) solver. Polynomial interpolation is not suitable for producing the solution at an arbitrary point mid-step, since efficient numerical methods based on the WKB approximation will step through multiple oscillations in a single step. Instead we construct the continuous solution by extending the numerical quadrature used in computing a WKB approximation of the solution with no additional evaluations of the differential equation or terms within, and provide an error estimate on this dense output. Finally, we draw attention to previous work on the continuous extension of Runge-Kutta formulae, and construct an extension to a RK method based on Gauss--Lobatto quadrature nodes, thus describing how to generate dense output from each of the methods underlying oscode., Comment: 10 pages, 5 figures, 4 tables. Submitted to PRResearch

We investigate the transformation of initial conditions for primordial curvature perturbations under two types of transformations of the associated action: simultaneous redefinition of time and the field to be quantised, and the addition of surface terms. The latter encompasses all canonical transformations, whilst the time- and field-redefinition is a distinct, non-canonical transformation since the initial and destination systems use different times. Actions related to each other via such transformations yield identical equations of motion and preserve the commutator structure. They further preserve the time-evolution of expectation values of quantum operators unless the vacuum state also changes under the transformation. These properties suggest that it is of interest to investigate vacuum prescriptions that also remain unchanged under canonical transformations. We find that initial conditions derived via minimising the vacuum expectation value of the Hamiltonian and those obtained using the Danielsson vacuum prescription are not invariant under these transformations, whereas those obtained by minimising the local energy density are. We derive the range of physically distinct initial conditions obtainable by Hamiltonian diagonalisation, and illustrate their effect on the scalar primordial power spectrum and the Cosmic Microwave Background under the just enough inflation model. We also generalise the analogy between the dynamics of a quantum scalar field on a curved, time-dependent spacetime and the gauge-invariant curvature perturbation. We argue that the invariance of the vacuum prescription obtained by minimising the renormalised stress--energy tensor should make it the preferred procedure for setting initial conditions for primordial perturbations. All other procedures reviewed in this work yield ambiguous initial conditions, which is problematic both in theory and practice., Comment: Revisions after review added. 17 pages, 8 figures. Accepted to PRD

We present a novel numerical routine (oscode) with a C++ and Python interface for the efficient solution of one-dimensional, second-order, ordinary differential equations with rapidly oscillating solutions. The method is based on a Runge-Kutta-like stepping procedure that makes use of the Wentzel-Kramers-Brillouin (WKB) approximation to skip regions of integration where the characteristic frequency varies slowly. In regions where this is not the case, the method is able to switch to a made-to-measure Runge-Kutta integrator that minimises the total number of function evaluations. We demonstrate the effectiveness of the method with example solutions of the Airy equation and an equation exhibiting a burst of oscillations, discussing the error properties of the method in detail. We then show the method applied to physical systems. First, the one-dimensional, time-independent Schr\"odinger equation is solved as part of a shooting method to search for the energy eigenvalues for a potential with quartic anharmonicity. Then, the method is used to solve the Mukhanov-Sasaki equation describing the evolution of cosmological perturbations, and the primordial power spectrum of the perturbations is computed in different cosmological scenarios. We compare the performance of our solver in calculating a primordial power spectrum of scalar perturbations to that of BINGO, an efficient code specifically designed for such applications., Comment: Physical Review Research, accepted. 23 pages, 15 figures. The associated code is available online at https://github.com/fruzsinaagocs/oscode

Data visualisation assists domain experts in understanding their data and helps them make critical decisions. Enhancing their cognitive insight essentially relies on the capability of combining domain-specific semantic information with concepts extracted out of the data and visualizing the resulting networks. Data scientists have the challenge of providing tools able to handle the overall network lifecycle. In this paper, we present how the combination of two powerful technologies namely the REST architecture style and JSON Schema/JSON Meta Schema enable data scientists to use a RESTful web service that permits the construction of knowledge graphs, one of the preferred representations of large and semantically rich networks., Comment: 5 pages; 6 figures

Human reasoning in visual analytics of data networks relies mainly on the quality of visual perception and the capability of interactively exploring the data from different facets. Visual quality strongly depends on networks' size and dimensional complexity while network exploration capability on the intuitiveness and expressiveness of user frontends. The approach taken in this paper aims at addressing the above by decomposing data networks into multiple networks of smaller dimensions and building an interactive graph query language that supports full navigation across the sub-networks. Within sub-networks of reduced dimensionality, structural abstraction and semantic techniques can then be used to enhance visual perception further., Comment: 7 pages, 9 figures

An experiment conducted in the framework of the EUHIT project and designed to characterize large scale structures in an adverse pressure gradient boundary layer flow is presented. Up to 16 sCMOS cameras were used in order to perform large scale turbulent boundary layer PIV measurements with a large field of view and appropriate spatial resolution. To access the span-wise / wall-normal signature of the structures as well, stereoscopic PIV measurements in span-wise/wall-normal planes were performed at specific stream-wise locations. To complement these large field of view measurements, long-range micro-PIV, time resolved near wall velocity profiles and film-based measurements were performed in order to determine the wall-shear stress and its fluctuations at some specific locations along the model., Comment: 50 pages