-
A Wigner molecule at extremely low densities: a numerically exact study
Authors:
Miguel Escobar Azor,
Léa Brooke,
Stefano Evangelisti,
Thierry Leininger,
Pierre-François Loos,
Nicolas Suaud,
J. A. Berger
Abstract:
In this work we investigate Wigner localization at very low densities by means of the exact diagonalization of the Hamiltonian. This yields numerically exact results. In particular, we study a quasi-one-dimensional system of two electrons that are confined to a ring by three-dimensional gaussians placed along the ring perimeter. To characterize the Wigner localization we study several appropriate…
▽ More
In this work we investigate Wigner localization at very low densities by means of the exact diagonalization of the Hamiltonian. This yields numerically exact results. In particular, we study a quasi-one-dimensional system of two electrons that are confined to a ring by three-dimensional gaussians placed along the ring perimeter. To characterize the Wigner localization we study several appropriate observables, namely the two-body reduced density matrix, the localization tensor and the particle-hole entropy. We show that the localization tensor is the most promising quantity to study Wigner localization since it accurately captures the transition from the delocalized to the localized state and it can be applied to systems of all sizes.
△ Less
Submitted 5 September, 2019; v1 submitted 4 July, 2019;
originally announced July 2019.
-
A phase of confined electroweak force in the early Universe
Authors:
Joshua Berger,
Andrew J. Long,
Jessica Turner
Abstract:
We consider a modified cosmological history in which the presence of beyond-the-Standard-Model physics causes the weak gauge sector, $\mathrm{SU}(2)_L$, to confine before it is Higgsed. Under the assumption of chiral symmetry breaking, quark and lepton weak doublets form condensates that break the global symmetries of the Standard Model, including baryon and lepton number, down to a…
▽ More
We consider a modified cosmological history in which the presence of beyond-the-Standard-Model physics causes the weak gauge sector, $\mathrm{SU}(2)_L$, to confine before it is Higgsed. Under the assumption of chiral symmetry breaking, quark and lepton weak doublets form condensates that break the global symmetries of the Standard Model, including baryon and lepton number, down to a $\mathrm{U}(1)$ subgroup under which only the weak singlet fermions and Higgs transform. The weakly-coupled gauge group $\mathrm{SU}(3)_c \times \mathrm{U}(1)_Y$ is also broken to an $\mathrm{SU}(2)_c \times \mathrm{U}(1)_Q$ gauge group. The light states include (pseudo-)Goldstone bosons of the global symmetry breaking, mostly-elementary fermions primarily composed of the weak singlet quarks and leptons, and the gauge bosons of the weakly-coupled gauge group. We discuss possible signatures from early Universe cosmology including gravitational wave radiation, topological defects, and baryogenesis.
△ Less
Submitted 5 September, 2019; v1 submitted 12 June, 2019;
originally announced June 2019.
-
On the comparison of three numerical methods applied to building simulation: finite-differences, RC circuit approximation and a spectral method
Authors:
Julien Berger,
Suelen Gasparin,
Denys Dutykh,
Nathan Mendes
Abstract:
Predictions of physical phenomena in buildings are carried out by using physical models formulated as a mathematical problem and solved by means of numerical methods, aiming at evaluating, for instance, the building thermal or hygrothermal performance by calculating distributions and fluxes of heat and moisture transfer. Therefore, the choice of the numerical method is crucial since it is a compro…
▽ More
Predictions of physical phenomena in buildings are carried out by using physical models formulated as a mathematical problem and solved by means of numerical methods, aiming at evaluating, for instance, the building thermal or hygrothermal performance by calculating distributions and fluxes of heat and moisture transfer. Therefore, the choice of the numerical method is crucial since it is a compromise among (i) the solution accuracy, (ii) the computational cost to obtain the solution and (iii) the complexity of the method implementation. An efficient numerical method enables to compute an accurate solution with a minimum computational run time (CPU). On that account, this article brings an investigation on the performance of three numerical methods. The first one is the standard and widely used finite-difference approach, while the second one is the so-called RC approach, which is a particular method brought to the building physics area by means of an analogy of electric circuits. The third numerical method is the spectral one, which has been recently proposed to solve nonlinear diffusive problems in building physics. The three methods are evaluated in terms of accuracy on the assessment of the dependent variable (temperature or vapor pressure) or of density of fluxes for three different cases: i) heat diffusion through a concrete slab, ii) moisture diffusion through an aerated concrete slab and iii) heat diffusion using measured temperatures as boundary conditions. Results highlight the spectral approach as the most accurate method. The RC based model with a few number of resistances does not provide accurate results for temperature and vapor pressure distributions neither to flux densities nor conduction loads.
△ Less
Submitted 29 May, 2019;
originally announced May 2019.
-
Evaluation of the reliability of building energy performance models for parameter estimation
Authors:
Julien Berger,
Denys Dutykh
Abstract:
The fidelity of a model relies both on its accuracy to predict the physical phenomena and its capability to estimate unknown parameters using observations. This article focuses on this second aspect by analyzing the reliability of two mathematical models proposed in the literature for the simulation of heat losses through building walls. The first one, named DuFort-Frankel (DF), is the classical h…
▽ More
The fidelity of a model relies both on its accuracy to predict the physical phenomena and its capability to estimate unknown parameters using observations. This article focuses on this second aspect by analyzing the reliability of two mathematical models proposed in the literature for the simulation of heat losses through building walls. The first one, named DuFort-Frankel (DF), is the classical heat diffusion equation combined with the DuFort-Frankel numerical scheme. The second is the so-called RC lumped approach, based on a simple ordinary differential equation to compute the temperature within the wall. The reliability is evaluated following a two stages method. First, samples of observations are generated using a pseudo-spectral numerical model for the heat diffusion equation with known input parameters. The results are then modified by adding a noise to simulate experimental measurements. Then, for each sample of observation, the parameter estimation problem is solved using one of the two mathematical models. The reliability is assessed based on the accuracy of the approach to recover the unknown parameter. Three case studies are considered for the estimation of (i) the heat capacity, (ii) the thermal conductivity or (iii) the heat transfer coefficient at the interface between the wall and the ambient air. For all cases, the DF mathematical model has a very satisfactory reliability to estimate the unknown parameters without any bias. However, the RC model lacks of fidelity and reliability. The error on the estimated parameter can reach 40% for the heat capacity, 80% for the thermal conductivity and 450% for the heat transfer coefficient.
△ Less
Submitted 29 May, 2019;
originally announced May 2019.
-
A geometric distance measurement to the Galactic Center black hole with 0.3% uncertainty
Authors:
R. Abuter,
A. Amorim,
M. Bauboeck,
J. P. Berger,
H. Bonnet,
W. Brandner,
Y. Clenet,
V. Coude du Foresto,
P. T. de Zeeuw,
J. Dexter,
G. Duvert,
A. Eckart,
F. Eisenhauer,
N. M. Foerster Schreiber,
P. Garcia,
F. Gao,
E. Gendron,
R. Genzel,
O. Gerhard,
S. Gillessen,
M. Habibi,
X. Haubois,
T. Henning,
S. Hippler,
M. Horrobin
, et al. (29 additional authors not shown)
Abstract:
We present a 0.16% precise and 0.27% accurate determination of R0, the distance to the Galactic Center. Our measurement uses the star S2 on its 16-year orbit around the massive black hole Sgr A* that we followed astrometrically and spectroscopically for 27 years. Since 2017, we added near-infrared interferometry with the VLTI beam combiner GRAVITY, yielding a direct measurement of the separation v…
▽ More
We present a 0.16% precise and 0.27% accurate determination of R0, the distance to the Galactic Center. Our measurement uses the star S2 on its 16-year orbit around the massive black hole Sgr A* that we followed astrometrically and spectroscopically for 27 years. Since 2017, we added near-infrared interferometry with the VLTI beam combiner GRAVITY, yielding a direct measurement of the separation vector between S2 and Sgr A* with an accuracy as good as 20 micro-arcsec in the best cases. S2 passed the pericenter of its highly eccentric orbit in May 2018, and we followed the passage with dense sampling throughout the year. Together with our spectroscopy, in the best cases with an error of 7 km/s, this yields a geometric distance estimate: R0 = 8178 +- 13(stat.) +- 22(sys.) pc. This work updates our previous publication in which we reported the first detection of the gravitational redshift in the S2 data. The redshift term is now detected with a significance level of 20 sigma with f_redshift = 1.04 +- 0.05.
△ Less
Submitted 11 April, 2019;
originally announced April 2019.
-
Neuralogram: A Deep Neural Network Based Representation for Audio Signals
Authors:
Prateek Verma,
Chris Chafe,
Jonathan Berger
Abstract:
We propose the Neuralogram -- a deep neural network based representation for understanding audio signals which, as the name suggests, transforms an audio signal to a dense, compact representation based upon embeddings learned via a neural architecture. Through a series of probing signals, we show how our representation can encapsulate pitch, timbre and rhythm-based information, and other attribute…
▽ More
We propose the Neuralogram -- a deep neural network based representation for understanding audio signals which, as the name suggests, transforms an audio signal to a dense, compact representation based upon embeddings learned via a neural architecture. Through a series of probing signals, we show how our representation can encapsulate pitch, timbre and rhythm-based information, and other attributes. This representation suggests a method for revealing meaningful relationships in arbitrarily long audio signals that are not readily represented by existing algorithms. This has the potential for numerous applications in audio understanding, music recommendation, meta-data extraction to name a few.
△ Less
Submitted 10 April, 2019;
originally announced April 2019.
-
First direct detection of an exoplanet by optical interferometry; Astrometry and K-band spectroscopy of HR8799 e
Authors:
S. Lacour,
M. Nowak,
J. Wang,
O. Pfuhl,
F. Eisenhauer,
R. Abuter,
A. Amorim,
N. Anugu,
M. Benisty,
J. P. Berger,
H. Beust,
N. Blind,
M. Bonnefoy,
H. Bonnet,
P. Bourget,
W. Brandner,
A. Buron,
C. Collin,
B. Charnay,
F. Chapron,
Y. Clénet,
V. Coudé du Foresto,
P. T. de Zeeuw,
C. Deen,
R. Dembet
, et al. (63 additional authors not shown)
Abstract:
To date, infrared interferometry at best achieved contrast ratios of a few times $10^{-4}$ on bright targets. GRAVITY, with its dual-field mode, is now capable of high contrast observations, enabling the direct observation of exoplanets. We demonstrate the technique on HR8799, a young planetary system composed of four known giant exoplanets. We used the GRAVITY fringe tracker to lock the fringes o…
▽ More
To date, infrared interferometry at best achieved contrast ratios of a few times $10^{-4}$ on bright targets. GRAVITY, with its dual-field mode, is now capable of high contrast observations, enabling the direct observation of exoplanets. We demonstrate the technique on HR8799, a young planetary system composed of four known giant exoplanets. We used the GRAVITY fringe tracker to lock the fringes on the central star, and integrated off-axis on the HR8799e planet situated at 390 mas from the star. Data reduction included post-processing to remove the flux leaking from the central star and to extract the coherent flux of the planet. The inferred K band spectrum of the planet has a spectral resolution of 500. We also derive the astrometric position of the planet relative to the star with a precision on the order of 100$\,μ$as. The GRAVITY astrometric measurement disfavors perfectly coplanar stable orbital solutions. A small adjustment of a few degrees to the orbital inclination of HR 8799 e can resolve the tension, implying that the orbits are close to, but not strictly coplanar. The spectrum, with a signal-to-noise ratio of $\approx 5$ per spectral channel, is compatible with a late-type L brown dwarf. Using Exo-REM synthetic spectra, we derive a temperature of $1150\pm50$\,K and a surface gravity of $10^{4.3\pm0.3}\,$cm/s$^{2}$. This corresponds to a radius of $1.17^{+0.13}_{-0.11}\,R_{\rm Jup}$ and a mass of $10^{+7}_{-4}\,M_{\rm Jup}$, which is an independent confirmation of mass estimates from evolutionary models. Our results demonstrate the power of interferometry for the direct detection and spectroscopic study of exoplanets at close angular separations from their stars.
△ Less
Submitted 28 March, 2019;
originally announced March 2019.
-
Many-Body Effective Energy Theory: photoemission at strong correlation
Authors:
S. Di Sabatino,
J. A. Berger,
P. Romaniello
Abstract:
In this work we explore the performance of a recently derived many-body effective energy theory for the calculation of photoemission spectra in the regime of strong electron correlation. We apply the theory to paramagnetic MnO, FeO, CoO, and NiO, which are typical examples of strongly correlated materials and, therefore, a challenge for standard theories. We show that our methods open a correlatio…
▽ More
In this work we explore the performance of a recently derived many-body effective energy theory for the calculation of photoemission spectra in the regime of strong electron correlation. We apply the theory to paramagnetic MnO, FeO, CoO, and NiO, which are typical examples of strongly correlated materials and, therefore, a challenge for standard theories. We show that our methods open a correlation gap in all the oxides studied without breaking the symmetry. Although the materials seem similar, we show that an analysis of the occupation numbers reveals that the nature of the gap is not the same for these materials. Overall the results are very promising, although improvements are clearly required, since the band gap is overestimated for all the systems studied. We indicate some possible strategies to further develop the theory.
△ Less
Submitted 22 March, 2019;
originally announced March 2019.
-
A new model for simulating heat, air and moisture transport in porous building materials
Authors:
Julien Berger,
Denys Dutykh,
Nathan Mendes,
Bolatbek Rysbaiuly
Abstract:
This work presents a detailed mathematical model combined with an innovative efficient numerical model to predict heat, air and moisture transfer through porous building materials. The model considers the transient effects of air transport and its impact on the heat and moisture transfer. The achievement of the mathematical model is detailed in the continuity of Luikov's work. A system composed of…
▽ More
This work presents a detailed mathematical model combined with an innovative efficient numerical model to predict heat, air and moisture transfer through porous building materials. The model considers the transient effects of air transport and its impact on the heat and moisture transfer. The achievement of the mathematical model is detailed in the continuity of Luikov's work. A system composed of two advection-diffusion differential equations plus one exclusively diffusion equation is derived. The main issue to take into account the transient air transfer arises in the very small characteristic time of the transfer, implying very fine discretisation. To circumvent these difficulties, the numerical model is based on the Du Fort-Frankel explicit and unconditionally stable scheme for the exclusively diffusion equation. It is combined with a two-step Runge-Kutta scheme in time with the Scharfetter-Gummel numerical scheme in space for the coupled advection-diffusion equations. At the end, the numerical model enables to relax the stability condition, and, therefore, to save important computational efforts. A validation case is considered to evaluate the efficiency of the model for a nonlinear problem. Results highlight a very accurate solution computed about 16 times faster than standard approaches. After this numerical validation, the reliability of the mathematical model is evaluated by comparing the numerical predictions to experimental observations. The latter is measured within a multi-layered wall submitted to a sudden increase of vapor pressure on the inner side and driven climate boundary conditions on the outer side. A very satisfactory agreement is noted between the numerical predictions and experimental observations indicating an overall good reliability of the proposed model.
△ Less
Submitted 7 March, 2019;
originally announced March 2019.
-
An innovative method to determine optimum insulation thickness based on non-uniform adaptive moving grid
Authors:
Suelen Gasparin,
Julien Berger,
Denys Dutykh,
Nathan Mendes
Abstract:
It is well known that thermal insulation is a leading strategy for reducing energy consumption associated to heating or cooling processes in buildings. Nevertheless, building insulation can generate high expenditures so that the selection of an optimum insulation thickness requires a detailed energy simulation as well as an economic analysis. In this way, the present study proposes an innovative n…
▽ More
It is well known that thermal insulation is a leading strategy for reducing energy consumption associated to heating or cooling processes in buildings. Nevertheless, building insulation can generate high expenditures so that the selection of an optimum insulation thickness requires a detailed energy simulation as well as an economic analysis. In this way, the present study proposes an innovative non-uniform adaptive method to determine the optimal insulation thickness of external walls. First, the method is compared with a reference solution to properly understand the features of the method, which can provide high accuracy with less spatial nodes. Then, the adaptive method is used to simulate the transient heat conduction through the building envelope of buildings located in Brazil, where there is a large potential of energy reduction. Simulations have been efficiently carried out for different wall and roof configurations, showing that the innovative method efficiently provides a gain of 25% on the computer run time.
△ Less
Submitted 26 February, 2019;
originally announced March 2019.
-
An efficient numerical model for liquid water uptake in porous material and its parameter estimation
Authors:
Ainagul Jumabekova,
Julien Berger,
Denys Dutykh,
Hervé Le Meur,
Aurélie Foucquier,
Mickaël Pailha,
Christophe Ménézo
Abstract:
The goal of this study is to propose an efficient numerical model for the predictions of capillary adsorption phenomena in a porous material. The Scharfetter-Gummel numerical scheme is proposed to solve an advection-diffusion equation with gravity flux. Its advantages such as accuracy, relaxed stability condition, and reduced computational cost are discussed along with the study of linear and nonl…
▽ More
The goal of this study is to propose an efficient numerical model for the predictions of capillary adsorption phenomena in a porous material. The Scharfetter-Gummel numerical scheme is proposed to solve an advection-diffusion equation with gravity flux. Its advantages such as accuracy, relaxed stability condition, and reduced computational cost are discussed along with the study of linear and nonlinear cases. The reliability of the numerical model is evaluated by comparing the numerical predictions with experimental observations of liquid uptake in bricks. A parameter estimation problem is solved to adjust the uncertain coefficients of moisture diffusivity and hydraulic conductivity.
△ Less
Submitted 27 February, 2019;
originally announced March 2019.
-
An adaptive simulation of nonlinear heat and moisture transfer as a boundary value problem
Authors:
Suelen Gasparin,
Julien Berger,
Denys Dutykh,
Nathan Mendes
Abstract:
This work presents an alternative view on the numerical simulation of diffusion processes applied to the heat and moisture transfer through porous building materials. Traditionally, by using the finite-difference approach, the discretization follows the Method Of Lines (MOL), when the problem is first discretized in space to obtain a large system of coupled Ordinary Differential Equations (ODEs).…
▽ More
This work presents an alternative view on the numerical simulation of diffusion processes applied to the heat and moisture transfer through porous building materials. Traditionally, by using the finite-difference approach, the discretization follows the Method Of Lines (MOL), when the problem is first discretized in space to obtain a large system of coupled Ordinary Differential Equations (ODEs). Thus, this paper proposes to change this viewpoint. First, we discretize in time to obtain a small system of coupled ODEs, which means instead of having a Cauchy (Initial Value) Problem (IVP), we have a Boundary Value Problem (BVP). Fortunately, BVPs can be solved efficiently today using adaptive collocation methods of high order. To demonstrate the benefits of this new approach, three case studies are presented, in which one of them is compared with experimental data. The first one considers nonlinear heat and moisture transfer through one material layer while the second one considers two material layers. Results show how the nonlinearities and the interface between materials are easily treated, by reasonably using a fourth-order adaptive method. Finally, the last case study compares numerical results with experimental measurements, showing a good agreement.
△ Less
Submitted 18 February, 2019;
originally announced February 2019.
-
Test of Einstein equivalence principle near the Galactic center supermassive black hole
Authors:
A. Amorim,
M. Bauböck,
J. P. Berger,
W. Brandner,
Y. Clénet,
V. Coudé du Foresto,
P. T. de Zeeuw,
J. Dexter,
G. Duvert,
M. Ebert,
A. Eckart,
F. Eisenhauer,
N. M. Förster Schreiber,
P. Garcia,
F. Gao,
E. Gendron,
R. Genzel,
S. Gillessen,
M. Habibi,
X. Haubois,
Th. Henning,
S. Hippler,
M. Horrobin,
Z. Hubert,
A. Jiménez Rosales
, et al. (27 additional authors not shown)
Abstract:
During its orbit around the four million solar mass black hole Sagittarius A* the star S2 experiences significant changes in gravitational potential. We use this change of potential to test one part of the Einstein equivalence principle: the local position invariance (LPI). We study the dependency of different atomic transitions on the gravitational potential to give an upper limit on violations o…
▽ More
During its orbit around the four million solar mass black hole Sagittarius A* the star S2 experiences significant changes in gravitational potential. We use this change of potential to test one part of the Einstein equivalence principle: the local position invariance (LPI). We study the dependency of different atomic transitions on the gravitational potential to give an upper limit on violations of the LPI. This is done by separately measuring the redshift from hydrogen and helium absorption lines in the stellar spectrum during its closest approach to the black hole. For this measurement we use radial velocity data from 2015 to 2018 and combine it with the gravitational potential at the position of S2, which is calculated from the precisely known orbit of S2 around the black hole. This results in a limit on a violation of the LPI of $|β_{He}-β_{H}| = (2.4 \pm 5.1) \cdot 10^{-2}$. The variation in potential that we probe with this measurement is six magnitudes larger than possible for measurements on Earth, and a factor ten larger than in experiments using white dwarfs. We are therefore testing the LPI in a regime where it has not been tested before.
△ Less
Submitted 25 March, 2019; v1 submitted 11 February, 2019;
originally announced February 2019.
-
A discontinuity in the $T_{\rm eff}$-radius relation of M-dwarfs
Authors:
Markus Rabus,
Régis Lachaume,
Andrés Jordán,
Rafael Brahm,
Tabetha Boyajian,
Kaspar von Braun,
Néstor Espinoza,
Jean-Philippe Berger,
Jean-Baptiste Le Bouquin,
Olivier Absil
Abstract:
We report on 13 new high-precision measurements of stellar diameters for low-mass dwarfs obtained by means of near-infrared long-baseline interferometry with PIONIER at the Very Large Telescope Interferometer. Together with accurate parallaxes from Gaia DR2, these measurements provide precise estimates for their linear radii, effective temperatures, masses, and luminosities. This allows us to refi…
▽ More
We report on 13 new high-precision measurements of stellar diameters for low-mass dwarfs obtained by means of near-infrared long-baseline interferometry with PIONIER at the Very Large Telescope Interferometer. Together with accurate parallaxes from Gaia DR2, these measurements provide precise estimates for their linear radii, effective temperatures, masses, and luminosities. This allows us to refine the effective temperature scale, in particular towards the coolest M-dwarfs. We measure for late-type stars with enhanced metallicity slightly inflated radii, whereas for stars with decreased metallicity we measure smaller radii. We further show that Gaia DR2 effective temperatures for M-dwarfs are underestimated by $\sim$ 8.2 % and give an empirical $M_{G}$-$T_{\rm eff}$ relation which is better suited for M-dwarfs with $T_{\rm eff}$ between 2600 and 4000 K. Most importantly, we are able to observationally identify a discontinuity in the $T_{\rm eff}$-radius plane, which is likely due to the transition from partially convective M-dwarfs to the fully convective regime. We found this transition to happen between 3200 K and 3340 K, or equivalently for stars with masses $\approx 0.23 M_{\odot}$. We find that in this transition region the stellar radii are in the range from 0.18 to 0.42$R_{\odot}$ for similar stellar effective temperatures.
△ Less
Submitted 23 January, 2019;
originally announced January 2019.
-
Towards reliable uncertainties in IR interferometry: The bootstrap for correlated statistical & systematic errors
Authors:
Régis Lachaume,
Markus Rabus,
Andrés Jordán,
Rafael Brahm,
Tabetha Boyajian,
Kaspar von Braun,
Jean-Philippe Berger
Abstract:
We propose a method to overcome the usual limitation of current data processing techniques in optical and infrared long-baseline interferometry: most reduction pipelines assume uncorrelated statistical errors and ignore systematics. We use the bootstrap method to sample the multivariate probability density function of the interferometric observables. It allows us to determine the correlations betw…
▽ More
We propose a method to overcome the usual limitation of current data processing techniques in optical and infrared long-baseline interferometry: most reduction pipelines assume uncorrelated statistical errors and ignore systematics. We use the bootstrap method to sample the multivariate probability density function of the interferometric observables. It allows us to determine the correlations between statistical error terms and their deviation from a Gaussian distribution. In addition, we introduce systematics as an additional, highly correlated error term whose magnitude is chosen to fit the data dispersion.
We have applied the method to obtain accurate measurements of stellar diameters for under-resolved stars, i.e. smaller than the angular resolution of the interferometer. We show that taking correlations and systematics has a significant impact on both the diameter estimate and its uncertainty. The robustness of our diameter determination comes at a price: we obtain 4 times larger uncertainties, of a few percent for most stars in our sample.
△ Less
Submitted 9 January, 2019;
originally announced January 2019.
-
Modified traces for quasi-Hopf algebras
Authors:
Johannes Berger,
Azat M. Gainutdinov,
Ingo Runkel
Abstract:
Let H be a finite-dimensional unimodular pivotal quasi-Hopf algebra over a field k, and let H-mod be the pivotal tensor category of finite-dimensional H-modules. We give a bijection between left (resp. right) modified traces on the tensor ideal H-pmod of projective modules and left (resp. right) cointegrals for H. The non-zero left/right modified traces are non-degenerate, and we show that non-deg…
▽ More
Let H be a finite-dimensional unimodular pivotal quasi-Hopf algebra over a field k, and let H-mod be the pivotal tensor category of finite-dimensional H-modules. We give a bijection between left (resp. right) modified traces on the tensor ideal H-pmod of projective modules and left (resp. right) cointegrals for H. The non-zero left/right modified traces are non-degenerate, and we show that non-degenerate left/right modified traces can only exist for unimodular H. This generalises results of Beliakova, Blanchet, and Gainutdinov from Hopf algebras to quasi-Hopf algebras. As an example we compute cointegrals and modified traces for the family of symplectic fermion quasi-Hopf algebras.
△ Less
Submitted 11 December, 2019; v1 submitted 26 December, 2018;
originally announced December 2018.
-
Dark matter amnesia in out-of-equilibrium scenarios
Authors:
Joshua Berger,
Djuna Croon,
Sonia El Hedri,
Karsten Jedamzik,
Ashley Perko,
Devin G. E. Walker
Abstract:
Models in which the dark matter is produced at extremely low rates from the annihilation of Standard Model particles in the early Universe allow us to explain the current dark matter relic density while easily evading the traditional experimental constraints. In scenarios where the dark matter interacts with the Standard Model via a new physics mediator, the early Universe dynamics of the dark sec…
▽ More
Models in which the dark matter is produced at extremely low rates from the annihilation of Standard Model particles in the early Universe allow us to explain the current dark matter relic density while easily evading the traditional experimental constraints. In scenarios where the dark matter interacts with the Standard Model via a new physics mediator, the early Universe dynamics of the dark sector can be particularly complex, as the dark matter and the mediator could be in thermal and chemical equilibrium with each other. This equilibration takes place via number-changing processes such as double Compton scattering and bremsstrahlung, whose amplitudes are cumbersome to calculate. In this paper, we show that in large regions of the parameter space, these equilibration mechanisms do not significantly affect the final dark matter relic density. In particular, for a model with a light dark photon mediator, the relic density can be reasonably estimated by considering that the dark matter is solely produced through the annihilation of Standard Model particles. This result considerably simplifies the treatment of a large class of dark matter theories, facilitating in particular the superimposition of the relic density constraints on the current and future experimental bounds.
△ Less
Submitted 20 December, 2018;
originally announced December 2018.
-
A Module For Boosted Dark Matter Event Generation in GENIE
Authors:
Joshua Berger
Abstract:
Models that produce a flux of semi-relativistic or relativistic boosted dark matter at large neutrino detectors are well-motivated extensions beyond the minimal weakly interacting massive particle (WIMP) paradigm. Current and upcoming liquid argon time projection chamber (LArTPC) based detectors will have improved sensitivity to such models, but also require improved theoretical modeling to better…
▽ More
Models that produce a flux of semi-relativistic or relativistic boosted dark matter at large neutrino detectors are well-motivated extensions beyond the minimal weakly interacting massive particle (WIMP) paradigm. Current and upcoming liquid argon time projection chamber (LArTPC) based detectors will have improved sensitivity to such models, but also require improved theoretical modeling to better understand their signals and optimize their analyses. I present the first full Monte Carlo tool for boosted dark matter interacting with nuclei in the energy regime accessible to LArTPC detectors, including the Deep Underground Neutrino Experiment (DUNE). The code uses the nuclear and strong physics modeling of the GENIE neutrino Monte Carlo event generator with particle physics modeling for dark matter. The code will be available in GENIE v3. In addition, I present a code for generating a GENIE-compatible flux of boosted dark matter coming from the Sun that is released independently.
△ Less
Submitted 13 December, 2018;
originally announced December 2018.
-
Identification of two-dimensional $FeO_2$ termination of hematite $α-Fe_2O_3(0001)$ surface
Authors:
Jesús Redondo,
Petr Lazar,
Pavel Procházka,
Stanislav Průša,
Jan Lachnitt,
Benjamín Mallada,
Aleš Cahlík,
Jan Berger,
Břetislav Šmíd,
Pavel Jelínek,
Jan Čechal,
Martin Švec
Abstract:
Iron oxides are among the most abundant compounds on Earth and have consequently been studied and used extensively in industrial processes. Despite these efforts, concrete understanding of some of their surface phase structures has remained elusive, in particular the oxidized $α-Fe_2O_3(0001)$ hematite surface. We detail an optimized recipe to produce this phase over the entire hematite surface an…
▽ More
Iron oxides are among the most abundant compounds on Earth and have consequently been studied and used extensively in industrial processes. Despite these efforts, concrete understanding of some of their surface phase structures has remained elusive, in particular the oxidized $α-Fe_2O_3(0001)$ hematite surface. We detail an optimized recipe to produce this phase over the entire hematite surface and study the geometrical parameters and composition of its complex structure by means of atomically resolved microscopy, electron diffraction and surface-sensitive spectroscopies. We conclude that the oxidized $α-Fe_2O_3(0001)$ surface is terminated by a two-dimensional iron oxide with structure, lattice parameters, and orientation different from the bulk substrate. Using total-energy density functional theory for simulation of a large-scale atomic model, we identify the structure of the surface layer as antiferromagnetic, conductive $1T-FeO_2$ attached on half-metal terminated bulk. The model succeeds in reproducing the characteristic modulations observed in the atomically resolved images and electron diffraction patterns.
△ Less
Submitted 25 November, 2018;
originally announced November 2018.
-
The principle underlying antiaromaticity
Authors:
Raphael J. F. Berger,
Alexandre Viel
Abstract:
Aromaticity is one of the most widely used chemical concepts. Current definitions are purely phenomenological and relate symmetry, reactive stability and the occurence of molecular diamagnetic response currents. The antithetical concept of antiaromaticity provides a connection between the contrary properties of structural instability or distortion out of higher symmetry, a small HOMO-LUMO gap, and…
▽ More
Aromaticity is one of the most widely used chemical concepts. Current definitions are purely phenomenological and relate symmetry, reactive stability and the occurence of molecular diamagnetic response currents. The antithetical concept of antiaromaticity provides a connection between the contrary properties of structural instability or distortion out of higher symmetry, a small HOMO-LUMO gap, and paramagnetic response currents. We reveal the principle that is underlying antiaromaticity in showing an intimate and strict symmetry induced relation between these properties. This principle can be proven and is formulated like: First order (and related) Jahn-Teller distorted molecules out of non-isometric point groups are prone to paramagnetic current susceptibility parallel to the main axis of symmetry. We show by the exemplary cases of cyclobutadiene, cylcooctatetraene, pentalene and manganese trifluoride how this principle works and discuss this new perspective on antiaromaticity.
△ Less
Submitted 22 February, 2019; v1 submitted 21 November, 2018;
originally announced November 2018.
-
A simple position operator for periodic systems
Authors:
Emilia Aragao Valença,
Diego Moreno,
Stefano Battaglia,
Gian Luigi Bendazzoli,
Stefano Evangelisti,
Thierry Leininger,
Nicolas Suaud,
J. A. Berger
Abstract:
We present a position operator that is compatible with periodic boundary conditions (PBC). It is a one-body operator that can be applied in calculations of correlated materials by simply replacing the traditional position vector by the new definition. We show that it satisfies important fundamental as well as practical constraints. To illustrate the usefulness of the PBC position operator we apply…
▽ More
We present a position operator that is compatible with periodic boundary conditions (PBC). It is a one-body operator that can be applied in calculations of correlated materials by simply replacing the traditional position vector by the new definition. We show that it satisfies important fundamental as well as practical constraints. To illustrate the usefulness of the PBC position operator we apply it to the localization tensor, a key quantity that is able to differentiate metallic from insulating states. In particular, we show that the localization tensor given in terms of the PBC position operator yields the correct expression in the thermodynamic limit. Moreover, we show that it correctly distinguishes between finite precursors of metals and insulators.
△ Less
Submitted 13 November, 2018;
originally announced November 2018.
-
Probing the Inner Disk Emission of the Herbig Ae Stars HD 163296 and HD 190073
Authors:
Benjamin R. Setterholm,
John D. Monnier,
Claire L. Davies,
Alexander Kreplin,
Stefan Kraus,
Fabien Baron,
Alicia Aarnio,
Jean-Philippe Berger,
Nuria Calvet,
Michel Curé,
Samer Kanaan,
Brian Kloppenborg,
Jean-Baptiste Le Bouquin,
Rafael Millan-Gabet,
Adam E. Rubinstein,
Michael L. Sitko,
Judit Sturmann,
Theo A. ten Brummelaar,
Yamina Touhami
Abstract:
The physical processes occurring within the inner few astronomical units of proto-planetary disks surrounding Herbig Ae stars are crucial to setting the environment in which the outer planet-forming disk evolves and put critical constraints on the processes of accretion and planet migration. We present the most complete published sample of high angular resolution H- and K-band observations of the…
▽ More
The physical processes occurring within the inner few astronomical units of proto-planetary disks surrounding Herbig Ae stars are crucial to setting the environment in which the outer planet-forming disk evolves and put critical constraints on the processes of accretion and planet migration. We present the most complete published sample of high angular resolution H- and K-band observations of the stars HD 163296 and HD 190073, including 30 previously unpublished nights of observations of the former and 45 nights of the latter with the CHARA long-baseline interferometer, in addition to archival VLTI data. We confirm previous observations suggesting significant near-infrared emission originates within the putative dust evaporation front of HD 163296 and show this is the case for HD 190073 as well. The H- and K-band sizes are the same within $(3 \pm 3)\%$ for HD 163296 and within $(6 \pm 10)\%$ for HD 190073. The radial surface brightness profiles for both disks are remarkably Gaussian-like with little or no sign of the sharp edge expected for a dust evaporation front. Coupled with spectral energy distribution analysis, our direct measurements of the stellar flux component at H and K bands suggest that HD 190073 is much younger (<400 kyr) and more massive (~5.6 M$_\odot$) than previously thought, mainly as a consequence of the new Gaia distance (891 pc).
△ Less
Submitted 9 November, 2018;
originally announced November 2018.
-
Detection of orbital motions near the last stable circular orbit of the massive black hole SgrA*
Authors:
GRAVITY Collaboration,
R. Abuter,
A. Amorim,
M. Bauböck,
J. P. Berger,
H. Bonnet,
W. Brandner,
Y. Clénet,
V. Coudé du Foresto,
P. T. de Zeeuw,
C. Deen,
J. Dexter,
G. Duvert,
A. Eckart,
F. Eisenhauer,
N. M. Förster Schreiber,
P. Garcia,
F. Gao,
E. Gendron,
R. Genzel,
S. Gillessen,
P. Guajardo,
M. Habibi,
X. Haubois,
Th. Henning
, et al. (35 additional authors not shown)
Abstract:
We report the detection of continuous positional and polarization changes of the compact source SgrA* in high states ('flares') of its variable near- infrared emission with the near-infrared GRAVITY-Very Large Telescope Interferometer (VLTI) beam-combining instrument. In three prominent bright flares, the position centroids exhibit clockwise looped motion on the sky, on scales of typically 150 mic…
▽ More
We report the detection of continuous positional and polarization changes of the compact source SgrA* in high states ('flares') of its variable near- infrared emission with the near-infrared GRAVITY-Very Large Telescope Interferometer (VLTI) beam-combining instrument. In three prominent bright flares, the position centroids exhibit clockwise looped motion on the sky, on scales of typically 150 micro-arcseconds over a few tens of minutes, corresponding to about 30% the speed of light. At the same time, the flares exhibit continuous rotation of the polarization angle, with about the same 45(+/-15)-minute period as that of the centroid motions. Modelling with relativistic ray tracing shows that these findings are all consistent with a near face-on, circular orbit of a compact polarized 'hot spot' of infrared synchrotron emission at approximately six to ten times the gravitational radius of a black hole of 4 million solar masses. This corresponds to the region just outside the innermost, stable, prograde circular orbit (ISCO) of a Schwarzschild-Kerr black hole, or near the retrograde ISCO of a highly spun-up Kerr hole. The polarization signature is consistent with orbital motion in a strong poloidal magnetic field.
△ Less
Submitted 30 October, 2018;
originally announced October 2018.
-
Multiple Star Systems in the Orion Nebula
Authors:
GRAVITY collaboration,
Martina Karl,
Oliver Pfuhl,
Frank Eisenhauer,
Reinhard Genzel,
Rebekka Grellmann,
Maryam Habibi,
Roberto Abuter,
Matteo Accardo,
António Amorim,
Narsireddy Anugu,
Gerardo Ávila,
Myriam Benisty,
Jean-Philippe Berger,
Nicolas Bland,
Henri Bonnet,
Pierre Bourget,
Wolfgang Brandner,
Roland Brast,
Alexander Buron,
Alessio Caratti o Garatti,
Frédéric Chapron,
Yann Clénet,
Claude Collin,
Vincent Coudé du Foresto
, et al. (111 additional authors not shown)
Abstract:
This work presents an interferometric study of the massive-binary fraction in the Orion Trapezium Cluster with the recently comissioned GRAVITY instrument. We observe a total of 16 stars of mainly OB spectral type. We find three previously unknown companions for $θ^1$ Ori B, $θ^2$ Ori B, and $θ^2$ Ori C. We determine a separation for the previously suspected companion of NU Ori. We confirm four co…
▽ More
This work presents an interferometric study of the massive-binary fraction in the Orion Trapezium Cluster with the recently comissioned GRAVITY instrument. We observe a total of 16 stars of mainly OB spectral type. We find three previously unknown companions for $θ^1$ Ori B, $θ^2$ Ori B, and $θ^2$ Ori C. We determine a separation for the previously suspected companion of NU Ori. We confirm four companions for $θ^1$ Ori A, $θ^1$ Ori C, $θ^1$ Ori D, and $θ^2$ Ori A, all with substantially improved astrometry and photometric mass estimates. We refine the orbit of the eccentric high-mass binary $θ^1$ Ori C and we are able to derive a new orbit for $θ^1$ Ori D. We find a system mass of 21.7 $M_{\odot}$ and a period of $53$ days. Together with other previously detected companions seen in spectroscopy or direct imaging, eleven of the 16 high-mass stars are multiple systems. We obtain a total number of 22 companions with separations up to 600 AU. The companion fraction of the early B and O stars in our sample is about 2, significantly higher than in earlier studies of mostly OB associations. The separation distribution hints towards a bimodality. Such a bimodality has been previously found in A stars, but rarely in OB binaries, which up to this point have been assumed to be mostly compact with a tail of wider companions. We also do not find a substantial population of equal-mass binaries. The observed distribution of mass ratios declines steeply with mass, and like the direct star counts, indicates that our companions follow a standard power law initial mass function. Again, this is in contrast to earlier findings of flat mass ratio distributions in OB associations. We exclude collision as a dominant formation mechanism but find no clear preference for core accretion or competitive accretion.
△ Less
Submitted 27 September, 2018;
originally announced September 2018.
-
GRAVITY chromatic imaging of Eta Car's core
Authors:
GRAVITY Collaboration,
J. Sanchez-Bermudez,
G. Weigelt,
J. M. Bestenlehner,
P. Kervella,
W. Brandner,
Th. Henning,
A. Müller,
G. Perrin,
J. -U. Pott,
M. Schöller,
R. van Boekel,
R. Abuter,
M. Accardo,
A. Amorim,
N. Anugu,
G. Ávila,
M. Benisty,
J. P. Berger,
N. Blind,
H. Bonnet,
P. Bourget,
R. Brast,
A. Buron,
F. Cantalloube
, et al. (110 additional authors not shown)
Abstract:
Eta Car is one of the most intriguing luminous blue variables in the Galaxy. Observations and models at different wavelengths suggest a central binary with a 5.54 yr period residing in its core. 2D and 3D radiative transfer and hydrodynamic simulations predict a primary with a dense and slow stellar wind that interacts with the faster and lower density wind of the secondary. The wind-wind collisio…
▽ More
Eta Car is one of the most intriguing luminous blue variables in the Galaxy. Observations and models at different wavelengths suggest a central binary with a 5.54 yr period residing in its core. 2D and 3D radiative transfer and hydrodynamic simulations predict a primary with a dense and slow stellar wind that interacts with the faster and lower density wind of the secondary. The wind-wind collision scenario suggests that the secondary's wind penetrates the primary's wind creating a low-density cavity in it, with dense walls where the two winds interact. We aim to trace the inner ~5-50 au structure of Eta Car's wind-wind interaction, as seen through BrG and, for the first time, through the He I 2s-2p line. We have used spectro-interferometric observations with GRAVITY at the VLTI. Our modeling of the continuum allows us to estimate its FWHM angular size close to 2 mas and an elongation ratio of 1.06 +/- 0.05 over a PA = 130 +/- 20 deg. Our CMFGEN modeling helped us to confirm that the role of the secondary should be taken into account to properly reproduce the observed BrG and He I lines. Chromatic images across BrG reveal a southeast arc-like feature, possibly associated to the hot post-shocked winds flowing along the cavity wall. The images of He I 2s-2p served to constrain the 20 mas structure of the line-emitting region. The observed morphology of He I suggests that the secondary is responsible for the ionized material that produces the line profile. Both the BrG and the He I 2s-2p maps are consistent with previous hydrodynamical models of the colliding wind scenario. Future dedicated simulations together with an extensive interferometric campaign are necessary to refine our constraints on the wind and stellar parameters of the binary, which finally will help us predict the evolutionary path of Eta Car.
△ Less
Submitted 6 August, 2018;
originally announced August 2018.
-
The Planet Formation Imager
Authors:
John D. Monnier,
Stefan Kraus,
Michael J. Ireland,
Fabien Baron,
Amelia Bayo,
Jean-Philippe Berger,
Michelle Creech-Eakman,
Ruobing Dong,
Gaspard Duchene,
Catherine Espaillat,
Chris Haniff,
Sebastian Honig,
Andrea Isella,
Attila Juhasz,
Lucas Labadie,
Sylvestre Lacour,
Stephanie Leifer,
Antoine Merand,
Ernest Michael,
Stefano Minardi,
Christoph Mordasini,
David Mozurkewich,
Johan Olofsson,
Claudia Paladini,
Romain Petrov
, et al. (14 additional authors not shown)
Abstract:
The Planet Formation Imager (PFI, www.planetformationimager.org) is a next-generation infrared interferometer array with the primary goal of imaging the active phases of planet formation in nearby star forming regions. PFI will be sensitive to warm dust emission using mid-infrared capabilities made possible by precise fringe tracking in the near-infrared. An L/M band combiner will be especially se…
▽ More
The Planet Formation Imager (PFI, www.planetformationimager.org) is a next-generation infrared interferometer array with the primary goal of imaging the active phases of planet formation in nearby star forming regions. PFI will be sensitive to warm dust emission using mid-infrared capabilities made possible by precise fringe tracking in the near-infrared. An L/M band combiner will be especially sensitive to thermal emission from young exoplanets (and their disks) with a high spectral resolution mode to probe the kinematics of CO and H2O gas. In this paper, we give an overview of the main science goals of PFI, define a baseline PFI architecture that can achieve those goals, point at remaining technical challenges, and suggest activities today that will help make the Planet Formation Imager facility a reality.
△ Less
Submitted 30 July, 2018;
originally announced July 2018.
-
Detection of the gravitational redshift in the orbit of the star S2 near the Galactic centre massive black hole
Authors:
GRAVITY Collaboration,
R. Abuter,
A. Amorim,
N. Anugu,
M. Bauböck,
M. Benisty,
J. P. Berger,
N. Blind,
H. Bonnet,
W. Brandner,
A. Buron,
C. Collin,
F. Chapron,
Y. Clénet,
V. Coudé du Foresto,
P. T. de Zeeuw,
C. Deen,
F. Delplancke-Ströbele,
R. Dembet,
J. Dexter,
G. Duvert,
A. Eckart,
F. Eisenhauer,
G. Finger,
N. M. Förster Schreiber
, et al. (73 additional authors not shown)
Abstract:
The highly elliptical, 16-year-period orbit of the star S2 around the massive black hole candidate Sgr A* is a sensitive probe of the gravitational field in the Galactic centre. Near pericentre at 120 AU, ~1400 Schwarzschild radii, the star has an orbital speed of ~7650 km/s, such that the first-order effects of Special and General Relativity have now become detectable with current capabilities. O…
▽ More
The highly elliptical, 16-year-period orbit of the star S2 around the massive black hole candidate Sgr A* is a sensitive probe of the gravitational field in the Galactic centre. Near pericentre at 120 AU, ~1400 Schwarzschild radii, the star has an orbital speed of ~7650 km/s, such that the first-order effects of Special and General Relativity have now become detectable with current capabilities. Over the past 26 years, we have monitored the radial velocity and motion on the sky of S2, mainly with the SINFONI and NACO adaptive optics instruments on the ESO Very Large Telescope, and since 2016 and leading up to the pericentre approach in May 2018, with the four-telescope interferometric beam-combiner instrument GRAVITY. From data up to and including pericentre, we robustly detect the combined gravitational redshift and relativistic transverse Doppler effect for S2 of z ~ 200 km/s / c with different statistical analysis methods. When parameterising the post-Newtonian contribution from these effects by a factor f, with f = 0 and f = 1 corresponding to the Newtonian and general relativistic limits, respectively, we find from posterior fitting with different weighting schemes f = 0.90 +/- 0.09 (stat) +\- 0.15 (sys). The S2 data are inconsistent with pure Newtonian dynamics.
△ Less
Submitted 24 July, 2018;
originally announced July 2018.
-
The Median Probability Model and Correlated Variables
Authors:
Marilena Barbieri,
James O. Berger,
Edward I. George,
Veronika Rockova
Abstract:
The median probability model (MPM) Barbieri and Berger (2004) is defined as the model consisting of those variables whose marginal posterior probability of inclusion is at least 0.5. The MPM rule yields the best single model for prediction in orthogonal and nested correlated designs. This result was originally conceived under a specific class of priors, such as the point mass mixtures of non-infor…
▽ More
The median probability model (MPM) Barbieri and Berger (2004) is defined as the model consisting of those variables whose marginal posterior probability of inclusion is at least 0.5. The MPM rule yields the best single model for prediction in orthogonal and nested correlated designs. This result was originally conceived under a specific class of priors, such as the point mass mixtures of non-informative and g-type priors. The MPM rule, however, has become so very popular that it is now being deployed for a wider variety of priors and under correlated designs, where the properties of MPM are not yet completely understood. The main thrust of this work is to shed light on properties of MPM in these contexts by (a) characterizing situations when MPM is still safe under correlated designs, (b) providing significant generalizations of MPM to a broader class of priors (such as continuous spike-and-slab priors). We also provide new supporting evidence for the suitability of g-priors, as opposed to independent product priors, using new predictive matching arguments. Furthermore, we emphasize the importance of prior model probabilities and highlight the merits of non-uniform prior probability assignments using the notion of model aggregates.
△ Less
Submitted 17 August, 2018; v1 submitted 22 July, 2018;
originally announced July 2018.
-
Unphysical Discontinuities in GW Methods
Authors:
Mickaël Véril,
Pina Romaniello,
J. A. Berger,
Pierre-François Loos
Abstract:
We report unphysical irregularities and discontinuities in some key experimentally-measurable quantities computed within the GW approximation of many-body perturbation theory applied to molecular systems. In particular, we show that the solution obtained with partially self-consistent GW schemes depends on the algorithm one uses to solve self-consistently the quasi-particle (QP) equation. The main…
▽ More
We report unphysical irregularities and discontinuities in some key experimentally-measurable quantities computed within the GW approximation of many-body perturbation theory applied to molecular systems. In particular, we show that the solution obtained with partially self-consistent GW schemes depends on the algorithm one uses to solve self-consistently the quasi-particle (QP) equation. The main observation of the present study is that each branch of the self-energy is associated with a distinct QP solution, and that each switch between solutions implies a significant discontinuity in the quasiparticle energy as a function of the internuclear distance. Moreover, we clearly observe "ripple" effects, i.e., a discontinuity in one of the QP energies induces (smaller) discontinuities in the other QP energies. Going from one branch to another implies a transfer of weight between two solutions of the QP equation. The case of occupied, virtual and frontier orbitals are separately discussed on distinct diatomics. In particular, we show that multisolution behavior in frontier orbitals is more likely if the HOMO-LUMO gap is small.
△ Less
Submitted 7 September, 2018; v1 submitted 19 July, 2018;
originally announced July 2018.
-
Hi-5: a potential high-contrast thermal near-infrared imager for the VLTI
Authors:
D. Defrère,
M. Ireland,
O. Absil,
J. -P. Berger,
W. C. Danchi,
S. Ertel,
A. Gallenne,
F. Hénault,
P. Hinz,
E. Huby,
S. Kraus,
L. Labadie,
J. -B. Le Bouquin,
G. Martin,
A. Matter,
B. Mennesson,
A. Mérand,
S. Minardi,
J. D. Monnier,
B. Norris,
G. Orban de Xivry,
E. Pedretti,
J. -U. Pott,
M. Reggiani,
E. Serabyn
, et al. (3 additional authors not shown)
Abstract:
Hi-5 is a high-contrast (or high dynamic range) infrared imager project for the VLTI. Its main goal is to characterize young extra-solar planetary systems and exozodiacal dust around southern main-sequence stars. In this paper, we present an update of the project and key technology pathways to improve the contrast achieved by the VLTI. In particular, we discuss the possibility to use integrated op…
▽ More
Hi-5 is a high-contrast (or high dynamic range) infrared imager project for the VLTI. Its main goal is to characterize young extra-solar planetary systems and exozodiacal dust around southern main-sequence stars. In this paper, we present an update of the project and key technology pathways to improve the contrast achieved by the VLTI. In particular, we discuss the possibility to use integrated optics, proven in the near-infrared, in the thermal near-infrared (L and M bands, 3-5~$μ$m) and advanced fringe tracking strategies. We also address the strong exoplanet science case (young exoplanets, planet formation, and exozodiacal disks) offered by this wavelength regime as well as other possible science cases such as stellar physics (fundamental parameters and multiplicity) and extragalactic astrophysics (active galactic nuclei and fundamental constants). Synergies and scientific preparation for other potential future instruments such as the Planet Formation Imager are also briefly discussed.
△ Less
Submitted 19 July, 2018;
originally announced July 2018.
-
Characterisation of ALPAO deformable mirrors for the NAOMI VLTI Auxiliary Telescopes Adaptive Optics
Authors:
Jean-Baptiste Le Bouquin,
Jean-Philippe Berger,
Jean-Luc Beuzit,
Eric Cottalorda,
Alain Delboulbe,
Sebastien E. Egner,
Frederic Yves Joseph Gonte,
Sylvain Guieu,
Pierre Haguenauer,
Laurent Jocou,
Yves Magnard,
Thibaut Moulin,
Sylvain Rochat,
Christophe Verinaud,
Julien Woillez
Abstract:
The Very Large Telescope Interferometer Auxiliary Telescopes will soon be equipped with an adaptive optics system called NAOMI. The corrective optics deformable mirror is the commercial DM241 from ALPAO. Being part of an interferometer operating from visible to mid-infrared, the DMs of NAOMI face several challenges (high level of reliability, open-loop chopping, piston-free control, WFS/DM pupil r…
▽ More
The Very Large Telescope Interferometer Auxiliary Telescopes will soon be equipped with an adaptive optics system called NAOMI. The corrective optics deformable mirror is the commercial DM241 from ALPAO. Being part of an interferometer operating from visible to mid-infrared, the DMs of NAOMI face several challenges (high level of reliability, open-loop chopping, piston-free control, WFS/DM pupil rotation, high desired bandwidth and stroke). We here describe our extensive characterization of the DMs through measurements and simulations. We summarize the operational scenario we have defined to handle the specific mirror properties. We conclude that the ALPAO DMs have overall excellent properties that fulfill most of the stringent requirements and that deviations from specifications are easily handled. To our knowledge, NAOMI will be the first astronomical system with a command in true Zernike modes (allowing software rotation), and the first astronomical system in which a chopping is performed with the deformable mirror (5'' sky, at 5~Hz).
△ Less
Submitted 27 June, 2018;
originally announced June 2018.
-
Dynamical masses of M-dwarf binaries in young moving groups: I - The case of TWA 22 and GJ 2060
Authors:
L. Rodet,
M. Bonnefoy,
S. Durkan,
H. Beust,
A-M Lagrange,
J. E. Schlieder,
M. Janson,
A. Grandjean,
G. Chauvin,
S. Messina,
A. -L. Maire,
W. Brandner,
J. Girard,
P. Delorme,
B. Biller,
C. Bergfors,
S. Lacour,
M. Feldt,
T. Henning,
A. Boccaletti,
J. -B. Le Bouquin,
J. -P. Berger,
J. -L. Monin,
S. Udry,
S. Peretti
, et al. (26 additional authors not shown)
Abstract:
Evolutionary models are widely used to infer the mass of stars, brown dwarfs, and giant planets. Their predictions are thought to be less reliable at young ages ($<$ 200 Myr) and in the low-mass regime ($\mathrm{<1~M_{\odot}}$). GJ 2060 AB and TWA 22 AB are two rare astrometric M-dwarf binaries respectively members of the AB Doradus and Beta Pictoris moving groups. As their dynamical mass can be m…
▽ More
Evolutionary models are widely used to infer the mass of stars, brown dwarfs, and giant planets. Their predictions are thought to be less reliable at young ages ($<$ 200 Myr) and in the low-mass regime ($\mathrm{<1~M_{\odot}}$). GJ 2060 AB and TWA 22 AB are two rare astrometric M-dwarf binaries respectively members of the AB Doradus and Beta Pictoris moving groups. As their dynamical mass can be measured within a few years, they can be used to calibrate the evolutionary tracks and set new constraints on the age of young moving groups. We find a total mass of $\mathrm{0.18\pm 0.02~M_\odot}$ for TWA 22. That mass is in good agreement with model predictions at the age of the Beta Pic moving group. We obtain a total mass of $\mathrm{1.09 \pm 0.10~M_{\odot}}$ for GJ 2060. We estimate a spectral type of M$1\pm0.5$, $\mathrm{L/L_{\odot}=-1.20\pm0.05}$ dex, and $\mathrm{T_{eff}=3700\pm100}$ K for GJ 2060 A. The B component is a M$3\pm0.5$ dwarf with $\mathrm{L/L_{\odot}=-1.63\pm0.05}$ dex and $\mathrm{T_{eff}=3400\pm100}$ K. The dynamical mass of GJ 2060 AB is inconsistent with the most recent models predictions (BCAH15, PARSEC) for an ABDor age in the range 50-150 Myr. It is 10 to 20\% (1-2 sigma, depending on the assumed age) above the models predictions, corresponding to an underestimation of $0.10$ to $0.20~\mathrm{M_\odot}$. Coevality suggests a young age for the system ($\sim$ 50 Myr) according to most evolutionary models. TWA 22 validates the predictions of recent evolutionary tracks at $\sim$20 Myr. On the other hand, we evidence a 1-2 sigma mismatch between the predicted and observed mass of GJ 2060 AB. This slight departure may indicate that one of the star hosts a tight companion. Alternatively, this would confirm the models tendency to underestimate the mass of young low-mass stars.
△ Less
Submitted 14 June, 2018;
originally announced June 2018.
-
On the solution of coupled heat and moisture transport in porous material
Authors:
Julien Berger,
Suelen Gasparin,
Denys Dutykh,
Nathan Mendes
Abstract:
Comparisons of experimental observation of heat and moisture transfer through porous building materials with numerical results have been presented in numerous studies reported in literature. However, some discrepancies have been observed, highlighting underestimation of sorption process and overestimation of desorption process. Some studies intend to explain the discrepancies by analysing the impo…
▽ More
Comparisons of experimental observation of heat and moisture transfer through porous building materials with numerical results have been presented in numerous studies reported in literature. However, some discrepancies have been observed, highlighting underestimation of sorption process and overestimation of desorption process. Some studies intend to explain the discrepancies by analysing the importance of hysteresis effects as well as carrying out sensitivity analyses on the input parameters as convective transfer coefficients. This article intends to investigate the accuracy and efficiency of the coupled solution by adding advective transfer of both heat and moisture in the physical model. In addition, the efficient Scharfetter and Gummel numerical scheme is proposed to solve the system of advection-diffusion equations, which has the advantages of being well-balanced and asymptotically preserving. Moreover, the scheme is particularly efficient in terms of accuracy and reduction of computational time when using large spatial discretisation parameters. Several linear and non-linear cases are studied to validate the method and highlight its specific features. At the end, an experimental benchmark from the literature is considered. The numerical results are compared to the experimental data for a pure diffusive model and also for the proposed model. The latter presents better agreement with the experimental data. The influence of the hysteresis effects on the moisture capacity is also studied, by adding a third differential equation.
△ Less
Submitted 14 May, 2018; v1 submitted 20 April, 2018;
originally announced May 2018.
-
Perfect tangles
Authors:
Johannes Berger,
Tobias J. Osborne
Abstract:
We introduce perfect tangles for modular tensor categories. These are intended to generalise the perfect tensors first introduced in the context of a toy model for the AdS/CFT correspondence. We construct perfect tangles for several categories of relevance for topological quantum computation, including the Temperley-Lieb, Fibonacci, Kuperberg spider, and Haagerup planar algebras. A general inducti…
▽ More
We introduce perfect tangles for modular tensor categories. These are intended to generalise the perfect tensors first introduced in the context of a toy model for the AdS/CFT correspondence. We construct perfect tangles for several categories of relevance for topological quantum computation, including the Temperley-Lieb, Fibonacci, Kuperberg spider, and Haagerup planar algebras. A general inductive construction proposed by Vaughan Jones for perfect tangles is also described.
△ Less
Submitted 9 April, 2018;
originally announced April 2018.
-
Analysis and improvement of the VTT mold growth model: application to bamboo fiberboard
Authors:
Julien Berger,
Hervé Le Meur,
Denys Dutykh,
Dang Mao Nguyen,
A. -C. Grillet
Abstract:
The reliability of a model is its accuracy in predicting the physical phenomena using the known input parameters. It also depends on the model's ability to estimate relevant parameters using observations of the physical phenomena. In this paper, the reliability of the VTT model is investigated under these two criteria for various given temperature and relative humidity constant in time. First of a…
▽ More
The reliability of a model is its accuracy in predicting the physical phenomena using the known input parameters. It also depends on the model's ability to estimate relevant parameters using observations of the physical phenomena. In this paper, the reliability of the VTT model is investigated under these two criteria for various given temperature and relative humidity constant in time. First of all, experiments are conducted on bamboo fiberboard. Using these data, five parameters of the VTT model, defining the mold vulnerability class of a material, are identified. The results highlight that the determined parameters are not within the range of the classes defined in the VTT model. In addition, the quality of the parameter estimation is not satisfactory. Then the sensitivity of the numerical results of the VTT model is analyzed by varying an input parameter. These investigations show that the VTT mathematical formulation of the physical model of mold growth is not reliable. An improved model is proposed with a new mathematical formulation. It is inspired by the logistic equation whose parameters are estimated using the experimental data obtained. The parameter estimation is very satisfactory. In the last parts of the paper, the numerical predictions of the improved model are compared to experimental data from the literature to prove its reliability.
△ Less
Submitted 10 April, 2018; v1 submitted 5 April, 2018;
originally announced April 2018.
-
Green functions and self-consistency: insights from the spherium model
Authors:
Pierre-François Loos,
Pina Romaniello,
J. A. Berger
Abstract:
We report an exhaustive study of the performance of different variants of Green function methods for the spherium model in which two electrons are confined to the surface of a sphere and interact via a genuine long-range Coulomb operator. We show that the spherium model provides a unique paradigm to study electronic correlation effects from the weakly correlated regime to the strongly correlated r…
▽ More
We report an exhaustive study of the performance of different variants of Green function methods for the spherium model in which two electrons are confined to the surface of a sphere and interact via a genuine long-range Coulomb operator. We show that the spherium model provides a unique paradigm to study electronic correlation effects from the weakly correlated regime to the strongly correlated regime, since the mathematics are simple while the physics is rich. We compare perturbative GW, partially self-consistent GW and second-order Green function (GF2) methods for the computation of ionization potentials, electron affinities, energy gaps, correlation energies as well as singlet and triplet neutral excitations by solving the Bethe-Salpeter equation (BSE). We discuss the problem of self-screening in GW and show that it can be partially solved with a second-order screened exchange correction (SOSEX). We find that, in general, self-consistency deteriorates the results with respect to those obtained within perturbative approaches with a Hartree-Fock starting point. Finally, we unveil an important problem of partial self-consistency in GW: in the weakly correlated regime, it can produce artificial discontinuities in the self-energy caused by satellite resonances with large weights.
△ Less
Submitted 12 March, 2018;
originally announced March 2018.
-
The path towards high-contrast imaging with the VLTI: the Hi-5 project
Authors:
D. Defrère,
O. Absil,
J. -P. Berger,
T. Boulet,
W. C. Danchi,
S. Ertel,
A. Gallenne,
F. Hénault,
P. Hinz,
E. Huby,
M. Ireland,
S. Kraus,
L. Labadie,
J. -B. Le Bouquin,
G. Martin,
A. Matter,
A. Mérand,
B. Mennesson,
S. Minardi,
J. Monnier,
B. Norris,
G. Orban de Xivry,
E. Pedretti,
J. -U. Pott,
M. Reggiani
, et al. (4 additional authors not shown)
Abstract:
The development of high-contrast capabilities has long been recognized as one of the top priorities for the VLTI. As of today, the VLTI routinely achieves contrasts of a few 10$^{-3}$ in the near-infrared with PIONIER (H band) and GRAVITY (K band). Nulling interferometers in the northern hemisphere and non-redundant aperture masking experiments have, however, demonstrated that contrasts of at leas…
▽ More
The development of high-contrast capabilities has long been recognized as one of the top priorities for the VLTI. As of today, the VLTI routinely achieves contrasts of a few 10$^{-3}$ in the near-infrared with PIONIER (H band) and GRAVITY (K band). Nulling interferometers in the northern hemisphere and non-redundant aperture masking experiments have, however, demonstrated that contrasts of at least a few 10$^{-4}$ are within reach using specific beam combination and data acquisition techniques. In this paper, we explore the possibility to reach similar or higher contrasts on the VLTI. After reviewing the state-of-the-art in high-contrast infrared interferometry, we discuss key features that made the success of other high-contrast interferometric instruments (e.g., integrated optics, nulling, closure phase, and statistical data reduction) and address possible avenues to improve the contrast of the VLTI by at least one order of magnitude. In particular, we discuss the possibility to use integrated optics, proven in the near-infrared, in the thermal near-infrared (L and M bands, 3-5 $μ$m), a sweet spot to image and characterize young extra-solar planetary systems. Finally, we address the science cases of a high-contrast VLTI imaging instrument and focus particularly on exoplanet science (young exoplanets, planet formation, and exozodiacal disks), stellar physics (fundamental parameters and multiplicity), and extragalactic astrophysics (active galactic nuclei and fundamental constants). Synergies and scientific preparation for other potential future instruments such as the Planet Formation Imager are also briefly discussed.
△ Less
Submitted 25 May, 2018; v1 submitted 12 January, 2018;
originally announced January 2018.
-
RobustGaSP: Robust Gaussian Stochastic Process Emulation in R
Authors:
Mengyang Gu,
Jesús Palomo,
James O. Berger
Abstract:
Gaussian stochastic process emulation is a powerful tool for approximating computationally intensive computer models. However, estimation of parameters in the GaSP emulator is a challenging task. No closed-form estimator is available and many numerical problems arise with standard estimates, e.g., the maximum likelihood estimator. In this package, we implement a marginal posterior mode estimator,…
▽ More
Gaussian stochastic process emulation is a powerful tool for approximating computationally intensive computer models. However, estimation of parameters in the GaSP emulator is a challenging task. No closed-form estimator is available and many numerical problems arise with standard estimates, e.g., the maximum likelihood estimator. In this package, we implement a marginal posterior mode estimator, for special priors and parameterizations, an estimation method that meets the robust parameter estimation criteria discussed in \cite{Gu2018robustness}; mathematical reasons are provided therein to explain why robust parameter estimation can greatly improve predictive performance of the emulator. In addition, inert inputs (inputs that almost have no effect on the variability of a function) can be identified from the marginal posterior mode estimation, at no extra computational cost. The package also implements the parallel partial Gaussian stochastic process (PP GaSP) emulator (\cite{gu2016parallel}) for the scenario where the computer model has multiple outputs on e.g. spatial-temporal coordinates. The package can be operated in a default mode, but also allows numerous user specifications, such as the capability of specifying trend functions and noise terms. Examples are studied herein to highlight the performance of the package in terms of out-of-sample prediction.}
△ Less
Submitted 14 June, 2019; v1 submitted 5 January, 2018;
originally announced January 2018.
-
A note on recent criticisms to Birnbaum's theorem
Authors:
Víctor Peña,
James O. Berger
Abstract:
In this note, we provide critical commentary on two articles that cast doubt on the validity and implications of Birnbaum's theorem: Evans (2013) and Mayo (2014). In our view, the proof is correct and the consequences of the theorem are alive and well.
In this note, we provide critical commentary on two articles that cast doubt on the validity and implications of Birnbaum's theorem: Evans (2013) and Mayo (2014). In our view, the proof is correct and the consequences of the theorem are alive and well.
△ Less
Submitted 21 November, 2017;
originally announced November 2017.
-
Restricted type II maximum likelihood priors on regression coefficients
Authors:
Víctor Peña,
James O. Berger
Abstract:
In Bayesian hypothesis testing and model selection, prior distributions must be chosen carefully. For example, setting arbitrarily large prior scales for location parameters, which is common practice in estimation problems, can lead to undesirable behavior in testing (Lindley's paradox). We study the properties of some restricted type II maximum likelihood (type II ML) priors on regression coeffic…
▽ More
In Bayesian hypothesis testing and model selection, prior distributions must be chosen carefully. For example, setting arbitrarily large prior scales for location parameters, which is common practice in estimation problems, can lead to undesirable behavior in testing (Lindley's paradox). We study the properties of some restricted type II maximum likelihood (type II ML) priors on regression coefficients. In type II ML, hyperparameters are "estimated" by maximizing the marginal likelihood of a model. In this article, we define priors by estimating their variances or covariance matrices, adding restrictions which ensure that the resulting priors are at least as vague as conventional proper priors for model uncertainty. We find that these type II ML priors typically yield results that are close to answers obtained with the Bayesian Information Criterion (BIC).
△ Less
Submitted 21 November, 2019; v1 submitted 21 November, 2017;
originally announced November 2017.
-
Determination of spin Hall effect and spin diffusion length of Pt from self-consistent fitting of damping enhancement and inverse spin-orbit torque measurements
Authors:
Andrew J. Berger,
Eric R. J. Edwards,
Hans T. Nembach,
Olof Karis,
Mathias Weiler,
T. J. Silva
Abstract:
Understanding the evolution of spin-orbit torque (SOT) with increasing heavy-metal thickness in ferromagnet/normal metal (FM/NM) bilayers is critical for the development of magnetic memory based on SOT. However, several experiments have revealed an apparent discrepancy between damping enhancement and damping-like SOT regarding their dependence on NM thickness. Here, using linewidth and phase-resol…
▽ More
Understanding the evolution of spin-orbit torque (SOT) with increasing heavy-metal thickness in ferromagnet/normal metal (FM/NM) bilayers is critical for the development of magnetic memory based on SOT. However, several experiments have revealed an apparent discrepancy between damping enhancement and damping-like SOT regarding their dependence on NM thickness. Here, using linewidth and phase-resolved amplitude analysis of vector network analyzer ferromagnetic resonance (VNA-FMR) measurements, we simultaneously extract damping enhancement and both field-like and damping-like inverse SOT in Ni$_{80}$Fe$_{20}$/Pt bilayers as a function of Pt thickness. By enforcing an interpretation of the data which satisfies Onsager reciprocity, we find that both the damping enhancement and damping-like inverse SOT can be described by a single spin diffusion length ($\approx$ 4 nm), and that we can separate the spin pumping and spin memory loss (SML) contributions to the total damping. This analysis indicates that less than 40% of the angular momentum pumped by FMR through the Ni$_{80}$Fe$_{20}$/Pt interface is transported as spin current into the Pt. On account of the SML and corresponding reduction in total spin current available for spin-charge transduction in the Pt, we determine the Pt spin Hall conductivity ($σ_\mathrm{SH} = (2.36 \pm 0.04)\times10^6 Ω^{-1} \mathrm{m}^{-1}$) and bulk spin Hall angle ($θ_\mathrm{SH}=0.387 \pm0.008$) to be larger than commonly-cited values. These results suggest that Pt can be an extremely useful source of SOT if the FM/NM interface can be engineered to minimize SML. Lastly, we find that self-consistent fitting of the damping and SOT data is best achieved by a model with Elliott-Yafet spin relaxation and extrinsic inverse spin Hall effect, such that both the spin diffusion length and spin Hall conductivity are proportional to the Pt charge conductivity.
△ Less
Submitted 21 March, 2018; v1 submitted 21 November, 2017;
originally announced November 2017.
-
On the Ubiquity of Information Inconsistency for Conjugate Priors
Authors:
Joris Mulder,
James O. Berger,
Víctor Peña,
M. J. Bayarri
Abstract:
Informally, "Information Inconsistency" is the property that has been observed in many Bayesian hypothesis testing and model selection procedures whereby the Bayesian conclusion does not become definitive when the data seems to become definitive. An example is that, when performing a t-test using standard conjugate priors, the Bayes factor of the alternative hypothesis to the null hypothesis remai…
▽ More
Informally, "Information Inconsistency" is the property that has been observed in many Bayesian hypothesis testing and model selection procedures whereby the Bayesian conclusion does not become definitive when the data seems to become definitive. An example is that, when performing a t-test using standard conjugate priors, the Bayes factor of the alternative hypothesis to the null hypothesis remains bounded as the t statistic grows to infinity. This paper shows that information inconsistency is ubiquitous in Bayesian hypothesis testing under conjugate priors. Yet the title does not fully describe the paper, since we also show that theoretically recommended priors, including scale mixtures of conjugate priors and adaptive priors, are information consistent. Hence the paper is simply a forceful warning that use of conjugate priors in testing and model selection is highly problematical, and should be replaced by the information consistent alternatives.
△ Less
Submitted 26 October, 2017;
originally announced October 2017.
-
Revisiting the origin of satellites in core level photoemission of transparent conducting oxides: the case of $n$-doped SnO$_2$
Authors:
F. Borgatti,
J. A. Berger,
D. Céolin,
J. S. Zhou,
J. J. Kas,
M. Guzzo,
C. F. McConville,
F. Offi,
G. Panaccione,
A. Regoutz,
D. J. Payne,
J. -P. Rueff,
O. Bierwagen,
M. E. White,
J. S. Speck,
M. Gatti,
R. G. Egdell
Abstract:
The longstanding problem of interpretation of satellite structures in core level photoemission spectra of metallic systems with a low density of conduction electrons is addressed using the specific example of Sb-doped SnO$_2$. Comparison of {\it ab initio} many-body calculations with experimental hard X-ray photoemission spectra of the Sn 4$d$ states shows that strong satellites are produced by co…
▽ More
The longstanding problem of interpretation of satellite structures in core level photoemission spectra of metallic systems with a low density of conduction electrons is addressed using the specific example of Sb-doped SnO$_2$. Comparison of {\it ab initio} many-body calculations with experimental hard X-ray photoemission spectra of the Sn 4$d$ states shows that strong satellites are produced by coupling of the Sn core hole to the plasma oscillations of the free electrons introduced by doping. Within the same theoretical framework, spectral changes of the valence band spectra are also related to dynamical screening effects. These results demonstrate that, for the interpretation of electron correlation features in the core level photoelectron spectra of such narrow-band materials, going beyond the homogeneous electron gas electron-plasmon coupling model is essential.
△ Less
Submitted 24 October, 2017;
originally announced October 2017.
-
Many-body perturbation theory and non-perturbative approaches: the screened interaction as key ingredient
Authors:
Walter Tarantino,
Bernardo S. Mendoza,
Pina Romaniello,
J. A. Berger,
Lucia Reining
Abstract:
Many-body perturbation theory is often formulated in terms of an expansion in the dressed instead of the bare Green's function, and in the screened instead of the bare Coulomb interaction. However, screening can be calculated on different levels of approximation, and it is important to define what is the most appropriate choice. We explore this question by studying a zero-dimensional model (so cal…
▽ More
Many-body perturbation theory is often formulated in terms of an expansion in the dressed instead of the bare Green's function, and in the screened instead of the bare Coulomb interaction. However, screening can be calculated on different levels of approximation, and it is important to define what is the most appropriate choice. We explore this question by studying a zero-dimensional model (so called 'one-point model') that retains the structure of the full equations. We study both linear and non-linear response approximations to the screening. We find that an expansion in terms of the screening in the random phase approximation is the most promising way for an application in real systems. Moreover, by making use of the nonperturbative features of the Kadanoff-Baym equation for the one-body Green's function, we obtain an approximate solution in our model that is very promising, although its applicability to real systems has still to be explored.
△ Less
Submitted 18 October, 2017;
originally announced October 2017.
-
Variable dynamics in the inner disk of HD 135344B revealed with multi-epoch scattered light imaging
Authors:
Tomas Stolker,
Mike Sitko,
Bernard Lazareff,
Myriam Benisty,
Carsten Dominik,
Rens Waters,
Michiel Min,
Sebastian Perez,
Julien Milli,
Antonio Garufi,
Jozua de Boer,
Christian Ginski,
Stefan Kraus,
Jean-Philippe Berger,
Henning Avenhaus
Abstract:
We present multi-epoch VLT/SPHERE observations of the protoplanetary disk around HD 135344B (SAO 206462). The $J$-band scattered light imagery reveal, with high spatial resolution ($\sim$41 mas, 6.4 au), the disk surface beyond $\sim$20 au. Temporal variations are identified in the azimuthal brightness distributions of all epochs, presumably related to the asymmetrically shading dust distribution…
▽ More
We present multi-epoch VLT/SPHERE observations of the protoplanetary disk around HD 135344B (SAO 206462). The $J$-band scattered light imagery reveal, with high spatial resolution ($\sim$41 mas, 6.4 au), the disk surface beyond $\sim$20 au. Temporal variations are identified in the azimuthal brightness distributions of all epochs, presumably related to the asymmetrically shading dust distribution in the inner disk. These shadows manifest themselves as narrow lanes, cast by localized density enhancements, and broader features which possibly trace the larger scale dynamics of the inner disk. We acquired visible and near-infrared photometry which shows variations up to 10% in the $JHK$ bands, possibly correlated with the presence of the shadows. Analysis of archival VLTI/PIONIER $H$-band visibilities constrain the orientation of the inner disk to $i = 18.2°^{+3.4}_{-4.1}$ and ${\rm PA} = 57.3°\pm 5.7°$, consistent with an alignment with the outer disk or a minor disk warp of several degrees. The latter scenario could explain the broad, quasi-stationary shadowing in N-NW direction in case the inclination of the outer disk is slightly larger. The correlation between the shadowing and the near-infrared excess is quantified with a grid of radiative transfer models. The variability of the scattered light contrast requires extended variations in the inner disk atmosphere ($H/r \lesssim 0.2$). Possible mechanisms that may cause asymmetric variations in the optical depth ($Δτ\lesssim1$) through the atmosphere of the inner disk include turbulent fluctuations, planetesimal collisions, or a dusty disk wind, possibly enhanced by a minor disk warp. A fine temporal sampling is required to follow day-to-day changes of the shadow patterns which may be a face-on variant of the UX Orionis phenomenon.
△ Less
Submitted 6 October, 2017;
originally announced October 2017.
-
The wind and the magnetospheric accretion onto the T Tauri star S Coronae Australis at sub-au resolution
Authors:
R. Garcia Lopez,
K. Perraut,
A. Caratti o Garatti,
B. Lazareff,
J. Sanchez-Bermudez,
M. Benisty,
C. Dougados,
L. Labadie,
W. Brandner,
P. J. V. Garcia,
Th. Henning,
T. P. Ray,
R. Abuter,
A. Amorim,
N. Anugu,
J. P. Berger,
H. Bonnet,
A. Buron,
P. Caselli,
Y. Clénet,
V. Coudé du Foresto,
W. de Wit,
C. Deen,
F. Delplancke-Ströbele,
J. Dexter
, et al. (48 additional authors not shown)
Abstract:
To investigate the inner regions of protoplanetary disks, we performed near-infrared interferometric observations of the classical TTauri binary system S CrA. We present the first VLTI-GRAVITY high spectral resolution ($R\sim$4000) observations of a classical TTauri binary, S CrA (composed of S CrA N and S CrA S and separated by $\sim$1.4"), combining the four 8-m telescopes in dual-field mode. Ou…
▽ More
To investigate the inner regions of protoplanetary disks, we performed near-infrared interferometric observations of the classical TTauri binary system S CrA. We present the first VLTI-GRAVITY high spectral resolution ($R\sim$4000) observations of a classical TTauri binary, S CrA (composed of S CrA N and S CrA S and separated by $\sim$1.4"), combining the four 8-m telescopes in dual-field mode. Our observations in the near-infrared K-band continuum reveal a disk around each binary component, with similar half-flux radii of about 0.1 au at d$\sim$130 pc, inclinations ($i=$28$\pm$3$^o$\ and $i=$22$\pm$6$^o$), and position angles (PA=0$^o\pm$6$^o$ and PA=-2$^o\pm$12$^o$), suggesting that they formed from the fragmentation of a common disk. The S CrA N spectrum shows bright HeI and Br$γ$ line emission exhibiting inverse P-Cygni profiles, typically associated with infalling gas. The continuum-compensated Br$γ$ line visibilities of S CrA N show the presence of a compact Br$γ$ emitting region the radius of which is about $\sim$0.06 au, which is twice as big as the truncation radius. This component is mostly tracing a wind. Moreover, a slight radius change between the blue- and red-shifted Br$γ$ line components is marginally detected. The presence of an inverse P-Cygni profile in the HeI and Br$γ$ lines, along with the tentative detection of a slightly larger size of the blue-shifted Br$γ$ line component, hint at the simultaneous presence of a wind and magnetospheric accretion in S CrA N.
△ Less
Submitted 5 September, 2017;
originally announced September 2017.
-
Robust Gaussian Stochastic Process Emulation
Authors:
Mengyang Gu,
Xiaojing Wang,
James O. Berger
Abstract:
We consider estimation of the parameters of a Gaussian Stochastic Process (GaSP), in the context of emulation (approximation) of computer models for which the outcomes are real-valued scalars. The main focus is on estimation of the GaSP parameters through various generalized maximum likelihood methods, mostly involving finding posterior modes; this is because full Bayesian analysis in computer mod…
▽ More
We consider estimation of the parameters of a Gaussian Stochastic Process (GaSP), in the context of emulation (approximation) of computer models for which the outcomes are real-valued scalars. The main focus is on estimation of the GaSP parameters through various generalized maximum likelihood methods, mostly involving finding posterior modes; this is because full Bayesian analysis in computer model emulation is typically prohibitively expensive. The posterior modes that are studied arise from objective priors, such as the reference prior. These priors have been studied in the literature for the situation of an isotropic covariance function or under the assumption of separability in the design of inputs for model runs used in the GaSP construction. In this paper, we consider more general designs (e.g., a Latin Hypercube Design) with a class of commonly used anisotropic correlation functions, which can be written as a product of isotropic correlation functions, each having an unknown range parameter and a fixed roughness parameter. We discuss properties of the objective priors and marginal likelihoods for the parameters of the GaSP and establish the posterior propriety of the GaSP parameters, but our main focus is to demonstrate that certain parameterizations result in more robust estimation of the GaSP parameters than others, and that some parameterizations that are in common use should clearly be avoided. These results are applicable to many frequently used covariance functions, e.g., power exponential, Mat{é}rn, rational quadratic and spherical covariance. We also generalize the results to the GaSP model with a nugget parameter. Both theoretical and numerical evidence is presented concerning the performance of the studied procedures.
△ Less
Submitted 15 August, 2017;
originally announced August 2017.
-
Supersymmetric Resonant Dark Matter: a Thermal Model for the AMS-02 Positron Excess
Authors:
Yang Bai,
Joshua Berger,
Sida Lu
Abstract:
We construct a thermal dark matter model with annihilation mediated by a resonance to explain the positron excess observed by PAMELA, Fermi-LAT and AMS-02, while satisfying constraints from cosmic microwave background (CMB) measurements. The challenging requirement is that the resonance has twice the dark matter mass to one part in a million. We achieve this by introducing an $SU(3)_f$ dark flavor…
▽ More
We construct a thermal dark matter model with annihilation mediated by a resonance to explain the positron excess observed by PAMELA, Fermi-LAT and AMS-02, while satisfying constraints from cosmic microwave background (CMB) measurements. The challenging requirement is that the resonance has twice the dark matter mass to one part in a million. We achieve this by introducing an $SU(3)_f$ dark flavor symmetry that is spontaneously broken to $SU(2)_f \times U(1)_f$. The resonance is the heaviest state in the dark matter flavor multiplet and the required mass relation is protected by the vacuum structure and supersymmetry from radiative corrections. The pseudo-Nambu Goldstone Bosons (PNGB's) from the dark flavor symmetry breaking can be slightly lighter than one GeV and dominantly decay into two muons just from kinematics, with subsequent decay into positrons. The PNGB's are produced in resonant dark matter semi-annihilation, where two dark matter particles annihilate into an anti-dark matter particle and a PNGB. The dark matter mass in our model is constrained to be below around 1.9 TeV from fitting thermal relic abundance, AMS-02 data and CMB constraints. The superpartners of Standard Model (SM) particles can cascade decay into a light PNGB along with SM particles, yielding a correlated signal of this model at colliders. One of the interesting signatures is a resonance of a SM Higgs boson plus two collimated muons, which has superb discovery potential at LHC Run 2.
△ Less
Submitted 29 June, 2017;
originally announced June 2017.
-
Advanced reduced-order models for moisture diffusion in porous media
Authors:
Suelen Gasparin,
Julien Berger,
Denys Dutykh,
Nathan Mendes
Abstract:
It is of great concern to produce numerically efficient methods for moisture diffusion through porous media, capable of accurately calculate moisture distribution with a reduced computational effort. In this way, model reduction methods are promising approaches to bring a solution to this issue since they do not degrade the physical model and provide a significant reduction of computational cost.…
▽ More
It is of great concern to produce numerically efficient methods for moisture diffusion through porous media, capable of accurately calculate moisture distribution with a reduced computational effort. In this way, model reduction methods are promising approaches to bring a solution to this issue since they do not degrade the physical model and provide a significant reduction of computational cost. Therefore, this article explores in details the capabilities of two model-reduction techniques - the Spectral Reduced-Order Model (Spectral-ROM) and the Proper Generalised Decomposition (PGD) - to numerically solve moisture diffusive transfer through porous materials. Both approaches are applied to three different problems to provide clear examples of the construction and use of these reduced-order models. The methodology of both approaches is explained extensively so that the article can be used as a numerical benchmark by anyone interested in building a reduced-order model for diffusion problems in porous materials. Linear and non-linear unsteady behaviors of unidimensional moisture diffusion are investigated. The last case focuses on solving a parametric problem in which the solution depends on space, time and the diffusivity properties. Results have highlighted that both methods provide accurate solutions and enable to reduce significantly the order of the model around ten times lower than the large original model. It also allows an efficient computation of the physical phenomena with an error lower than 10^{-2} when compared to a reference solution.
△ Less
Submitted 2 July, 2018; v1 submitted 16 May, 2017;
originally announced May 2017.
-
On the estimation of moisture permeability and advection coefficients of a wood fibre material using the optimal experiment design approach
Authors:
Julien Berger,
Thomas Busser,
Denys Dutykh,
Nathan Mendes
Abstract:
This paper presents a practical application of the concept of Optimal Experiment Design (OED) for the determination of properties of porous materials with in-situ measurements and an identification method. First, an experimental set-up was presented and used for the measurement of relative humidity within a wood fibre material submitted to single and multiple steps of relative humidity variation.…
▽ More
This paper presents a practical application of the concept of Optimal Experiment Design (OED) for the determination of properties of porous materials with in-situ measurements and an identification method. First, an experimental set-up was presented and used for the measurement of relative humidity within a wood fibre material submitted to single and multiple steps of relative humidity variation. Then, the application of OED enabled to plan the experimental conditions in terms of sensor positioning and boundary conditions out of 20 possible designs. The OED search was performed using the Fisher information matrix and a priori knowledge of the parameters. It ensures to provide the best accuracy of the identification method and thus the estimated parameter. Optimal design results have been found for single steps from the relative humidity phi = 10 to 75%, with one sensor located at the position X between 4 and 6 cm, for the estimation of moisture permeability coefficients, while from phi = 75% to phi = 33%, with one sensor located at X° = 3 cm, for the estimation of the advection coefficient. The OED has also been applied for the identification of couples of parameters. A sample submitted to multiple relative humidity steps (phi = 10-75-33-75%) with a sensor placed at X° = 5 cm was found as the best option for determining both properties with the same experiment. These OED parameters have then been used for the determination of moisture permeability and advection coefficients. The estimated moisture permeability coefficients are twice higher than the a priori values obtained using standard methods. The advection parameter corresponds to the mass average velocity of the order of v = 0.01 mm/s within the material and may play an important role on the simulation of moisture front.
△ Less
Submitted 9 July, 2017; v1 submitted 8 April, 2017;
originally announced May 2017.