Showing 201–250 of 485 results for author: Berger, J

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.

Comments: Resubmission to SciPost 8 pages, 3 figures

Journal ref: SciPost Phys. Core 1, 001 (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.

Comments: version 2: 15 pages, 1 figure, matches PRD version

Journal ref: Phys. Rev. D 100, 055005 (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.

Comments: 35 pages, 19 figures, 2 tables, 42 references. Other author's papers can be downloaded at http://www.denys-dutykh.com/

Journal ref: Building Simulation (2020), Vol. 13, pp. 1-18

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.

Comments: 29 pages, 20 figures, 4 tables, 1 appendix, 29 references. Other author's papers can be downloaded at http://www.denys-dutykh.com/

Journal ref: Computational Technologies (2019), Vol. 24, Number 3, pp. 4-32

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.

Comments: 9 pages, 7 figures, submitted to A&A

Journal ref: 2019, A&A, 625, L10

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.

Comments: Submitted to DAFx 2019, the 22nd International Conference on Digital Audio Effects, Birmingham, United Kingdom

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.

Comments: published in A&A

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.

Comments: 46 pages, 12 figures, 3 tables, 46 references. Other author's papers can be downloaded at http://www.denys-dutykh.com/

Journal ref: International Journal of Heat and Mass Transfer (2019), Vol. 134, pp.1041-1060

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.

Comments: 36 pages, 18 figures, 6 tables, 42 references. Other author's papers can be downloaded at http://www.denys-dutykh.com/

Journal ref: Journal of the Brazilian Society of Mechanical Sciences and Engineering (2019), Vol. 41, Article number: 173

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.

Comments: 41 pages, 19 figures, 3 tables, 35 references. Other author's papers can be downloaded at http://www.denys-dutykh.com/

Journal ref: Numerical Heat Transfer, Part A: Applications (2019), Vol. 75, Issue 2, pp. 110-136

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.

Comments: 45 pages, 21 figures, 5 tables, 47 references. Other author's papers can be downloaded at http://www.denys-dutykh.com/

Journal ref: International Journal of Thermal Sciences (2018), Vol. 133, pp.120-139

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.

Comments: Accepted for publication in Physical Review Letters

Journal ref: Phys. Rev. Lett. 122, 101102 (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.

Comments: 11 pages, 9 figures, accepted in MNRAS

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.

Comments: Accepted by MNRAS

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.

Comments: 21 pages; generalized Lemma 3.7 to the non-unimodular setting

Journal ref: J. Algebra 548 (2020) 96-119

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.

Comments: 16 pages, 2 figures

Report number: PITT-PACC-1822

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.

Comments: 34 pages

Report number: PITT-PACC-1818

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.

Journal ref: Zeitschrift für Naturforschung B, Volume 75: Issue 4, 2020, 327--393

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.

Comments: 5 pages, 1 figure

Journal ref: Phys. Rev. B 99, 205144 (2019)

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.

Comments: 19 pages, 6 figures

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.

Comments: accepted by A&A; 16 pages

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.

Comments: Accepted for publication in A&A

Journal ref: A&A 620, A116 (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.

Comments: 32 pages, 27 figures (14 in the main body, 13 in the appendices), accepted for publication in AA

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.

Comments: Published in Experimental Astronomy as part of topical collection "Future of Optical-infrared Interferometry in Europe"

Journal ref: Monnier, J.D., Kraus, S., Ireland, M.J. et al. Exp Astron (2018). https://doi.org/10.1007/s10686-018-9594-1

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.

Comments: Accepted for publication in A&A Letters, 29 June 2018, 10 pages, 6 figures, corresponding author: F. Eisenhauer <[email protected]>

Journal ref: A&A 615, L15 (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.

Comments: 7 pages, 5 figures

Journal ref: J. Chem. Theory Comput. 14, 5220 (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.

Comments: 15 pages, 3 figures, SPIE proceedings, based on arXiv:1801.04148

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.

Comments: Accepted for publication in A&A. 26 pages, 22 figures, 12 tables

Journal ref: A&A 618, A23 (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.

Comments: 51 pages, 21 figures, 3 tables, 44 references. Other author's publications can be downloaded at http://www.denys-dutykh.com/

Journal ref: Transport in Porous Media (2018) Vol. 121, pp. 665-702

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.

Comments: 5 + 19 pages

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.

Comments: 33 pages, 10 figures, 9 tables, 29 references. Other author's papers can be downloaded at http://www.denys-dutykh.com/

Journal ref: Building and Environment (2018), Vol. 138, pp. 262-274

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.

Comments: 11 pages, 7 figures

Journal ref: J. Chem. Theory Comput. 14, 3071 (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.

Comments: 24 pages, 2 figures, accepted for publication in Experimental Astronomy

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.

Journal ref: Phys. Rev. B 98, 024402 (2018)

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.

Comments: 6 PAGES, 4 FIGURES

Journal ref: Phys. Rev. B 97, 155102 (2018)

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.

Journal ref: J. Phys.: Condens. Matter 30 135602 (2018)

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.

Comments: 19 pages, 8 figures, accepted for publication in ApJ

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.

Comments: Accepted by A&A

Journal ref: A&A 608, A78 (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.

Comments: 34 pages, 11 figures

Journal ref: Phys. Rev. D 97, 115012 (2018)

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.

Comments: 42 pages, 14 figures, 1 table, 69 references. Other author's papers can be downloaded at http://www.denys-dutykh.com/. arXiv admin note: text overlap with arXiv:1704.07607

Journal ref: Transport in Porous Media (2018), Vol. 124, pp. 965-994

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.

Comments: 34 pages, 14 figures, 6 tables, 43 references. Other author's papers can be downloaded at http://www.denys-dutykh.com/. arXiv admin note: text overlap with arXiv:1704.06565

MSC Class: 35R30 (primary); 35K05; 80A20; 65M32 (secondary)

Journal ref: Experimental Thermal and Fluid Science (2018), Vol. 90, pp. 246-259