39 research outputs found
DISPATCH: A Numerical Simulation Framework for the Exa-scale Era. I. Fundamentals
We introduce a high-performance simulation framework that permits the
semi-independent, task-based solution of sets of partial differential
equations, typically manifesting as updates to a collection of `patches' in
space-time. A hybrid MPI/OpenMP execution model is adopted, where work tasks
are controlled by a rank-local `dispatcher' which selects, from a set of tasks
generally much larger than the number of physical cores (or hardware threads),
tasks that are ready for updating. The definition of a task can vary, for
example, with some solving the equations of ideal magnetohydrodynamics (MHD),
others non-ideal MHD, radiative transfer, or particle motion, and yet others
applying particle-in-cell (PIC) methods. Tasks do not have to be grid-based,
while tasks that are, may use either Cartesian or orthogonal curvilinear
meshes. Patches may be stationary or moving. Mesh refinement can be static or
dynamic. A feature of decisive importance for the overall performance of the
framework is that time steps are determined and applied locally; this allows
potentially large reductions in the total number of updates required in cases
when the signal speed varies greatly across the computational domain, and
therefore a corresponding reduction in computing time. Another feature is a
load balancing algorithm that operates `locally' and aims to simultaneously
minimise load and communication imbalance. The framework generally relies on
already existing solvers, whose performance is augmented when run under the
framework, due to more efficient cache usage, vectorisation, local
time-stepping, plus near-linear and, in principle, unlimited OpenMP and MPI
scaling.Comment: 17 pages, 8 figures. Accepted by MNRA
Ionization: a possible explanation for the difference of mean disk sizes in star-forming regions
Surveys of protoplanetary disks in star-forming regions of similar age
revealed significant variations in average disk mass between some regions. For
instance, disks in the Orion Nebular Cluster (ONC) and Corona Australis (CrA)
are on average smaller than disks observed in Lupus, Taurus, Chamaeleon I or
Ophiuchus. In contrast to previous models that study truncation of disks at a
late stage of their evolution, we investigate whether disks may already be born
with systematically smaller disk sizes in more massive star-forming regions as
a consequence of enhanced ionization rates. Assuming various cosmic-ray
ionization rates, we compute the resistivities for ambipolar diffusion and
Ohmic dissipation with a chemical network, and perform 2D non-ideal
magnetohydrodynamical protostellar collapse simulations. A higher ionization
rate leads to stronger magnetic braking, and hence to the formation of smaller
disks. Accounting for recent findings that protostars act as forges of cosmic
rays and considering only mild attenuation during the collapse phase, we show
that a high average cosmic-ray ionization rate in star-forming regions like the
ONC or CrA can explain the detection of smaller disks in these regions. Our
results show that on average a higher ionization rate leads to the formation of
smaller disks. Therefore, smaller disks in regions of similar age can be the
consequence of different levels of ionization, and may not exclusively be
caused by disk truncation via external photoevaporation. We strongly encourage
observations that allow measuring the cosmic-ray ionization degrees in
different star-forming regions to test our hypothesis.Comment: 9 pages, 5 figures, accepted for publication in A&
Late encounter-events as a source of disks and spiral structures -- Forming second generation disks
Both observations of arc-like structures and luminosity bursts of stars > 1
Myr in age indicate that at least some stars undergo late infall events. We
investigate scenarios of replenishing the mass reservoir around a star via
capturing and infalling events of cloudlets. We carry out altogether 24
three-dimensional hydrodynamical simulations of cloudlet encounters with a
Herbig star of mass 2.5 solar mass using the moving mesh code AREPO. To account
for the two possibilities of a star or a cloudlet traveling through the
interstellar medium (ISM), we put either the star or the cloudlet at rest with
respect to the background gas. For absent cooling in the adiabatic runs, almost
none of the cloudlet gas is captured due to high thermal pressure. However,
second-generation disks easily form when accounting for cooling of the gas. The
disk radii range from several 100 au to about 1000 au and associated arc-like
structures up to 10 000 au in length form around the star for runs with and
without stellar irradiation. Consistent with angular momentum conservation, the
arcs and disks are larger for larger impact parameters. Accounting for
turbulence in the cloudlet only mildly changes the model outcome. In the case
of the star being at rest with the background gas, the disk formation and mass
replenishment process is more pronounced and the associated arc-shaped
streamers are longer-lived. The results of our models confirm that late
encounter events lead to the formation of transitional disks associated with
arc-shaped structures such as observed for AB Aurigae or HD 100546. In
addition, we find that second-generation disks and their associated filamentary
arms are longer lived (>100 000 yrs) in infall events, when the star is at rest
with the background gas.Comment: 14 pages, 7 figures in main text and 2 figures in appendices A and B,
accepted for publication in A&
Linear dust polarization during the embedded phase of protostar formation
Measuring polarization from thermal dust emission can provide constraints on
the magnetic field structure around embedded protostars. However, interpreting
the observations is challenging without models that consistently account for
both the complexity of the protostellar birth environment and polarization
mechanisms. We aim to provide a better understanding with a focus on
bridge-like structures such as that observed towards the protostellar multiple
IRAS 16293--2422 by comparing synthetic polarization maps of thermal reemission
with observations. We analyze the magnetic field properties associated with the
formation of a protostellar multiple based on ideal MHD 3D zoom-in simulations
carried out with the RAMSES code. To compare with observations, we post-process
a snapshot of a bridge-like structure that is associated with a forming triple
star system with the radiative transfer code POLARIS and produce
multi-wavelength dust polarization maps. In the most prominent bridge of our
sample, the typical density is about 10^(-16) g cm^(-3), and the magnetic field
strength is about 1 to 2 mG. The magnetic field structure has an elongated
toroidal morphology and the dust polarization maps trace the complex
morphology. In contrast, the magnetic field strength associated with the
launching of asymmetric bipolar outflows is significantly more magnetized (~100
mG). At {\lambda}=1.3 mm, the orientation of grains in the bridge is similar
for the case accounting for radiative alignment torques (RATs) compared to
perfect alignment with magnetic field lines. However, the polarization fraction
in the bridge is three times smaller for the RAT scenario compared to assuming
perfect alignment. At shorter wavelengths ({\lambda} < 200 {\mu}m), dust
polarization does not trace the magnetic field because other effects such as
self-scattering and dichroic extinction dominate the orientation of the
polarization.Comment: 18 pages, 12 figures plus 3 figures in the appendix, accepted for
publication in A&
High-precision photometry by telescope defocussing - VIII.WASP-22, WASP-41,WASP-42 andWASP-55
We present 13 high-precision and four additional light curves of four bright southernhemisphere transiting planetary systems: WASP-22, WASP-41, WASP-42 and WASP-55. In the cases of WASP-42 and WASP-55, these are the first follow-up observations since their discovery papers. We present refined measurements of the physical properties and orbital ephemerides of all four systems. No indications of transit timing variations were seen. All four planets have radii inflated above those expected from theoretical models of gas-giant planets; WASP-55 b is the most discrepant with a mass of 0.63MJup and a radius of 1.34 RJup. WASP-41 shows brightness anomalies during transit due to the planet occulting spots on the stellar surface. Two anomalies observed 3.1 d apart are very likely due to the same spot. We measure its change in position and determine a rotation period for the host star of 18.6 ± 1.5 d, in good agreement with a published measurement from spot-induced brightness modulation, and a sky-projected orbital obliquity of λ = 6 ± 11°. We conclude with a compilation of obliquity measurements from spot-tracking analyses and a discussion of this technique in the study of the orbital configurations of hot Jupiters
An analysis of binary microlensing event OGLE-2015-BLG-0060
We present the analysis of stellar binary microlensing event OGLE-2015-BLG-0060 based on observations obtained from 13 different telescopes. Intensive coverage of the anomalous parts of the light curve was achieved by automated follow-up observations from the robotic telescopes of the Las Cumbres Observatory. We show that, for the first time, all main features of an anomalous microlensing event are well covered by follow-up data, allowing us to estimate the physical parameters of the lens. The strong detection of second-order effects in the event light curve necessitates the inclusion of longer-baseline survey data in order to constrain the parallax vector. We find that the event was most likely caused by a stellar binary-lens with masses M⋆1=0.87±0.12M⊙ and M⋆2=0.77±0.11M⊙. The distance to the lensing system is 6.41 ± 0.14 kpc and the projected separation between the two components is 13.85 ± 0.16 AU. Alternative interpretations are also considered
Spitzer Observations Of Ogle-2015-blg-1212 Reveal A New Path Toward Breaking Strong Microlens Degeneracies
Spitzer microlensing parallax observations of OGLE-2015-BLG-1212 decisively break a degeneracy between planetary and binary solutions that is somewhat ambiguous when only ground-based data are considered. Only eight viable models survive out of an initial set of 32 local minima in the parameter space. These models clearly indicate that the lens is a stellar binary system possibly located within the bulge of our Galaxy, ruling out the planetary alternative. We argue that several types of discrete degeneracies can be broken via such space-based parallax observations
The Spitzer Microlensing Program As A Probe For Globular Cluster Planets: Analysis Of Ogle-2015-BLG-0448
The microlensing event OGLE-2015-BLG-0448 was observed by Spitzer and lay within the tidal radius of the globular cluster NGC 6558. The event had moderate magnification and was intensively observed, hence it had the potential to probe the distribution of planets in globular clusters. We measure the proper motion of NGC 6558 () as well as the source and show that the lens is not a cluster member. Even though this particular event does not probe the distribution of planets in globular clusters, other potential cluster lens events can be verified using our methodology. Additionally, we find that microlens parallax measured using Optical Gravitational Lens Experiment (OGLE) photometry is consistent with the value found based on the light curve displacement between the Earth and Spitzer
High-resolution Imaging of Transiting Extrasolar Planetary systems (HITEP). II. Lucky Imaging results from 2015 and 2016
Context. The formation and dynamical history of hot Jupiters is currently debated, with wide stellar binaries having been suggested as a potential formation pathway. Additionally, contaminating light from both binary companions and unassociated stars can significantly bias the results of planet characterisation studies, but can be corrected for if the properties of the contaminating star are known. Aims. We search for binary companions to known transiting exoplanet host stars, in order to determine the multiplicity properties of hot Jupiter host stars. We also characterise unassociated stars along the line of sight, allowing photometric and spectroscopic observations of the planetary system to be corrected for contaminating light. Methods. We analyse lucky imaging observations of 97 Southern hemisphere exoplanet host stars, using the Two Colour Instrument on the Danish 1.54m telescope. For each detected companion star, we determine flux ratios relative to the planet host star in two passbands, and measure the relative position of the companion. The probability of each companion being physically associated was determined using our two-colour photometry. Results. A catalogue of close companion stars is presented, including flux ratios, position measurements, and estimated companion star temperature. For companions that are potential binary companions, we review archival and catalogue data for further evidence. For WASP-77AB and WASP-85AB, we combine our data with historical measurements to determine the binary orbits, showing them to be moderately eccentric and inclined to the line of sight and planetary orbital axis. Combining our survey with the similar Friends of Hot Jupiters survey, we conclude that known hot Jupiter host stars show a deficit of high mass stellar companions compared to the field star population; however, this may be a result of the biases in detection and target selection by ground-based surveys.PostprintPeer reviewe
High-precision photometry by telescope defocussing - VIII. WASP-22, WASP-41, WASP-42 and WASP-55
We present 13 high-precision and four additional light curves of four bright southern-hemisphere transiting planetary systems: WASP-22, WASP-41, WASP-42 and WASP-55. In the cases of WASP-42 and WASP-55, these are the first follow-up observations since their discovery papers. We present refined measurements of the physical properties and orbital ephemerides of all four systems. No indications of transit timing variations were seen. All four planets have radii inflated above those expected from theoretical models of gas-giant planets; WASP-55 b is the most discrepant with a mass of 0.63M(Jup) and a radius of 1.34 R-Jup. WASP-41 shows brightness anomalies during transit due to the planet occulting spots on the stellar surface. Two anomalies observed 3.1 d apart are very likely due to the same spot. We measure its change in position and determine a rotation period for the host star of 18.6 +/- 1.5 d, in good agreement with a published measurement from spot-induced brightness modulation, and a sky-projected orbital obliquity of lambda = 6 +/- 11 degrees. We conclude with a compilation of obliquity measurements from spot-tracking analyses and a discussion of this technique in the study of the orbital configurations of hot Jupiters