Publications

We investigate the interplay of electronic correlations and spin-orbit coupling (SOC) for a one-band and a two-band honeycomb lattice model. The main difference between the two models concerning SOC is that in the one-band case the SOC is a purely non-local term in the basis of the p

Physical systems obey strict symmetry principles. We expect that machine learning methods that intrinsically respect these symmetries should have higher prediction accuracy and better generalization in prediction of physical dynamics. In this work we implement a principled model based on invariant scalars, and release open-source code. We apply this Scalars method to a simple chaotic dynamical system, the springy double pendulum. We show that the Scalars method outperforms state-of-the-art approaches for learning the properties of physical systems with symmetries, both in terms of accuracy and speed. Because the method incorporates the fundamental symmetries, we expect it to generalize to different settings, such as changes in the force laws in the system.

The BISOU (Balloon Interferometer for Spectral Observations of the Universe) project aims to study the viability and prospects of a balloon-borne spectrometer, pathfinder of a future space mission dedicated to the measurements of the CMB spectral distortions. We present here a preliminary concept based on previous space mission proposals, together with some sensitivity calculation results for the observation goals, showing that a 5-sigma measurement of the y-distortions is achievable.

Gas giant planets are expected to accrete most of their mass via a circumplanetary disk. If the planet is unmagnetized and initially slowly rotating, it will accrete gas via a radially narrow boundary layer and rapidly spin up. Radial broadening of the boundary layer as the planet spins up reduces the specific angular momentum of accreted gas, allowing the planet to find a terminal rotation rate short of the breakup rate. Here, we use axisymmetric viscous hydrodynamic simulations to quantify the terminal rotation rate of planets accreting from their circumplanetary disks. For an isothermal planet-disk system with a disk scale height h/r=0.1 near the planetary surface, spin up switches to spin down at between 70\% and 80\% of the planet's breakup angular velocity. In a qualitative difference from vertically-averaged models -- where spin down can co-exist with mass accretion -- we observe \emph{decretion} accompanying solutions where angular momentum is being lost. The critical spin rate depends upon the disk thickness near the planet. For an isothermal system with a disk scale height of h/r=0.15 near the planet, the critical spin rate drops to between 60\% and 70\% of the planet's breakup angular velocity. In the disk outside the boundary layer, we identify meridional circulation flows, which are unsteady and instantaneously asymmetric across the mid-plane. The simulated flows are strong enough to vertically redistribute solid material in early-stage satellite formation. We discuss how extrasolar planetary rotation measurements, when combined with spectroscopic and variability studies of protoplanets with circumplanetary disks, could determine the role of magnetic and non-magnetic processes in setting giant planet spins.

Gas giant planets are expected to accrete most of their mass via a circumplanetary disk. If the planet is unmagnetized and initially slowly rotating, it will accrete gas via a radially narrow boundary layer and rapidly spin up. Radial broadening of the boundary layer as the planet spins up reduces the specific angular momentum of accreted gas, allowing the planet to find a terminal rotation rate short of the breakup rate. Here, we use axisymmetric viscous hydrodynamic simulations to quantify the terminal rotation rate of planets accreting from their circumplanetary disks. For an isothermal planet-disk system with a disk scale height h/r = 0.1 near the planetary surface, spin-up switches to spin-down at between 70% and 80% of the planet's breakup rate. In a qualitative difference from vertically averaged models—where spin-down can coexist with mass accretion—we observe decretion accompanying solutions where angular momentum is being lost. The critical spin rate depends upon the disk thickness near the planet. For a disk scale height of h/r = 0.15, the critical spin rate drops to between 60% and 70% of the planet's breakup rate. In the disk outside the boundary layer, we identify meridional circulation flows, which are unsteady and instantaneously asymmetric across the midplane. The simulated flows are strong enough to vertically redistribute solid material in early stage satellite formation. We discuss how exoplanetary rotation measurements, when combined with spectroscopic and variability studies of protoplanets with circumplanetary disks, could determine the role of magnetic and nonmagnetic processes in setting planet spins.

The largest fluctuation in the CMB sky is the CMB dipole, which is believed to be caused by the motion of our observation frame with respect to the CMB rest frame. This motion accounts for the known motion of the Solar System barycentre with a best-fit amplitude of 369 km/s, in the direction (ℓ=264∘, b=48∘) in galactic coordinates. Along with the CMB dipole signal, this motion also causes an inevitable signature of statistical anisotropy in the higher multipoles due to the modulation and aberration of the CMB temperature and polarization fields. This leads to a correlation between adjacent CMB multipoles causing a non-zero value of the off-diagonal terms in the covariance matrix which can be captured in terms of the dipolar spectra of the bipolar spherical harmonics (BipoSH). In our work, we jointly infer the CMB power spectrum and the BipoSH spectrum in a Bayesian framework using the Planck-2018 𝚂𝙼𝙸𝙲𝙰 temperature map. We detect amplitude and direction of the local motion consistent with the canonical value v=369 km/s inferred from CMB dipole with a statistical significance of 4.54σ, 4.97σ and 5.23σ respectively from the masked temperature map with the available sky fraction 40.1%, 59.1%, and 72.2%, confirming the common origin of both the signals. The Bayes factor in favor of the canonical value is between 7 to 8 depending on the choice of mask. But it strongly disagrees (by a value of the Bayes factor about 10−10−10−11) with a higher value of local motion which one can infer from the amplitude of the dipole signal obtained from the CatWISE2020 quasar catalog using the WISE and NEOWISE data set.

The covariance matrix Σ of non-linear clustering statistics that are measured in current and upcoming surveys is of fundamental interest for comparing cosmological theory and data and a crucial ingredient for the likelihood approximations underlying widely used parameter inference and forecasting methods. The extreme number of simulations needed to estimate Σ to sufficient accuracy poses a severe challenge. Approximating Σ using inexpensive but biased surrogates introduces model error with respect to full simulations, especially in the non-linear regime of structure growth. To address this problem we develop a matrix generalization of Convergence Acceleration by Regression and Pooling (CARPool) to combine a small number of simulations with fast surrogates and obtain low-noise estimates of Σ that are unbiased by construction. Our numerical examples use CARPool to combine GADGET-III N-body simulations with fast surrogates computed using COmoving Lagrangian Acceleration (COLA). Even at the challenging redshift z=0.5, we find variance reductions of at least (101) and up to (104) for the elements of the matter power spectrum covariance matrix on scales 8.9×10−3<kmax<1.0 hMpc−1. We demonstrate comparable performance for the covariance of the matter bispectrum, the matter correlation function and probability density function of the matter density field. We compare eigenvalues, likelihoods, and Fisher matrices computed using the CARPool covariance estimate with the standard sample covariance estimators and generally find considerable improvement except in cases where Σ is severely ill-conditioned.

The angular size of the broad line region (BLR) of the nearby active galactic nucleus (AGN) NGC 3783 has been spatially resolved by recent observations with VLTI/GRAVITY. A reverberation mapping (RM) campaign has also recently obtained high quality light curves and measured the linear size of the BLR in a way that is complementary to the GRAVITY measurement. The size and kinematics of the BLR can be better constrained by a joint analysis that combines both GRAVITY and RM data. This, in turn, allows us to obtain the mass of the supermassive black hole in NGC3783 with an accuracy that is about a factor of two better than that inferred from GRAVITY data alone. We derive MBH=2.54+0.90−0.72×107M⊙. Finally, and perhaps most notably, we are able to measure a geometric distance to NGC 3783 of 39.9+14.5−11.9 Mpc. We are able to test the robustness of the BLR-based geometric distance with measurements based on the Tully-Fisher relation and other indirect methods. We find the geometric distance is consistent with other methods within their scatter. We explore the potential of BLR-based geometric distances to directly constrain the Hubble constant, H0, and identify differential phase uncertainties as the current dominant limitation to the H0 measurement precision for individual sources.

Disk accretion onto weakly magnetized objects possessing a material surface must proceed via the so-called boundary layer (BL) - a region at the inner edge of the disk, in which the velocity of accreting material abruptly decreases from its Keplerian value. Supersonic shear arising in the BL is known to be conducive to excitation of acoustic waves that propagate into both the accretor and the disk, enabling angular momentum and mass transport across the BL. We carry out a numerical exploration of different wave modes that operate near the BL, focusing on their morphological characteristics in the innermost parts of accretion disk. Using a large suite of simulations covering a broad range of Mach numbers (of the supersonic shear flow in the BL), we provide accurate characterization of the different types of modes, verifying their properties against analytical results, when available. We discover new types of modes, in particular, global spiral density waves launched by vortices forming in the disk near the BL as a result of the Rossby wave instability; this instability is triggered by the vortensity production in that region caused by the nonlinear damping of acoustic waves. Azimuthal wavenumbers of the dominant modes that we observe appear to increase monotonically with the Mach number of the runs, but a particular mix of modes found in a simulation is mildly stochastic. Our results provide a basis for better understanding of the angular momentum and mass transport across the BL as well as the emission variability in accreting objects.

Many sensory-driven behaviors rely on predictions about future states of the environment. Visual input typically evolves along complex temporal trajectories that are difficult to extrapolate. We test the hypothesis that spatial processing mechanisms in the early visual system facilitate prediction by constructing neural representations that follow straighter temporal trajectories. We recorded V1 population activity in anesthetized macaques while presenting static frames taken from brief video clips, and developed a procedure to measure the curvature of the associated neural population trajectory. We found that V1 populations straighten naturally occurring image sequences, but entangle artificial sequences that contain unnatural temporal transformations. We show that these effects arise in part from computational mechanisms that underlie the stimulus selectivity of V1 cells. Together, our findings reveal that the early visual system uses a set of specialized computations to build representations that can support prediction in the natural environment.