Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
25works
0followers
16topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

25 published item(s)

preprint2022arXiv

An Efficient ADER-DG Local Time Stepping Scheme for 3D HPC Simulation of Seismic Waves in Poroelastic Media

Many applications from geosciences require simulations of seismic waves in porous media. Biot's theory of poroelasticity describes the coupling between solid and fluid phases and introduces a stiff source term, thereby increasing computational cost and motivating efficient methods utilising High-Performance Computing. We present a novel realisation of the discontinuous Galerkin scheme with Arbitrary DERivative time stepping (ADER-DG) that copes with stiff source terms. To integrate this source term with a reasonable time step size, we use an element-local space-time predictor, which needs to solve medium-sized linear systems - with 1000 to 10000 unknowns - in each element update (i.e., billions of times). We present a novel block-wise back-substitution algorithm for solving these systems efficiently. In comparison to LU decomposition, we reduce the number of floating-point operations by a factor of up to 25. The block-wise back-substitution is mapped to a sequence of small matrix-matrix multiplications, for which code generators are available to generate highly optimised code. We verify the new solver thoroughly in problems of increasing complexity. We demonstrate high-order convergence for 3D problems. We verify the correct treatment of point sources, material interfaces and traction-free boundary conditions. In addition, we compare against a finite difference code for a newly defined layer over half-space problem. We find that extremely high accuracy is required to resolve the slow P-wave at a free surface, while solid particle velocities are not affected by coarser resolutions. By using a clustered local time stepping scheme, we reduce time to solution by a factor of 6 to 10 compared to global time stepping. We conclude our study with a scaling and performance analysis, demonstrating our implementation's efficiency and its potential for extreme-scale simulations.

preprint2022arXiv

Collective photon emission patterns from two atoms in free space

Modification of spontaneous decay in space and time is a central topic of quantum physics. It has been predominantly investigated in the context of cavity quantum electrodynamics (QED), gaining new interest recently in the domain of nano-optics. Beyond cavity-QED, spontaneous emission may be modified also in free space due to correlations among the photon emitters, a phenomenon known as super- and sub-radiance. Correlations may stem either from direct interactions between the particles, from long-range exchange of photons, or by measuring single photons in a common mode. Yet, the genuine spatial spontaneous emission pattern of an atomic ensemble in an entangled quantum state has not been observed so far, due to the lack of ultra-fast cameras with high spatial resolution suited for recording single photons from single atoms. Preparing two trapped ions in free space in entangled Dicke states via photon detection, we study the resulting collective spontaneous emission patterns. Depending on the symmetry of the Dicke states, associated with the direction of detection of the first state-determining photon, we observe fundamentally different emission patterns for the subsequently scattered photon, including super- and sub-radiance. Our results demonstrate that the detection of a single photon can profoundly modify the collective emission of an atomic array, here represented by its most elementary building block of two atoms in free space.

preprint2022arXiv

High-sensitivity quantum sensing with pump-enhanced spontaneous parametric down-conversion

Recent years have seen the development of quantum sensing concepts utilizing nonlinear interferometers based on correlated photon pairs generated by spontaneous parametric down-conversion (SPDC). Using SPDC far from frequency degeneracy allows a 'division of labor' between the mid-infrared photon for strongest sample interaction and the correlated near-infrared photon for low-noise detection. The small number of photons provided by SPDC and the resulting inferior signal-to-noise ratio are, however, a limiting factor preventing wide applicability of the novel sensing concept. Here, we demonstrate a nonlinear interferometer based on pump-enhanced SPDC with strongly improved emission rates, but maintaining broadband, spontaneous emission. For validation of the concept, we demonstrate high-resolution mid-infrared spectroscopy with near-infrared detection, showcasing the improved accuracy. Although the number of mid-infrared photons is about five orders of magnitude smaller than in classical spectrometers, the sensitivity of the quantum spectrometer becomes comparable, marking an essential step toward real-world applications.

preprint2022arXiv

Locating dust and molecules in the inner circumstellar environment of R~Sculptoris with MATISSE

AGB stars are one of the main sources of dust production in the Galaxy. However, it is not clear what this process looks like and where the dust is condensing in the circumstellar environment. By characterizing the location of the dust and the molecules in the close environment of an AGB star, we aim to achieve a better understanding the history of the dust formation process. We observed the carbon star R Scl with the VLTI-MATISSE instrument in L- and N-bands. The high angular resolution of the VLTI observations, combined with a large uv-plane coverage allowed us to use image reconstruction methods. To constrain the dust and molecules' location, we used two different methods: MIRA image reconstruction and the 1D code RHAPSODY. We found evidence of C2H2 and HCN molecules between 1 and 3.4 Rstar which is much closer to the star than the location of the dust (between 3.8 and 17.0 Rstar). We also estimated a mass-loss rate of 1.2+-0.4x10-6 Msun per yr. In the meantime, we confirmed the previously published characteristics of a thin dust shell, composed of amorphous carbon (amC) and silicon carbide (SiC). However, no clear SiC feature has been detected in the MATISSE visibilities. This might be caused by molecular absorption that can affect the shape of the SiC band at 11.3 micron. The appearance of the molecular shells is in good agreement with predictions from dynamical atmosphere models. For the first time, we co-located dust and molecules in the environment of an AGB star. We confirm that the molecules are located closer to the star than the dust. The MIRA images unveil the presence of a clumpy environment in the fuzzy emission region beyond 4.0 Rstar. Furthermore, with the available dynamic range and angular resolution, we did not detect the presence of a binary companion. Additional observations combining MATISSE and SAM-VISIR instrument should enable this detection in future studies.

preprint2022arXiv

Magnetic fields and outflows in the large Bok globule CB 54

We have observed the large Bok globule CB 54 in 850$μ$m polarised light using the POL-2 polarimeter on the James Clerk Maxwell Telescope (JCMT). We find that the magnetic field in the periphery of the globule shows significant, ordered deviation from the mean field direction in the globule centre. This deviation appears to correspond with the extended but relatively weak $^{12}$CO outflow emanating from the Class 0 sources at the centre of the globule. Energetics analysis suggests that if the outflow is reshaping the magnetic field in the globule&#39;s periphery, then we can place an upper limit of $<27\,μ$G on the magnetic field strength in the globule&#39;s periphery. Comparison with archival Planck and CARMA measurements shows that the field in the centre of the globule is consistent over several orders of magnitude in size scale, and oriented parallel to the density structure in the region in projection. We thus hypothesise that while non-thermal motions in the region may be sub-Alfvénic, the magnetic field is subdominant to gravity over a wide range of size scales. Our results suggest that even a relatively weak outflow may be able to significantly reshape magnetic fields in star-forming regions on scales $> 0.1$ pc.

preprint2022arXiv

Polarimetric investigation of selected cloud compositions in exoplanetary atmospheres

We investigated the impact of selected cloud condensates in exoplanetary atmospheres on the polarization of scattered stellar radiation. We considered a selection of 25 cloud condensates that are expected to be present in extrasolar planetary atmospheres. Using the three-dimensional Monte Carlo radiative transfer code POLARIS and assuming Mie scattering theory, we calculated and studied the net polarization of scattered radiation as a function of planetary phase angle at optical to near-infrared wavelengths. In addition to the well-known characteristics in the state of polarization, such as the rainbow determined by the real part of the refractive index, the behavior of the underlying imaginary part of the refractive index causes an increase or decrease in the degree of polarization and a change of sign in the polarization at a characteristic wavelength. In contrast to Al$_2$O$_3$ and MgFeSiO$_4$, clouds composed of SiO, MnS, Na$_2$S, or ZnS produce a rapidly decreasing degree of polarization with increasing wavelength in the context of an exoplanetary atmosphere. Furthermore, the sign of the polarization changes at a wavelength of about 0.5 $μ$m to 0.6 $μ$m, depending on the specific cloud condensate. The resulting net polarization is mainly positive for cloud compositions with large imaginary parts of the refractive index, such as Fe, FeS, and FeO. In addition, for Fe and FeS clouds, the maximum degree of polarization at long wavelengths is shifted to larger phase angles than for FeO. We found that most of these cloud condensates are distinguishable from each other due to their unique wavelength-dependent complex refractive index. In particular, an increase or decrease of the net polarization as a function of wavelength and a change of sign in the polarization at specific wavelengths are important features for characterizing cloud compositions in exoplanetary atmospheres.

preprint2022arXiv

Radiative transfer simulations for in-situ particle size diagnostic in reactive, particle growing plasmas

When considering particles produced in reactive plasmas, their basic properties, such as refractive index and grain size often need to be known. They can be constrained both ex-situ, e.g., by microscopy, and in-situ by polarimetry, i.e., analyzing the polarization state of scattered light. Polarimetry has the advantage of temporal resolution and real-time measurement, but the analysis is often limited by the assumption of single scattering and thus optically thin dust clouds. This limits the investigation of the growth process typically to grain sizes smaller than about 200 nm. Using 3D polarized radiative transfer simulations, however, it is possible to consider multiple scattering and to analyze the properties of dense particle clouds. We study the impact of various properties of dust clouds on the scattering polarization, namely the optical depth of the cloud, the spatial density distribution of the particles, their refractive index as well as the particle size dispersion. We find that ambiguities can occur regarding optical depth and spatial density distribution as well as regarding refractive index and particle size dispersion. Determining the refractive index correctly is especially important as it has a strong impact on the derived particle sizes. With this knowledge, we are able to design an in-situ diagnostics strategy for the investigation of the particle growth process based on radiative transfer simulations which are used to model the polarization over the whole growth process. The application of this strategy allows us for the first time to analyze the polarization measured during a growth experiment in a reactive argon-acetylene plasma for particle radii up to 280 nm.

preprint2022arXiv

Robust Polarization Gradient Cooling of Trapped Ions

We implement three-dimensional polarization gradient cooling of trapped ions. Counter-propagating laser beams near $393\,$nm impinge in lin$\,\perp\,$lin configuration, at a frequency below the S$_{1/2}$ to P$_{3/2}$ resonance in $^{40}$Ca$^+$. We demonstrate mean phonon numbers of $5.4(4)$ at a trap frequency of $2π\times 285\,$kHz and $3.3(4)$ at $2π\times480\,$kHz, in the axial and radial directions, respectively. Our measurements demonstrate that cooling with laser beams detuned to lower frequencies from the resonance is robust against an elevated phonon occupation number, and thus works well for an initial ion motion far out of the Lamb-Dicke regime, for up to four ions, and for a micromotion modulation index $β\leq 0.1$. Still, we find that the spectral impurity of the laser field influences both, cooling rates and cooling limits. Thus, a Fabry-Pérot cavity filter is employed for efficiently suppressing amplified spontaneous emission of the diode laser.

preprint2022arXiv

Synthetic observations using POLARIS: an application to simulations of massive prestellar cores

Young massive stars are usually found embedded in dense and massive molecular clumps which are known for being highly obscured and distant. During their formation process, the degree of deuteration can be used as a potential indicator of the very early formation stages. This is particularly effective when employing the abundance of H$_2$D$^+$. However, its low abundances and large distances make detections in massive sources hard to achieve. We present an application of the radiative transfer code POLARIS, with the goal to test the observability of the ortho-H$_2$D$^+$ transition $1_{10}$-$1_{11}$ (372.42 GHz) using simulations of high-mass collapsing cores that include deuteration chemistry. We analyzed an early and a late stage of the collapse of a 60 M$_{\odot}$ core, testing different source distances. For all cases, we generated synthetic single-dish and interferometric observations and studied the differences in both techniques. The column densities we derive are comparable to values reported for similar sources. These estimates depend on the extent over which they are averaged, and sources with compact emission they can be highly affected by beam dilution. Combined ALMA-ACA observations improve in signal-to-noise ratio and lead to better column density estimates as compared to ALMA alone. We confirm the feasibility to study ortho-H$_2$D$^+$ emission up to distances of 7 kpc. We provide a proof-of-concept of our framework for synthetic observations and highlight its importance when comparing numerical simulations with real observations. This work also proves how relevant it is to combine single-dish and interferometric measurements to derive appropriate source column densities.

preprint2022arXiv

The Twisted Magnetic Field of the Protobinary L483

We present H-band (1.65 $μ$m) and SOFIA HAWC+ 154 $μ$m polarization observations of the low-mass core L483. Our H-band observations reveal a magnetic field that is overwhelmingly in the E-W direction, which is approximately parallel to the bipolar outflow that is observed in scattered IR light and in single-dish $^{12}$CO observations. From our 154 $μ$m data, we infer a $\sim$ 45$^{\circ}$ twist in the magnetic field within the inner 5&#34; (1000 au) of L483. We compare these new observations with published single-dish 350 $μ$m polarimetry and find that the 10,000 au scale H-band data match the smaller scale 350 $μ$m data, indicating that the collapse of L483 is magnetically regulated on these larger scales. We also present high-resolution 1.3 mm ALMA data of L483 which reveals it is a close binary star with a separation of 34 au. The plane of the binary of L483 is observed to be approximately parallel to the twisted field in the inner 1000 au. Comparing this result to the $\sim$ 1000 au protostellar envelope, we find that the envelope is roughly perpendicular to the 1000 au HAWC+ field. Using the data presented, we speculate that L483 initially formed as a wide binary and the companion star migrated to its current position, causing an extreme shift in angular momentum thereby producing the twisted magnetic field morphology observed. More observations are needed to further test this scenario.

preprint2022arXiv

Topological superconductivity on the honeycomb lattice: Effect of normal state topology

The search for topological superconductors is one of the most pressing and challenging questions in condensed matter and material research. Despite some early suggestions that doping a topological insulator might be a successful recipe to find topological superconductors, until today there is no general understanding of the relationship of the topology of the superconductor and the topology of its underlying normal state system. One of the major obstacles is the strong effect of the Fermi surface and its subsequent pairing tendencies within the Hubbard model framework, usually preventing a detailed comparison between different topological superconducting systems. Here we present an analysis of doped insulators - topological and trivial - where the dominant Fermi surface effects have been equalized. Our approach allows us to study and compare superconducting instabilities of different normal state systems and present rigorous results about the influence of the normal state system&#39;s topology.

preprint2022arXiv

Triplet Superconductivity from Nonlocal Coulomb Repulsion in an Atomic Sn Layer Deposited onto a Si(111) Substrate

Atomic layers deposited on semiconductor substrates introduce a platform for the realization of the extended electronic Hubbard model, where the consideration of electronic repulsion beyond the onsite term is paramount. Recently, the onset of superconductivity at 4.7K has been reported in the hole-doped triangular lattice of tin atoms on a silicon substrate. Through renormalization group methods designed for weak and intermediate coupling, we investigate the nature of the superconducting instability in hole-doped Sn/Si(111). We find that the extended Hubbard nature of interactions is crucial to yield triplet pairing, which is f-wave (p-wave) for moderate (higher) hole doping. In light of persisting challenges to tailor triplet pairing in an electronic material, our finding promises to pave unprecedented ways for engineering unconventional triplet superconductivity.

preprint2021arXiv

Magnetic field structure of OMC-3 in the far infrared revealed by SOFIA/HAWC+

We report the SOFIA/HAWC+ band D (154$\,μ$m) and E (214$\,μ$m) polarimetric observations of the filamentary structure OMC-3 that is part of the Orion molecular cloud. The polarization pattern is uniform for both bands and parallel to the filament structure. The polarization degree decreases toward regions with high intensity for both bands, revealing a so called &#34;polarization hole.&#34; We identified an optical depth effect in which polarized emission and extinction act as counteracting mechanisms as a potential contributor to this phenomenon. Assuming that the detected polarization is caused by the emission of magnetically aligned non-spherical dust grains, the inferred magnetic field is uniform and oriented perpendicular to the filament. The magnetic field strength derived from the polarization patterns at 154$\,μ$m and 214$\,μ$m amounts to 202$\,μ$G and 261$\,μ$G, respectively. The derived magnetic field direction is consistent with that derived from previous polarimetric observations in the far infrared and submillimeter (submm) wavelength range. Investigating the far-infrared polarization spectrum derived from the SOFIA/HAWC+ observations, we do not find a clear correlation between the polarization spectrum and cloud properties, namely, the column density, $N(H_2$), and temperature, $T$.

preprint2021arXiv

Nonlinear interferometer for Fourier-transform mid-infrared gas spectroscopy using near-infrared detection

Nonlinear interferometers allow for mid-infrared spectroscopy with near-infrared detection using correlated photons. Previous implementations have demonstrated a spectral resolution limited by spectrally selective detection. In our work, we demonstrate mid-infrared transmission spectroscopy in a nonlinear interferometer using single-pixel near-infrared detection and Fourier-transform analysis. A sub-wavenumber spectral resolution allows for rotational-line-resolving spectroscopy of gaseous samples in a spectral bandwidth of over 700$\,$cm$^{-1}$. We use methane transmission spectra around 3.3$\,μ$m wavelength to characterize the spectral resolution, noise limitations and transmission accuracy of our device. The combination of nonlinear interferometry and Fourier-transform analysis paves the way towards performant and efficient mid-infrared spectroscopy with near-infrared detection.

preprint2021arXiv

The impact of the presence of water ice on the analysis of debris disk observations

The analysis of debris disk observations is often based on the assumption of a dust phase composed of compact spherical grains consisting of astronomical silicate. Instead, observations indicate the existence of water ice in debris disks. We quantify the impact of water ice as a potential grain constituent in debris disks on the disk parameter values estimated from photometric and spatially resolved observations in the mid- and far-infrared. We simulated photometric measurements and radial profiles of debris disks containing water ice and analyzed them by applying a disk model purely consisting of astronomical silicate. Subsequently, we quantified the deviations between the derived and the true parameter values. As stars in central positions we discuss a $β$ Pic sibling and main-sequence stars with spectral types ranging from A0 to K5. To simulate observable quantities we employed selected observational scenarios regarding the choice of wavelengths and instrument characteristics. For the $β$ Pic stellar model and ice fractions $\geq 50\ \%$ the derived inner disk radius is biased by ice sublimation toward higher values. However, the derived slope of the radial density profile is mostly unaffected. Along with an increasing ice fraction, the slope of the grain size distribution is overestimated by up to a median factor of $\sim 1.2$ for an ice fraction of $90\ \%$ while the total disk mass is underestimated by a factor of $\sim 0.4$. The reliability of the derived minimum grain size strongly depends on the spectral type of the central star. For an A0-type star the minimum grain size can be underestimated by a factor of $\sim 0.2$, while for solar-like stars it is overestimated by up to a factor of $\sim 4 - 5$. Neglecting radial profile measurements and using solely photometric measurements, the factor of overestimation increases for solar-like stars up to $\sim 7 - 14$.

preprint2021arXiv

Thermal imaging of dust hiding the black hole in the Active Galaxy NGC 1068

In the widely accepted &#39;Unified Model&#39; solution of the classification puzzle of Active Galactic Nuclei, the orientation of a dusty accretion torus around the central black hole dominates their appearance. In &#39;type-1&#39; systems, the bright nucleus is visible at the centre of a face-on torus. In &#39;type-2&#39; systems the thick, nearly edge-on torus hides the central engine. Later studies suggested evolutionary effects and added dusty clumps and polar winds but left the basic picture intact. However, recent high-resolution images of the archetypal type-2 galaxy NGC 1068 suggested a more radical revision. They displayed a ring-like emission feature which the authors advocated to be hot dust surrounding the black hole at the radius where the radiation from the central engine evaporates the dust. That ring is too thin and too far tilted from edge-on to hide the central engine, and ad hoc foreground extinction is needed to explain the type-2 classification. These images quickly generated reinterpretations of the type 1-2 dichotomy. Here we present new multi-band mid-infrared images of NGC1068 that detail the dust temperature distribution and reaffirm the original model. Combined with radio data, our maps locate the central engine below the previously reported ring and obscured by a thick, nearly edge-on disk, as predicted by the Unified Model. We also identify emission from polar flows and absorbing dust that is mineralogically distinct from that towards the Milky Way centre.

preprint2020arXiv

Constraining the magnetic field properties of Bok globule B335 using SOFIA/HAWC+

Thanks to their well-defined shape and mostly isolated locations, Bok globules are suitable objects for studying the physics of low-mass star formation. To study the magnetic field of the prototypical Bok globule B335, we obtained a spatially resolved polarization map with SOFIA/HAWC+ at a wavelength of 214$\,μ$m. For the first time, these observations reveal that polarization holes in Bok globules, that is, the decrease in polarization degree towards their dense centers, also occur in the far-infrared wavelength regime. The observed polarization pattern is uniform with a mean polarization angle of 48$^\circ\pm $26$^\circ$ and a magnetic field strength of $\sim$ 142$\,μ$G. Moreover, we use complementary polarimetic data for B335 obtained at near-infrared to millimeter wavelengths to analyze and constrain the magnetic field across different scales. By applying the 3D Monte-Carlo radiative transfer code POLARIS (Reissl et al. 2016), we developed a model for the density and magnetic field structure as well as for the dust properties of this globule. We conclude that the column density towards the center of B335 is too low to cause the observed polarization hole in B335 via dichroic absorption (Brauer et al. 2016). Furthermore, we conclude that the effect of self-scattering has no significant impact on the observed polarization. Adopting dust-grain alignment via the radiative torque mechanism, a combination of the interstellar radiation field and the central star as radiation sources is consistent with the decrease in polarization degree at the outer regions of B335 ($\approx\,$10$^4\,$au from the core). However, the model fails to explain the low polarization degree within the inner 5000 au.

preprint2020arXiv

Fourier transform infrared spectroscopy with visible light

Nonlinear interferometers allow spectroscopy in the mid-infrared range by detecting correlated visible light, for which non-cooled detectors with higher specific detectivity and lower dark count rates are available. We present a new approach for the registration of spectral information, which combines a nonlinear interferometer using non-degenerate spontaneous parametric down-conversion (SPDC) with a Fourier-transform spectroscopy concept. In order to increase the spectral coverage, we use broadband non-collinear SPDC in periodically poled LiNbO$_3$. Without the need for spectrally selective detection, continuous spectra with a spectral bandwidth of more than 100$\,$cm$^{-1}$ are achieved. We demonstrate transmission spectra of a polypropylene sample measured with ~ 6$\,$cm$^{-1}$ resolution in the spectral range between 3.2$\,μ$m to 3.9$\,μ$m.

preprint2020arXiv

Self-scattering in protoplanetary disks with dust settling

Scattering of re-emitted flux is considered to be at least partially responsible for the observed polarisation in the (sub-)millimetre wavelength range of several protoplanetary disks. Although the degree of polarisation produced by scattering is highly dependent on the dust model, early studies investigating this mechanism relied on the assumption of single grain sizes and simple density distribution of the dust. However, in the dense inner regions where this mechanism is usually most efficient, the existence of dust grains with sizes ranging from nanometres to millimetres has been confirmed. Additionally, the presence of gas forces larger grains to migrate vertically towards the disk midplane, introducing a dust segregation in the vertical direction. Using polarisation radiative transfer simulations, we analyse the dependence of the resulting scattered light polarisation at 350 $μ$m, 850 $μ$m, 1.3 mm, and 2 mm on various parameters describing protoplanetary disks, including the effect of dust grain settling. We find that the different disk parameters change the degree of polarisation mostly by affecting the anisotropy of the radiation field, the optical depth, or both. It is therefore very challenging to deduce certain disk parameter values directly from polarisation measurements alone. However, assuming a high dust albedo, it is possible to trace the transition from optically thick to optically thin disk regions. The degree of polarisation in most of the considered disk configurations is lower than what is found observationally, implying the necessity to revisit models that describe the dust properties and disk structure.

preprint2020arXiv

Synthetic observations of deuterated molecules in massive prestellar cores

Young massive stars are usually found embedded in dense massive molecular clumps and are known for being highly obscured and distant. During their formation process, deuteration is regarded as a potentially good indicator of the very early formation stages. In this work, we test the observability of the ground-state transition of ortho-H$_2$D$^+$ $J_{\rm {K_a, K_c}} = 1_{10}$-$1_{11} $ by performing interferometric and single-dish synthetic observations using magneto-hydrodynamic simulations of high-mass collapsing molecular cores, including deuteration chemistry. We studied different evolutionary times and source distances (from 1 to 7 kpc) to estimate the information loss when comparing the column densities inferred from the synthetic observations to the column densities in the model. We mimicked single-dish observations considering an APEX-like beam and interferometric observations using CASA and assuming the most compact configuration for the ALMA antennas. We found that, for centrally concentrated density distributions, the column densities are underestimated by about 51% in the case of high-resolution ALMA observations ($\leqslant$1&#34;) and up to 90% for APEX observations (17&#34;). Interferometers retrieve values closer to the real ones, however, their finite spatial sampling results in the loss of contribution from large-scale structures due to the lack of short baselines. We conclude that, the emission of o-H$_2$D$^+$ in distant massive dense cores is faint and would require from $\sim$1 to $\sim$7 hours of observation at distances of 1 and 7 kpc, respectively, to achieve a 14$σ$ detection in the best case scenario. Additionally, the column densities derived from such observations will certainly be affected by beam dilution in the case of single-dishes and spatial filtering in the case of interferometers.

preprint2020arXiv

Three-dimensional continuum radiative transfer of polarized radiation in exoplanetary atmospheres

Polarimetry is about to become a powerful tool for determining the atmospheric properties of exoplanets. To provide the basis for the interpretation of observational results and for predictive studies to guide future observations, sophisticated analysis tools are required. Our goal is to develop a radiative transfer tool that contains all the relevant continuum polarization mechanisms for the comprehensive analysis of the polarized flux resulting from the scattering in the atmosphere of, on the surface of, and in the local planetary environment (e.g., planetary rings, exomoons) of extra-solar planets. Furthermore, our goal is to avoid common simplifications such as locally plane-parallel planetary atmospheres, the missing cross-talk between latitudinal and longitudinal regions, or the assumption of either a point-like star or plane-parallel illumination. As a platform for the newly developed numerical algorithms, we use the 3D Monte Carlo radiative transfer code POLARIS. The code is extended and optimized for the radiative transfer in exoplanetary atmospheres. We investigate the reflected flux and its degree of polarization for different phase angles for a homogeneous cloud-free atmosphere and an inhomogeneous cloudy atmosphere. To take advantage of the 3D radiative transfer and to demonstrate the potential of the code, the impact of an additional circumplanetary ring on the reflected polarized flux is studied. The presence of a circumplanetary ring consisting of small water-ice particles has a noticeable impact on the reflected polarized radiation. In particular, the reflected flux strongly increases at larger phase angles if the planetary orbit is seen edge-on because the considered particles tend to scatter forwards. In contrast, the degree of polarization decreases at these phase angles.

preprint2019arXiv

Characterization of mid-infrared polarization due to scattering in protoplanetary disks

It is generally assumed that magnetic fields play an important role in the formation and evolution of protoplanetary disks. One way of observationally constraining magnetic fields is to measure polarized emission and absorption produced by magnetically aligned elongated dust grains. The fact that radiation also becomes linearly polarized by light scattering at optical to millimeter wavelengths complicates magnetic field studies. We characterize the linear polarization of mid-infrared radiation due to scattering of the stellar radiation and dust thermal re-emission radiation (self-scattering). We find that the thermal re-emission radiation is stronger than the scattered stellar radiation for disks with inner holes smaller than 10 au within the considered parameter range. The mid-infrared polarization due to scattering shows several clear trends: For scattered stellar radiation only, the linear polarization degree decreases slightly with increasing radial distance, while it increases with radial distance for thermal re-emission radiation only and for a combination of scattered stellar radiation and thermal re-emission radiation. The linear polarization degree decreases with increasing disk flaring and luminosity of the central star. An increasing inner radius shifts the increase of the linear polarization degree further outside, while a larger scale height increases the linear polarization degree for small radial distances and decreases this degree further outside. For longer wavelengths, the linear polarization degree converges more slowly.

preprint2019arXiv

Laser modulation of superconductivity in a cryogenic widefield nitrogen-vacancy microscope

Microscopic imaging based on nitrogen-vacancy (NV) centres in diamond, a tool increasingly used for room-temperature studies of condensed matter systems, has recently been extended to cryogenic conditions. However, it remains unclear whether the technique is viable for imaging temperature-sensitive phenomena below 10 K given the inherent laser illumination requirements, especially in a widefield configuration. Here we realise a widefield NV microscope with a field of view of 100 $μ$m and a base temperature of 4 K, and use it to image Abrikosov vortices and transport currents in a superconducting Nb film. We observe the disappearance of vortices upon increase of laser power and their clustering about hot spots upon decrease, indicating that laser powers as low as 1 mW (4 orders of magnitude below the NV saturation) are sufficient to locally quench the superconductivity of the film ($T_c = 9$ K). This significant local heating is confirmed by resistance measurements, which reveal the presence of large temperature gradients (several K) across the film. We then investigate the effect of such gradients on transport currents, where the current path is seen to correlate with the temperature profile even in the fully superconducting phase. In addition to highlighting the role of temperature inhomogeneities in superconductivity phenomena, this work establishes that, under sufficiently low laser power conditions, widefield NV microscopy enables imaging over mesoscopic scales down to 4 K with a submicrometer spatial resolution, providing a new platform for real-space investigations of a range of systems from topological insulators to van der Waals ferromagnets.

preprint2019arXiv

Light from an ion crystal: bunching or antibunching?

Photon statistics divides light sources into three different categories, characterized by bunched, antibunched or uncorrelated photon arrival times. Single atoms, ions, molecules, or solid state emitters display antibunching of photons, while classical thermal sources exhibit photon bunching. Here we demonstrate a light source in free space, where the photon statistics depends on the direction of observation, undergoing a continuous crossover between photon bunching and antibunching. We employ two trapped ions, observe their fluorescence under continuous laser light excitation, and record the spatially resolved autocorrelation function $g^{(2)}(τ)$ with a movable Hanbury Brown and Twiss detector. Varying the detector position we find a minimum value for antibunching, $g^{(2)}(0) = 0.60(5)$ and a maximum of $g^{(2)}(0)=1.46(8)$ for bunching, demonstrating that this source radiates fundamentally different types of light alike. The observed variation of the autocorrelation function is understood in the Dicke model of heralded entangled states and the observed maximum and minimum values are modeled, taking independently measured experimental parameters into account.

preprint2018arXiv

Variational Monte Carlo - Bridging Concepts of Machine Learning and High Dimensional Partial Differential Equations

A statistical learning approach for parametric PDEs related to Uncertainty Quantification is derived. The method is based on the minimization of an empirical risk on a selected model class and it is shown to be applicable to a broad range of problems. A general unified convergence analysis is derived, which takes into account the approximation and the statistical errors. By this, a combination of theoretical results from numerical analysis and statistics is obtained. Numerical experiments illustrate the performance of the method with the model class of hierarchical tensors.