Catalog footprint

What is connected

66works
22topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

66 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's periphery, then we can place an upper limit of $<27\,μ$G on the magnetic field strength in the globule'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" (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'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 "polarization hole." 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 'Unified Model' 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 'type-1' systems, the bright nucleus is visible at the centre of a face-on torus. In 'type-2' 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") and up to 90% for APEX observations (17"). 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.

preprint2016arXiv

Magnetic fields in protoplanetary disks: from MHD simulations to ALMA observations

Magnetic fields significantly influence the evolution of protoplanetary disks and the formation of planets, following the predictions of numerous magnetohydrodynamic (MHD) simulations. However, these predictions are yet observationally unconstrained. To validate the predictions on the influence of magnetic fields on protoplanetary disks, we apply 3D radiative transfer simulations of the polarized emission of aligned aspherical dust grains that directly link 3D global non-ideal MHD simulations to ALMA observations. Our simulations show that it is feasible to observe the predicted toroidal large-scale magnetic field structures, not only in the ideal observations but also with high-angular resolution ALMA observations. Our results show further that high angular resolution observations by ALMA are able to identify vortices embedded in outer magnetized disk regions.

preprint2016arXiv

Origin and evolution of two-component debris discs and an application to the q$^1$ Eridani system

Many debris discs reveal a two-component structure, with an outer Kuiper-belt analogue and a warm inner component whose origin is still a matter of debate. One possibility is that warm emission stems from an "asteroid belt" closer in to the star. We consider a scenario in which a set of giant planets is formed in an initially extended planetesimal disc. These planets carve a broad gap around their orbits, splitting up the disc into the outer and the inner belts. After the gas dispersal, both belts undergo collisional evolution in a steady-state regime. This scenario is explored with detailed collisional simulations involving realistic physics to describe a long-term collisional depletion of the two-component disc. We find that the inner disc may be able to retain larger amounts of material at older ages than thought before on the basis of simplified analytic models. We show that the proposed scenario is consistent with a suite of thermal emission and scattered light observational data for a bright two-temperature debris disc around a nearby solar-type star q$^1$ Eridani. This implies a Solar System-like architecture of the system, with an outer massive "Kuiper belt", an inner "asteroid belt", and a few Neptune- to Jupiter-mass planets in between.

preprint2016arXiv

Planet Formation Imager (PFI): science vision and key requirements

The Planet Formation Imager (PFI) project aims to provide a strong scientific vision for ground-based optical astronomy beyond the upcoming generation of Extremely Large Telescopes. We make the case that a breakthrough in angular resolution imaging capabilities is required in order to unravel the processes involved in planet formation. PFI will be optimised to provide a complete census of the protoplanet population at all stellocentric radii and over the age range from 0.1 to about 100 Myr. Within this age period, planetary systems undergo dramatic changes and the final architecture of planetary systems is determined. Our goal is to study the planetary birth on the natural spatial scale where the material is assembled, which is the "Hill Sphere" of the forming planet, and to characterise the protoplanetary cores by measuring their masses and physical properties. Our science working group has investigated the observational characteristics of these young protoplanets as well as the migration mechanisms that might alter the system architecture. We simulated the imprints that the planets leave in the disk and study how PFI could revolutionise areas ranging from exoplanet to extragalactic science. In this contribution we outline the key science drivers of PFI and discuss the requirements that will guide the technology choices, the site selection, and potential science/technology tradeoffs.

preprint2016arXiv

Radiative transfer with POLARIS: I. Analysis of magnetic fields through synthetic dust continuum polarization measurements

Aims: We present POLARIS (POLArized RadIation Simulator), a newly developed three-dimensional Monte-Carlo radiative transfer code. POLARIS was designed to calculate dust temperature, polarization maps, and spectral energy distributions. It is optimized to handle data that results from sophisticated magneto-hydrodynamic simulations. The main purpose of the code is to prepare and analyze multi-wavelength continuum polarization measurements in the context of magnetic field studies in the interstellar medium. An exemplary application is the investigation of the role of magnetic fields in star formation processes. Methods: We combine currently discussed state-of-the-art grain alignment theories with existing dust heating and polarization algorithms. We test the POLARIS code on multiple scales in complex astrophysical systems that are associated with different stages of star formation. POLARIS uses the full spectrum of dust polarization mechanisms to trace the underlying magnetic field morphology. Results: Resulting temperature distributions are consistent with the density and position of radiation sources resulting from magneto-hydrodynamic (MHD) - collapse simulations. The calculated layers of aligned dust grains in the considered cirumstellar disk models are in excellent agreement with theoretical predictions. Finally, we compute unique patterns in synthetic multi-wavelength polarization maps that are dependent on applied dust-model and grain-alignment theory in analytical cloud models.

preprint2016arXiv

RetiNet: Automatic AMD identification in OCT volumetric data

Optical Coherence Tomography (OCT) provides a unique ability to image the eye retina in 3D at micrometer resolution and gives ophthalmologist the ability to visualize retinal diseases such as Age-Related Macular Degeneration (AMD). While visual inspection of OCT volumes remains the main method for AMD identification, doing so is time consuming as each cross-section within the volume must be inspected individually by the clinician. In much the same way, acquiring ground truth information for each cross-section is expensive and time consuming. This fact heavily limits the ability to acquire large amounts of ground truth, which subsequently impacts the performance of learning-based methods geared at automatic pathology identification. To avoid this burden, we propose a novel strategy for automatic analysis of OCT volumes where only volume labels are needed. That is, we train a classifier in a semi-supervised manner to conduct this task. Our approach uses a novel Convolutional Neural Network (CNN) architecture, that only needs volume-level labels to be trained to automatically asses whether an OCT volume is healthy or contains AMD. Our architecture involves first learning a cross-section pathology classifier using pseudo-labels that could be corrupted and then leverage these towards a more accurate volume-level classification. We then show that our approach provides excellent performances on a publicly available dataset and outperforms a number of existing automatic techniques.

preprint2016arXiv

The circumstellar disk of FS Tau B - A self-consistent model based on observations in the mid-infrared with NACO -

Protoplanetary disks are a byproduct of the star formation process. In the dense mid-plane of these disks, planetesimals and planets are expected to form. The first step in planet formation is the growth of dust particles from submicrometer-sized grains to macroscopic mm-sized aggregates. The grain growth is accompanied by radial drift and vertical segregation of the particles within the disk. To understand this essential evolutionary step, spatially resolved multi-wavelength observations as well as photometric data are necessary which reflect the properties of both disk and dust. We present the first spatially resolved image obtained with NACO at the VLT in the L$_\text{p}$ band of the near edge-on protoplanetary disk FS Tau B. Based on this new image, a previously published Hubble image in H band and the spectral energy distribution from optical to millimeter wavelengths, we derive constraints on the spatial dust distribution and the progress of grain growth. For this purpose we perform a disk modeling using the radiative transfer code MC3D. Radial drift and vertical sedimentation of the dust are not considered. We find a best-fit model which features a disk extending from $2\,\text{AU}$ to several hundreds AU with a moderately decreasing surface density and $M_\text{disk}=2.8\,\times\,10^{-2}\,\text{M}_\odot$. The inclination amounts to $i=80^\circ$. Our findings indicate that substantial dust grain growth has taken place and that grains of a size equal to or larger than $1\,\text{mm}$ are present in the disk. In conclusion, the parameters describing the vertical density distribution are better constrained than those describing the radial disk structure.

preprint2016arXiv

The influence of dust grain porosity on the analysis of debris disc observations

Debris discs are often modelled assuming compact dust grains, but more and more evidence for the presence of porous grains is found. We aim at quantifying the systematic errors introduced when modelling debris discs composed of porous dust with a disc model assuming spherical, compact grains. We calculate the optical dust properties derived via the fast, but simple effective medium theory. The theoretical lower boundary of the size distribution -- the so-called 'blowout size' -- is compared in the cases of compact and porous grains. Finally, we simulate observations of hypothetical debris discs with different porosities and feed them into a fitting procedure using only compact grains. The deviations of the results for compact grains from the original model based on porous grains are analysed. We find that the blowout size increases with increasing grain porosity up to a factor of two. An analytical approximation function for the blowout size as a function of porosity and stellar luminosity is derived. The analysis of the geometrical disc set-up, when constrained by radial profiles, are barely affected by the porosity. However, the determined minimum grain size and the slope of the grain size distribution derived using compact grains are significantly overestimated. Thus, the unexpectedly high ratio of minimum grain size to blowout size found by previous studies using compact grains can be partially described by dust grain porosity, although the effect is not strong enough to completely explain the trend.

preprint2016arXiv

Visibility of Young's interference fringes: Scattered light from small ion crystals

We observe interference in the light scattered from trapped $^{40}$Ca$^+$ ion crystals. By varying the intensity of the excitation laser, we study the influence of elastic and inelastic scattering on the visibility of the fringe pattern and discriminate its effect from that of the ion temperature and wave-packet localization. In this way we determine the complex degree of coherence and the mutual coherence of light fields produced by individual atoms. We obtain interference fringes from crystals consisting of two, three and four ions in a harmonic trap. Control of the trapping potential allows for the adjustment of the interatomic distances and thus the formation of linear arrays of atoms serving as a regular grating of microscopic scatterers.

preprint2015arXiv

Microscopy with a Deterministic Single Ion Source

We realize a single particle microscope by using deterministically extracted laser cooled $^{40}$Ca$^+$ ions from a Paul trap as probe particles for transmission imaging. We demonstrate focusing of the ions with a resolution of 5.8$\;\pm\;$1.0$\,$nm and a minimum two-sample deviation of the beam position of 1.5$\,$nm in the focal plane. The deterministic source, even when used in combination with an imperfect detector, gives rise to much higher signal to noise ratios as compared with conventional Poissonian sources. Gating of the detector signal by the extraction event suppresses dark counts by 6 orders of magnitude. We implement a Bayes experimental design approach to microscopy in order to maximize the gain in spatial information. We demonstrate this method by determining the position of a 1$\,μ$m circular hole structure to an accuracy of 2.7$\,$nm using only 579 probe particles.

preprint2015arXiv

Structures in circumbinary disks: Prospects for observability

During the past decade circumbinary disks have been discovered around various young binary stars. Hydrodynamical calculations indicate that the gravitational interaction between the central binary star and the surrounding disk results in global perturbations of the disk density profile. We study the observability of characteristic large-scale disk structures resulting from the binary-disk interaction in the case of close binary systems. We derived the structure of circumbinary disks from smoothed-particle hydrodynamic simulations. Subsequently, we performed radiative transfer simulations to obtain scattered light and thermal reemission maps. We investigated the influence of the binary mass ratio, the inclination of the binary orbit relative to the disk midplane, and the eccentricity of the binary orbit on observational quantities. We find that ALMA will allow tracing asymmetries of the inner edge of the disk and potentially resolving spiral arms if the disk is seen face-on. For an edge-on orientation, ALMA will allow detecting perturbations in the disk density distribution through asymmetries in the radial brightness profile. Through the asymmetric structure of the disks, areas are formed with a temperature $2.6$ times higher than at the same location in equivalent unperturbed disks. The time-dependent appearance of the density waves and spiral arms in the disk affects the total re-emission flux of the object by a few percent.

preprint2014arXiv

HD139614: the interferometric case for a group-Ib pre-transitional young disk

The Herbig Ae star HD 139614 is a group-Ib object, which featureless SED indicates disk flaring and a possible pre-transitional evolutionary stage. We present mid- and near-IR interferometric results collected with MIDI, AMBER and PIONIER with the aim of constraining the spatial structure of the 0.1-10 AU disk region and assess its possible multi-component structure. A two-component disk model composed of an optically thin 2-AU wide inner disk and an outer temperature-gradient disk starting at 5.6 AU reproduces well the observations. This is an additional argument to the idea that group-I HAeBe inner disks could be already in the disk-clearing transient stage. HD 139614 will become a prime target for mid-IR interferometric imaging with the second-generation instrument MATISSE of the VLTI.

preprint2014arXiv

Large-scale magnetic fields in Bok globules

Context: The role of magnetic fields in the star formation process is a contentious matter of debate. In particular, no clear observational proof exists of a general influence by magnetic fields during the initial collapse of molecular clouds. Aims: Our aim is to examine magnetic fields and their influence on a wide range of spatial scales in low-mass star-forming regions. Method: We trace the large-scale magnetic field structure on scales of 10^3-10^5 AU in the local environment of Bok globules through optical and near-infrared polarimetry and combine these measurements with existing submillimeter measurements, thereby characterizing the small-scale magnetic field structure on scales of 10^2-10^3 AU. Results: For the first time, we present polarimetric observations in the optical and near-infrared of the three Bok globules B335, CB68, and CB54, combined with archival observations in the submillimeter and the optical. We find a significant polarization signal (P>=2%, P/sigma(P)>3) in the optical and near-infrared for all three globules. Additionally, we detect a connection between the structure on scales of 10^2-10^3 AU to 10^3-10^4 AU for both B335 and CB68. Furthermore, for CB54, we trace ordered polarization vectors on scales of ~10^5 AU. We determine a magnetic field orientation that is aligned with the CO outflow in the case of CB54, but nearly perpendicular to the CO outflow for CB68. For B335 we find a change in the magnetic field oriented toward the outflow direction, from the inner core to the outer regions. Conclusion: We find strongly aligned polarization vectors that indicate dominant magnetic fields on a wide range of spatial scales.

preprint2014arXiv

Planet-induced disk structures: A comparison between (sub)mm and infrared radiation

Young giant planets, which are embedded in a circumstellar disk, will significantly perturb the disk density distribution. This effect can potentially be used as an indirect tracer for planets. We investigate the feasibility of observing planet-induced gaps in circumstellar disks in scattered light. We perform 3D hydrodynamical disk simulations combined with subsequent radiative transfer calculations in scattered light for different star, disk, and planet configurations. The results are compared to those of a corresponding study for the (sub)mm thermal re-emission. The feasibility of detecting planet-induced gaps in scattered light is mainly influenced by the optical depth of the disk and therefore by the disk size and mass. Planet-induced gaps are in general only detectable if the photosphere of the disks is sufficiently disturbed. Within the limitations given by the parameter space here considered, we find that gap detection is possible in the case of disks with masses below $\sim 10^{-4\dots-3} \, \rm M_\odot$. Compared to the disk mass that marks the lower Atacama Large (Sub)Millimeter Array (ALMA) detection limit for the thermal radiation re-emitted by the disk, it is possible to detect the same gap both in re-emission and scattered light only in a narrow range of disk masses around $\sim 10^{-4} \, \rm M_\odot$, corresponding to $16\%$ of cases considered in our study.

preprint2014arXiv

Shadows and cavities in protoplanetary disks: HD163296, HD141569A, and HD150193A in polarized light

The morphological evolution of dusty disks around young (few Myr-old) stars is pivotal to better understand planet formation. Since both dust grains and the global disk geometry evolve on short timescale, high-resolution imaging of a sample of objects may provide important hints towards such an evolution. We enlarge the sample of protoplanetary disks imaged in polarized light with high-resolution by observing the Herbig Ae/Be stars HD163296, HD141569A, and HD150193A. We integrate our data with previous datasets to paint a larger picture of their morphology. We report a weak detection of the disk around HD163296 in both H and Ks band. The disk is resolved as a broken ring structure with a significan surface brightness drop inward of 0.6 arcsec. No sign of extended polarized emission is detected from the disk around HD141569A and HD150193A. We propose that the absence of scattered light in the inner 0.6 arcsec around HD163296 and the non-detection of the disk around HD150193A may be due to similar geometric factors. Since these disks are known to be flat or only moderately flared, self-shadowing by the disk inner wall is the favored explanation. We show that the polarized brightness of a number of disks is indeed related to their flaring angle. Other scenarios (such as dust grain growth or interaction with icy molecules) are also discussed. On the other hand, the non-detection of HD141569A is consistent with previous datasets revealing the presence of a huge cavity in the dusty disk.

preprint2014arXiv

Single particle microscopy with nanometer resolution

We experimentally demonstrate nanoscopic transmission microscopy relying on a deterministic single particle source. This increases the signal-to-noise ratio with respect to conventional microscopy methods, which employ Poissonian particle sources. We use laser-cooled ions extracted from a Paul trap, and demonstrate remote imaging of transmissive objects with a resolution of 8.6 $\pm$ 2.0nm and a minimum two-sample deviation of the beam position of 1.5nm. Detector dark counts can be suppressed by 6 orders of magnitudes through gating by the extraction event. The deterministic nature of our source enables an information-gain driven approach to imaging. We demonstrate this by performing efficient beam characterization based on a Bayes experiment design method.

preprint2013arXiv

A diversity of dusty AGN tori: Data release for the VLTI/MIDI AGN Large Program and first results for 23 galaxies

The AGN-heated dust distribution (the "torus") is increasingly recognized not only as the absorber required in unifying models, but as a tracer for the reservoir that feeds the nuclear Super-Massive Black Hole. Yet, even its most basic structural properties (such as its extent, geometry and elongation) are unknown for all but a few archetypal objects. Since most AGNs are unresolved in the mid-infrared, we utilize the MID-infrared interferometric Instrument (MIDI) at the Very Large Telescope Interferometer (VLTI) that is sensitive to structures as small as a few milli-arcseconds (mas). We present here an extensive amount of new interferometric observations from the MIDI AGN Large Program (2009 - 2011) and add data from the archive to give a complete view of the existing MIDI observations of AGNs. Additionally, we have obtained high-quality mid-infrared spectra from VLT/VISIR. We present correlated and total flux spectra for 23 AGNs and derive flux and size estimates at 12 micron using simple axisymmetric geometrical models. Perhaps the most surprising result is the relatively high level of unresolved flux and its large scatter: The median "point source fraction" is 70 % for type 1 and 47 % for type 2 AGNs meaning that a large part of the flux is concentrated on scales smaller than about 5 mas (0.1 - 10 pc). Among sources observed with similar spatial resolution, it varies from 20 % - 100 %. For 18 of the sources, two nuclear components can be distinguished in the radial fits. While these models provide good fits to all but the brightest sources, significant elongations are detected in eight sources. The half-light radii of the fainter sources are smaller than expected from the size ~ L^0.5 scaling of the bright sources and show a large scatter, especially when compared to the relatively tight size--luminosity relation in the near-infrared.

preprint2013arXiv

Evidence of a discontinuous disk structure around the Herbig Ae star HD 139 614

A new class of pre-main sequence objects has been recently identified as pre-transitional disks. They present near-infrared excess coupled to a flux deficit at about 10 microns and a rising mid-infrared and far-infrared spectrum. These features suggest a disk structure with inner and outer dust components, separated by a dust-depleted region (or gap). We here report on the first interferometric observations of the disk around the Herbig Ae star HD 139614. Its infrared spectrum suggests a flared disk, and presents pre-transitional features,namely a substantial near-infrared excess accompanied by a dip around 6 microns and a rising mid-infrared part. In this framework, we performed a study of the spectral energy distribution (SED) and the mid-infrared VLTI/MIDI interferometric data to constrain thespatial structure of the inner dust disk region and assess its possibly multi-component structure. We based our work on a temperature-gradient disk model that includes dust opacity. While we could not reproduce the SED and interferometric visibilities with a one-component disk, a better agreement was obtained with a two-component disk model composed of an optically thin inner disk extending from 0.22 to 2.3 au, a gap, and an outer temperature-gradient disk starting at 5.6 au. Therefore, our modeling favors an extended and optically thin inner dust component and in principle rules out the possibility that the near-infrared excess originates only from a spatially confined region. Moreover, the outer disk is characterized by a very steep temperature profile and a temperature higher than 300 K at its inner edge. This suggests the existence of a warm component corresponding to a scenario where the inner edge of the outer disk is directly illuminated by the central star. This is an expected consequence of the presence of a gap, thus indicative of a pre-transitional structure.

preprint2013arXiv

Gaps in the HD169142 protoplanetary disk revealed by polarimetric imaging: Signs of ongoing planet formation?

We present H-band VLT/NACO polarized light images of the Herbig Ae/Be star HD169142 probing its protoplanetary disk as close as ~0.1" to the star. Our images trace the face-on disk out to ~1.7" (~250 AU) and reveal distinct sub-structures for the first time: 1) the inner disk (<20 AU) appears to be depleted in scattering dust grains; 2) an unresolved disk rim is imaged at ~25 AU; 3) an annular gap extends from ~40 - 70 AU; 4) local brightness asymmetries are found on opposite sides of the annular gap. We discuss different explanations for the observed morphology among which ongoing planet formation is a tempting - but yet to be proven - one. Outside of ~85 AU the surface brightness drops off roughly r^{-3.3}, but describing the disk regions between 85-120 AU / 120-250 AU separately with power-laws r^{-2.6} / r^{-3.9} provides a better fit hinting towards another discontinuity in the disk surface. The flux ratio between the disk integrated polarized light and the central star is ~4.1 * 10^{-3}. Finally, combining our results with those from the literature, ~40% of the scattered light in the H-band appears to be polarized. Our results emphasize that HD169142 is an interesting system for future planet formation or disk evolution studies.

preprint2013arXiv

Mid-infrared observations of the circumstellar disks around PDS 66 and CRBR 2422.8-3423

Aims: We present mid-infrared observations and photometry of the circumstellar disks around PDS 66 and CRBR 2422.8-3423, obtained with VISIR/VLT in the N band and for the latter also in the Q band. Our aim is to resolve the inner regions of these protoplanetary disks, which carry potential signatures of intermediate or later stages of disk evolution and ongoing planet formation. Methods: We determined the radial brightness profiles of our target objects and the corresponding PSF reference that were observed before and after our target objects. Background standard deviations, the standard errors, and the seeing variations during the observations were considered. Adopting a simple radiative transfer model based on parameters taken from previous studies, we derived constraints on the inner-disk hole radius of the dust disk. Results: Neither of the circumstellar disks around our science targets are spatially resolved in our observations. However, we are able to constrain the inner-disk hole radius to <15.0(+0.5,-0.5) AU and <10.5(+0.5,-1.0) AU for PDS 66 and CRBR 2422.8-3423, respectively. The photometry we performed yields N-band flux densities of 599+-8 mJy for PDS 66 and 130+-14 mJy for CRBR 2422.8-3423, as well as a Q-band flux density of 858+-109 mJy for CRBR 2422.8-3423.

preprint2013arXiv

Porous dust grains in debris disks

When modeling the density and grain size distribution in debris disks, the minimum particle size is often significantly larger than the corresponding blowout size. While the dust particles are usually modeled as compact, homogenous spheres, we instead investigate the impact of porosity. The optical properties of porous particles are determined, and the influences of porosity on the blowout size and dust temperatures investigated. Using the method of discrete dipole approximation, we calculate the scattering and absorption cross sections of porous particles and derive the blowout size and the behavior of the dust temperature. We investigate the influence on the beta-ratio. Blowout sizes are calculated for various stellar luminosities and porosities, and an approximation equation is derived to estimate the blowout size as a function of these parameters. Furthermore, we investigate the influence of the porosity on the dust equilibrium temperature. The blowout size increases with the particle porosity and stellar luminosity. The dust temperature of porous particles is lower than the one of the compact spheres, in particular the temperature of blowout grains decreases for porous particles.

preprint2013arXiv

Small vs large dust grains in transitional disks: do different cavity sizes indicate a planet?

Transitional disks represent a short stage of the evolution of circumstellar material. Studies of dust grains in these objects can provide pivotal information on the mechanisms of planet formation. Dissimilarities in the spatial distribution of small (micron-size) and large (millimeter-size) dust grains have recently been pointed out. Constraints on the small dust grains can be obtained by imaging the distribution of scattered light at near-infrared wavelengths. We aim at resolving structures in the surface layer of transitional disks (with particular emphasis on the inner 10 - 50 AU), thus increasing the scarce sample of high resolution images of these objects. We obtained VLT/NACO near-IR high-resolution polarimetric differential imaging observations of SAO 206462 (HD135344B). This technique allows one to image the polarized scattered light from the disk without any occulting mask and to reach an inner working angle of 0.1''. A face-on disk is detected in H and Ks bands between 0.1'' and 0.9''. No significant differences are seen between the H and Ks images. In addition to the spiral arms, these new data allow us to resolve for the first time an inner cavity for small dust grains. The cavity size (about 28 AU) is much smaller than what is inferred for large dust grains from (sub)mm observations (39 to 50 AU). The interaction between the disk and potential orbiting companion(s) can explain both the spiral arm structure and the discrepant cavity sizes for small and large dust grains. One planet may be carving out the gas (and, thus, the small grains) at 28 AU, and generating a pressure bump at larger radii (39 AU), which holds back the large grains. We analytically estimate that, in this scenario, a single giant planet (with a mass between 5 and 15 Jupiter masses) at 17 to 20 AU from the star is consistent with the observed cavity sizes.

preprint2013arXiv

Structures in the protoplanetary disk of HD142527 seen in polarized scattered light

We present H- and Ks-band polarized differential images (PDI) of the Herbig Ae/Be star HD142527, revealing its optically thick outer disk and the nearly empty gap. The very small inner working angle (~0.1") and high resolution achievable with an 8m-class telescope, together with a careful polarimetric calibration strategy, allow us to achieve images that surpass the quality of previous scattered light images. Previously known substructures are resolved more clearly and new structures are seen. Specifically, we are able to resolve 1) half a dozen spiral structures in the disk, including previously known outer-disk spirals as well as new spiral arms and arcs close to the inner rim of the disk; 2) peculiar holes in the polarized surface brightness at position angles of ~0° and ~160°; 3) the inner rim on the eastern side of the disk; 4) the gap between the outer and inner disk, ranging from the inner working angle of 0.1" out to between 0.7 and 1.0", which is nearly devoid of dust. We then use a Markov-chain Monte-Carlo algorithm to determine several structural parameters of the disk, using very simple assumptions, including its inclination, eccentricity, and the scale height of the inner rim. We compare our results to previous work on this object, and try to produce a consistent picture of the system and its transition disk.

preprint2013arXiv

Vertical settling and radial segregation of large dust grains in the circumstellar disk of the Butterfly Star

Context: Circumstellar disks are considered to be the environment for the formation of planets. The growth of dust grains in these disks is the first step in the core accretion-gas capture planet formation scenario. Indicators and evidence of disk evolution can be traced in spatially resolved images and the spectral energy distribution (SED) of these objects. Aims: We develop a model for the dust phase of the edge-on oriented circumstellar disk of the Butterfly Star which allows one to fit observed multi-wavelength images and the SED simultaneously. Methods: Our model is based on spatially resolved high angular resolution observations at 1.3 mm, 894 micron, 2.07 micron, 1.87 micron, 1.60 micron, and 1.13 micron and an extensively covered SED ranging from 12 micron to 2.7 mm, including a detailed spectrum obtained with the Spitzer Space Telescope in the range from 12 micron to 38 micron. A parameter study based on a grid search method involving the detailed analysis of every parameter was performed to constrain the disk parameters and find the best-fit model for the independent observations. The individual observations were modeled simultaneously, using our continuum radiative transfer code. Results: We derived a model that is capable of reproducing all of the observations of the disk at the same time. We find quantitative evidence for grain growth up to ~100 micron-sized particles, vertical settling of larger dust grains toward the disk midplane, and radial segregation of the latter toward the central star. Within our best-fit model the large grains have a distribution with a scale height of 3.7 AU at 100 AU and a radial extent of 175 AU compared to a hydrostatic scale height of 6.9 AU at 100 AU and an outer disk radius of 300 AU. Our results are consistent with current theoretical models for the evolution of circumstellar disks and the early stages of planet formation.

preprint2012arXiv

A Comparison of Approaches in Fitting Continuum SEDs

We present a detailed comparison of two approaches, the use of a pre-calculated database and simulated annealing (SA), for fitting the continuum spectral energy distribution (SED) of astrophysical objects whose appearance is dominated by surrounding dust. While pre-calculated databases are commonly used to model SED data, only few studies to date employed SA due to its unclear accuracy and convergence time for this specific problem. From a methodological point of view, different approaches lead to different fitting quality, demand on computational resources and calculation time. We compare the fitting quality and computational costs of these two approaches for the task of SED fitting to provide a guide to the practitioner to find a compromise between desired accuracy and available resources. To reduce uncertainties inherent to real datasets, we introduce a reference model resembling a typical circumstellar system with 10 free parameters. We derive the SED of the reference model with our code MC3D at 78 logarithmically distributed wavelengths in the range [0.3um, 1.3mm] and use this setup to simulate SEDs for the database and SA. Our result shows directly the applicability of SA in the field of SED modeling, since the algorithm regularly finds better solutions to the optimization problem than a pre-calculated database. As both methods have advantages and shortcomings, a hybrid approach is preferable. While the database provides an approximate fit and overall probability distributions for all parameters deduced using Bayesian analysis, SA can be used to improve upon the results returned by the model grid.

preprint2012arXiv

A Herschel Survey of Cold Dust in Disks Around Brown Dwarfs and Low-Mass Stars

We report the complete photometric results from our Herschel study which is the first comprehensive program to search for far-infrared emission from cold dust around young brown dwarfs. We surveyed 50 fields containing 51 known or suspected brown dwarfs and very low mass stars that have evidence of circumstellar disks based on Spitzer photometry and/or spectroscopy. The objects with known spectral types range from M3 to M9.5. Four of the candidates were subsequently identified as extragalactic objects. Of the remaining 47 we have successfully detected 36 at 70micron and 14 at 160micron with S/N greater than 3, as well as several additional possible detections with low S/N. The objects exhibit a range of [24]--[70] micron colors suggesting a range in mass and/or structure of the outer disk. We present modeling of the spectral energy distributions of the sample and discuss trends visible in the data. Using two Monte Carlo radiative transfer codes we investigate disk masses and geometry. We find a very wide range in modeled total disk masses from less than 1e-6 solar masses up to 1e-3 solar masses with a median disk mass of order 3e-5 solar masses, suggesting that the median ratio of disk mass to central object mass may be lower than for T Tauri stars. The disk scale heights and flaring angles, however, cover a range consistent with those seen around T Tauri stars. The host clouds in which the young brown dwarfs and low-mass stars are located span a range in estimated age from ~1-3 Myr to ~10 Myr and represent a variety of star-forming environments. No obvious dependence on cloud location or age is seen in the disk properties, though the statistical significance of this conclusion is not strong.

preprint2012arXiv

Circumstellar disks and planets. Science cases for next-generation optical/infrared long-baseline interferometers

We present a review of the interplay between the evolution of circumstellar disks and the formation of planets, both from the perspective of theoretical models and dedicated observations. Based on this, we identify and discuss fundamental questions concerning the formation and evolution of circumstellar disks and planets which can be addressed in the near future with optical and infrared long-baseline interferometers. Furthermore, the importance of complementary observations with long-baseline (sub)millimeter interferometers and high-sensitivity infrared observatories is outlined.

preprint2012arXiv

Observing planet-disk interaction in debris disks

Context. Structures in debris disks induced by planetdisk interaction are promising to provide valuable constraints on the existence and properties of embedded planets. Aims. We investigate the observability of structures in debris disks induced by planet-disk interaction. Methods. The observability of debris disks with the Atacama Large Millimeter/submillimeter Array (ALMA) is studied on the basis of a simple analytical disk model. Furthermore, N-body simulations are used to model the spatial dust distribution in debris disks under the influence of planet-disk interaction. Images at optical scattered light to millimeter thermal re-emission are computed. Available information about the expected capabilities of ALMA and the James Webb Space Telescope (JWST) are used to investigate the observability of characteristic disk structures through spatially resolved imaging. Results. Planet-disk interaction can result in prominent structures. This provides the opportunity of detecting and characterizing extrasolar planets in a range of masses and radial distances from the star that is not accessible to other techniques. Facilities that will be available in the near future are shown to provide the capabilities to spatially resolve and characterize structures in debris disks. Limitations are revealed and suggestions for possible instrument setups and observing strategies are given. In particular, ALMA is limited by its sensitivity to surface brightness, which requires a trade-off between sensitivity and spatial resolution. Space-based midinfrared observations will be able to detect and spatially resolve regions in debris disks even at a distance of several tens of AU from the star, where the emission from debris disks in this wavelength range is expected to be low. [Abridged]

preprint2012arXiv

The Debris Disk Candidates: Eleven 24$μ$m excess stars in Spitzer SWIRE Fields

We present the optical to mid-infrared SEDs of 11 debris disk candidates from $Spitzer$ SWIRE fields. All these candidates are selected from SWIRE 24$μ$m sources matched with both the SDSS star catalog and the 2MASS point source catalog. They show an excess in the mid-infrared at 24$μ$m ($K_S$-[24]$_{Vega}$ $\ge$ 0.44) indicating the presence of a circumstellar dust disk. The observed optical spectra show that they are all late type main-sequence stars covering the spectral types of FGKM. Their fractional luminosities are well above 5$\times10^{-5}$, even up to the high fractional luminosity of 1$\times10^{-3}$. The high galactic latitudes of SWIRE fields indicate that most of these candidates could belong to the oldest stars in the thick disk. Our results indicate that the high fractional luminosity debris disks could exist in the old solar-like star systems, though they are now still quite rare. Their discoveries at high-galactic latitudes will also provide us an excellent opportunity to the further studies of properties and evolution of the debris disk in the ISM poor environments.

preprint2012arXiv

Tracing large-scale structures in circumstellar disks with ALMA

Planets are supposed to form in circumstellar disks. The gravitational potential of a planet perturbs the disk and leads to characteristic structures, i.e. spiral waves and gaps, in the disk's density profile. We perform a large-scale parameter study of the observability of these planet-induced structures in circumstellar disks with ALMA. On the basis of HD and MHD simulations, we calculated the disk temperature structure and (sub)mm images of these systems. These were used to derive simulated ALMA images. Because appropriate objects are frequent in Taurus, we focused on a distance of 140pc and a declination of 20°. The explored range of star-disk-planet configurations consists of 6 HD simulations (including magnetic fields and different planet masses), 9 disk sizes, 15 total disk masses, 6 different central stars, and two different grain size distributions. On almost all scales and in particular down to a scale of a few AU, ALMA is able to trace disk structures induced by planet-disk interaction or by the influence of magnetic fields on the wavelength range between 0.4 and 2.0mm. In most cases, the optimum angular resolution is limited by the sensitivity. However, within the range of typical masses of protoplanetary disks (0.1-0.001Msun) the disk mass has a minor impact on the observability. It is possible to resolve disks down to 2.67e-6Msun and trace gaps induced by a planet with M_p/M_s = 0.001 in disks with 2.67e-4Msun with a signal-to-noise ratio greater than three. The central star has a major impact on the observability of gaps, as well as the considered maximum grainsize of the dust in the disk. In general, it is more likely to trace planet-induced gaps in our MHD models, because gaps are wider in the presence of magnetic fields. We also find that zonal flows resulting from MRI create gap-like structures in the disk's re-emission radiation, which are observable with ALMA.

preprint2011arXiv

A Herschel Search For Cold Dust in Brown Dwarf Disks: First Results

We report initial results from a {\it Herschel} program to search for far-infrared emission from cold dust around a statistically significant sample of young brown dwarfs. The first three objects in our survey are all detected at 70\micron, and we report the first detection of a brown dwarf at 160\micron. The flux densities are consistent with the presence of substantial amounts of cold dust in the outer disks around these objects. We modeled the SED's with two different radiative transfer codes. We find that a broad range of model parameters provides a reasonable fit to the SED's, but that the addition of our 70\micron, and especially the 160\micron\ detection enables strong lower limits to be placed on the disk masses since most of the mass is in the outer disk. We find likely disk masses in the range of a few $\times 10^{-6}$ to $10^{-4}$ \msun. Our models provide a good fit to the SED's and do not require dust settling.

preprint2011arXiv

A High Mass Dusty Disk Candidate: The Case of IRAS 18151-1208

Many questions remain regarding the properties of disks around massive prototstars. Here we present the observations of a high mass protostellar object including an elongated dust continuum structure perpendicular to the outflow. Submillimeter Array 230 GHz line and continuum observations of the high mass protostellar object IRAS 18151-1208 along with single dish IRAM 30m observations afford us high spatial resolution (0.8") as well as recovery of the extended emission that gets filtered out by the interferometer. The observations of 12CO confirm the outflow direction to be in the southeast-northwest direction, and the 1.3 mm continuum exhibits an elongation in the direction perpendicular to the outflow. We model the physical parameters of the elongated structure by simultaneously fitting the observed spectral energy distribution (SED) and the brightness profile along the major axis using the 3D Radiative Transfer code MC3D. Assuming a density profile similar to that of a low mass disk, we can also reproduce the observations of this high mass protostellar object. This is achieved by using the same density distribution and flaring parameters as were used in the low mass case, and scaling up the size parameters that successfully modeled the circumstellar disk of several T Tauri stars. We also calculate that a region within the inner 30 AU of such a high mass disk is stable under the Toomre criterion. While we do not rule out other scenarios, we show here that the observations in the high mass regime are consistent with a scaled up version of a low mass disk. Implications on high mass star formation are discussed.

preprint2011arXiv

Mid-infrared observations of the transitional disks around DH Tau, DM Tau, and GM Aur

Aims: We present mid-infrared observations and photometry of the transitional disks around the young stellar objects DH Tau, DM Tau, and GM Aur, obtained with VISIR/VLT in N band. Our aim is to resolve the inner region and the large-scale structures of these transitional disks, carrying potential signatures of intermediate or later stages of disk evolution and ongoing planet formation. Methods: We use the simultaneously observed standard-stars as PSF reference to constrain the radial flux profiles of our target objects. Subtracting the obtained standard-star profile from the corresponding science object profile yields the flux residuals produced by the star-disk system. A detection threshold takes into account the background standard deviation and also the seeing variations during the observations to evaluate the significance of these flux residuals. On the basis of a simple model for the dust re-emission, we derive constraints on the inner radius of the dust disk. Results: We spatially resolve the transitional disk around GM Aur and determine an inner-disk hole radius of 20.5(+1.0,-0.5) AU. The circumstellar disks around DH Tau and DM Tau are not spatially resolved but we are able to constrain the inner-disk hole radius to <15.5(+9.0,-2.0) AU and <15.5(+0.5,-0.5) AU, respectively. The performed photometry yields fluxes of 178+-31 mJy for DH Tau, 56+-6 mJy for DM Tau, and 229+-14 mJy for GM Aur.

preprint2011arXiv

Multi-wavelength modeling of the spatially resolved debris disk of HD 107146

(abridged) We aim to constrain the location, composition, and dynamical state of planetesimal populations and dust around the young, sun-like (G2V) star HD 107146}. We consider coronagraphic observations obtained with the Advanced Camera for Surveys (HST/ACS) onboard the HST in broad V and broad I filters, a resolved 1.3mm map obtained with the Combined Array for Research in Millimeter-Wave Astronomy (CARMA), Spitzer/IRS low resolution spectra, and the spectral energy distribution (SED) of the object at wavelengths ranging from 3.5micron to 3.1mm. We complement these data with new coronagraphic high resolution observations of the debris disk using the Near Infrared Camera and Multi-Object Spectrometer (HST/NICMOS) aboard the HST in the F110W filter. The SED and images of the disk in scattered light as well as in thermal reemission are combined in our modeling using a parameterized model for the disk density distribution and optical properties of the dust. A detailed analytical model of the debris disk around HD 107146 is presented that allows us to reproduce the almost entire set of spatially resolved and unresolved multi-wavelength observations. Considering the variety of complementary observational data, we are able to break the degeneracies produced by modeling SED data alone. We find the disk to be an extended ring with a peak surface density at 131AU. Furthermore, we find evidence for an additional, inner disk probably composed of small grains released at the inner edge of the outer disk and moving inwards due to Poynting-Robertson drag. A birth ring scenario (i.e., a more or less broad ring of planetesimals creating the dust disk trough collisions) is found to be the most likely explanation of the ringlike shape of the disk.

preprint2011arXiv

Resolving the inner regions of the HD97048 circumstellar disk with VLT/NACO polarimetric differential imaging

Circumstellar disks are the cradles of planetary systems and their physical and chemical properties directly influence the planet formation process. As most planets supposedly form in the inner disk regions, i.e., within a few tens of AU, it is crucial to study circumstellar disk on these scales to constrain the conditions for planet formation. Our aims are to characterize the inner regions of the circumstellar disk around the young Herbig Ae/Be star HD97048 in polarized light. We use VLT/NACO to observe HD97048 in polarimetric differential imaging (PDI) mode in the H and Ks band. We spatially resolve the disk around HD97048 in polarized flux in both filters on scales between ~0.1"-1.0" corresponding to the inner ~16-160 AU. Fitting isophots to the flux calibrated H-band image between 13 - 14 mag/arcsec^2 and 14 - 15 mag/arcsec^2 we derive a apparent disk inclination angle of 34+-5 deg and 47+-2 deg, respectively. The disk position angle in both brightness regimes is almost identical and roughly 80 deg. Along the disk major axis the surface brightness of the polarized flux drops from ~11 mag/arcsec^2 at ~0.1" (~16 AU) to ~15.3 mag/arcsec^2 at ~1.0" (~160 AU). The brightness profiles along the major axis are fitted with power-laws falling off as ~r^(-1.78+-0.02) in H and ~r^(-2.34+-0.04) in Ks. As the surface brightness drops off more rapidly in Ks compared to H the disks becomes relatively bluer at larger separations possibly indicating changing dust grain properties as a function of radius. For the first time the inner ~0.1"-1.0" (~16-160 AU) of the surface layer of the HD97048 circumstellar disk have been imaged in scattered light demonstrating the power of ground-based imaging polarimetry. Our data fill an important gap in a large collection of existing data including resolved thermal dust and PAH emission images as well as resolved gas emission lines.

preprint2011arXiv

Transition disk chemistry and future prospects with ALMA

We explore the chemical structure of a disk that contains a large central gap of R ~ 45 AU, as is commonly seen in transitional disk systems. In our chemical model of a disk with a cleared inner void, the midplane becomes revealed to the central star so that it is directly irradiated. The midplane material at the truncation radius is permissive to reprocessed optical heating radiation, but opaque to the photo-dissociating ultraviolet, creating an environment abundant in gas-phase molecules. Thus the disk midplane, which would otherwise for a full disk be dominated by near complete heavy element freeze-out, should become observable in molecular emission. If this prediction is correct this has exciting prospects for observations with the Atacama Large Millimeter/Submillimeter Array (ALMA), as the inner transition region should thus be readily detected and resolved, especially using high-J rotational transitions excited in the high density midplane gas. Therefore such observations will potentially provide us with a direct probe of the physics and chemistry at this actively evolving interface.

preprint2011arXiv

VLT/NACO Polarimetric Differential Imaging of HD100546 - Disk Structure and Dust Grain Properties between 10-140 AU

We present polarimetric differential imaging (PDI) data of the circumstellar disk around the Herbig Ae/Be star HD100546 obtained with VLT/NACO. We resolve the disk in polarized light in the H and Ks filter between ~0.1-1.4" (i.e., ~10-140 AU). The innermost disk regions are directly imaged for the first time and the mean apparent disk inclination and position angle are derived. The surface brightness along the disk major axis drops off roughly with S(r) ~ r^-3 but has a maximum around 0.15" suggesting a marginal detection of the main disk inner rim at ~15 AU. We find a significant brightness asymmetry along the disk minor axis in both filters with the far side of the disk appearing brighter than the near side. This enhanced backward scattering and a low total polarization degree of the scattered disk flux of 14%(+19%/-8%) suggests that the dust grains on the disk surface are larger than typical ISM grains. Empirical scattering functions reveal the backward scattering peak at the largest scattering angles and a second maximum for the smallest scattering angles. This indicates a second dust grain population preferably forward scattering and smaller in size. It shows that, relatively, in the inner disk regions (40-50 AU) a higher fraction of larger grains is found compared to the outer disk regions (100-110 AU). Finally, our images reveal distinct substructures between 25-35 AU physical separation from the star and we discuss the possible origin for the two features in the context of ongoing planet formation.

preprint2010arXiv

Observing dust settling and coagulation in circumstellar discs: Selected constraints from high resolution imaging

Circumstellar discs are expected to be the nursery of planets. Grain growth within such discs is the first step in the planet formation process in the core-accretion gas-capture scenario. We aim at providing selected criteria on observational quantities derived from multi-wavelength imaging observations that allow to identify dust grain growth and settling. We define a wide-ranged parameter space of discs in various states of their evolution. Using a parametrised model set-up and radiative transfer techniques we compute multi-wavelength images of discs at different inclinations. Using millimetre and sub-millimetre images we are in the position to constrain the process of dust grain growth and sedimentation. However, the degeneracy between parameters prohibit the same achievement using near- to mid-infrared images. Using face-on observations in the N and Q Band, the sedimentation height can be constrained.

preprint2009arXiv

Dust amorphization in protoplanetary disks

High-energy irradiation of the circumstellar material might impact the structure and the composition of a protoplanetary disk and hence the process of planet formation. In this paper, we present a study on the possible influence of the stellar irradiation, indicated by X-ray emission, on the crystalline structure of the circumstellar dust. The dust crystallinity is measured for 42 class II T Tauri stars in the Taurus star-forming region using a decomposition fit of the 10 micron silicate feature, measured with the Spitzer IRS instrument. Since the sample includes objects with disks of various evolutionary stages, we further confine the target selection, using the age of the objects as a selection parameter. We correlate the X-ray luminosity and the X-ray hardness of the central object with the crystalline mass fraction of the circumstellar dust and find a significant anti-correlation for 20 objects within an age range of approx. 1 to 4.5 Myr. We postulate that X-rays represent the stellar activity and consequently the energetic ions of the stellar winds which interact with the circumstellar disk. We show that the fluxes around 1 AU and ion energies of the present solar wind are sufficient to amorphize the upper layer of dust grains very efficiently, leading to an observable reduction of the crystalline mass fraction of the circumstellar, sub-micron sized dust. This effect could also erase other relations between crystallinity and disk/star parameters such as age or spectral type.

preprint2006arXiv

A close look into the carbon disk at the core of the planetary nebula CPD-568032

We present high spatial resolution observations of the dusty core of the Planetary Nebula with Wolf-Rayet central star CPD-568032. These observations were taken with the mid-infrared interferometer VLTI/MIDI in imaging mode providing a typical 300 mas resolution and in interferometric mode using UT2-UT3 47m baseline providing a typical spatial resolution of 20 mas. The visible HST images exhibit a complex multilobal geometry dominated by faint lobes. The farthest structures are located at 7" from the star. The mid-IR environment of CPD-568032 is dominated by a compact source, barely resolved by a single UT telescope in a 8.7 micron filter. The infrared core is almost fully resolved with the three 40-45m projected baselines ranging from -5 to 51 degree but smooth oscillating fringes at low level have been detected in spectrally dispersed visibilities. This clear signal is interpreted in terms of a ring structure which would define the bright inner rim of the equatorial disk. Geometric models allowed us to derive the main geometrical parameters of the disk. For instance, a reasonably good fit is reached with an achromatic and elliptical truncated Gaussian with a radius of 97+/-11 AU, an inclination of 28+/-7 degree and a PA for the major axis at 345+/-7 degree. Furthermore, we performed some radiative transfer modeling aimed at further constraining the geometry and mass content of the disk, by taking into account the MIDI dispersed visibilities, spectra, and the large aperture SED of the source. These models show that the disk is mostly optically thin in the N band and highly flared.

preprint2006arXiv

Are Debris Disks and Massive Planets Correlated?

Using data from the Spitzer Space Telescope Legacy Science Program ``Formation and Evolution of Planetary Systems'' (FEPS), we have searched for debris disks around 9 FGK stars (2-10 Gyr), known from radial velocity (RV) studies to have one or more massive planets. Only one of the sources, HD 38529, has excess emission above the stellar photosphere; at 70 micron the signal-to-noise ratio in the excess is 4.7 while at wavelengths < 30 micron there is no evidence of excess. The remaining sources show no excesses at any Spitzer wavelengths. Applying survival tests to the FEPS sample and the results for the FGK survey published in Bryden et al. (2006), we do not find a significant correlation between the frequency and properties of debris disks and the presence of close-in planets. We discuss possible reasons for the lack of a correlation.

preprint2005arXiv

VITRUV - Science Cases

VITRUV is a second generation spectro-imager for the PRIMA enabled Very Large Telescope Interferometer. By combining simultaneously up to 8 telescopes VITRUV makes the VLTI up to 6 times more efficient. This operational gain allows two novel scientific methodologies: 1) massive surveys of sizes; 2) routine interferometric imaging. The science cases presented concentrate on the qualitatively new routine interferometric imaging methodology. The science cases are not exhaustive but complementary to the PRIMA reference mission. The focus is on: a) the close environment of young stars probing for the initial conditions of planet formation and disk evolution; b) the surfaces of stars tackling dynamos, activity, pulsation, mass-loss and evolution; c) revealing the origin of the extraordinary morphologies of Planetary Nebulae and related stars; d) studying the accretion-ejection structures of stellar black-holes (microquasars) in our galaxy; e) unveiling the different interacting components (torus, jets, BLRs) of Active Galactic Nuclei; and f) probing the environment of nearby supermassive black-holes and relativistic effects in the Galactic Center black-hole.