Tuesdays 10:30 - 11:30 | Fridays 11:30 - 12:30
Showing votes from 2023-12-01 12:30 to 2023-12-05 11:30 | Next meeting is Tuesday Oct 29th, 10:30 am.
The intergalactic helium became fully ionized by the end of cosmic noon ($z\sim2$). Similarly to the reionization of hydrogen, helium reionization is expected to be patchy, driven by luminous quasars that ionize the intergalactic gas in their surrounding environment. Probing the morphology of ionized electrons during this epoch can provide crucial information about early structure formation, including the clustering and luminosities of quasars, the accretion rates, variability, and lifetimes of active galactic nuclei, as well as the growth and evolution of supermassive black holes. In this study, we present how measurements of the cosmic microwave background (CMB) can be used to reconstruct the optical-depth fluctuations resulting from patchy helium reionization. As helium reionization occurred at lower redshifts, upcoming probes of large-scale structure surveys will present a significant opportunity to enhance the prospects of probing this epoch by their combined analysis with the CMB. Using a joint information-matrix analysis of hydrogen and helium reionization, we show that near-future galaxy and CMB surveys will have enough statistical power to detect optical-depth fluctuations due to doubly-ionized helium, providing a way of measuring the redshift and duration of helium reionization to high significance. We also show that modeling uncertainties in helium reionization can impact the measurement precision of parameters characterizing hydrogen reionization.
Temporary matter domination and late entropy dilution, injected by a "long-lived" particle in the early universe, serves as a standard mechanism for yielding the correct dark matter relic density. We recently pointed out the cosmological significance of diluting particle's partial decay into dark matter. When repopulated in such a way, dark matter carries higher momentum than its thermal counterpart, resulting in a suppression of the linear matter power spectrum that is constrained by the large scale structure observations. In this work, we study the impact of such constraints on the minimal left-right symmetric model that accounts for the origin of neutrino mass. We map out a systematic anatomy of possible dilution scenarios with viable parameter spaces, allowed by cosmology and various astrophysical and terrestrial constraints. We show that to accommodate the observed dark matter relic abundance the spontaneous left-right symmetry breaking scale must be above PeV and cosmology will continue to provide the most sensitive probes of it. In case the dilutor is one of the heavier right-handed neutrinos, it can be much lighter and lie near the electroweak scale.
The Scalar Field Dark Matter model has been known in various ways throughout its history; Fuzzy, BEC, Wave, Ultralight, Axion-like Dark Matter, etc. All of them consist in proposing that the dark matter of the universe is a spinless field $\Phi$ that follows the Klein-Gordon (KG) equation of motion $\Box\Phi-dV/d\Phi=0$, for a given scalar field potential $V$. The difference between different models is sometimes the choice of the scalar field potential $V$. In the literature we find that people usually work in the nonrelativistic, weak-field limit of the KG equation where it transforms into the Schr\"odinger equation and the Einstein equations into the Poisson equation, reducing the KG-Einstein system, to the Schr\"odinger-Poisson system. In this paper, we review some of the most interesting achievements of this model from the historical point of view and its comparison with observations, showing that this model could be the last answer to the question about the nature of dark matter in the universe.
This paper presents analytic results on the anatomy of nested sampling, from which a technique is developed to estimate the run-time of the algorithm that works for any nested sampling implementation. We test these methods on both toy models and true cosmological nested sampling runs. The method gives an order-of-magnitude prediction of the end point at all times, forecasting the true endpoint within standard error around the halfway point.
Creating radio galaxy catalogues from next-generation deep surveys requires automated identification of associated components of extended sources and their corresponding infrared hosts. In this paper, we introduce RadioGalaxyNET, a multimodal dataset, and a suite of novel computer vision algorithms designed to automate the detection and localization of multi-component extended radio galaxies and their corresponding infrared hosts. The dataset comprises 4,155 instances of galaxies in 2,800 images with both radio and infrared channels. Each instance provides information about the extended radio galaxy class, its corresponding bounding box encompassing all components, the pixel-level segmentation mask, and the keypoint position of its corresponding infrared host galaxy. RadioGalaxyNET is the first dataset to include images from the highly sensitive Australian Square Kilometre Array Pathfinder (ASKAP) radio telescope, corresponding infrared images, and instance-level annotations for galaxy detection. We benchmark several object detection algorithms on the dataset and propose a novel multimodal approach to simultaneously detect radio galaxies and the positions of infrared hosts.
Continuing work presented in Li et al. (2021), we performed a series of tests to our high-resolution three-dimensional mass map reconstruction algorithm \splinv{}. We test the mass reconstruction accuracy against realistic mock catalogs generated using shear field produced by triaxial halos with the inner density profile of $\rho \propto r^{-1}$ and of $\rho \propto r^{-1.5}$. The galaxy shape noise is modeled based on the Year-1 Subaru Hyper Suprime-Cam (HSC) Survey. After reviewing mathematical details of our algorithm and dark matter halo models, we determine an optimal value of the coefficient of the adaptive LASSO regression penalty term for single halo reconstruction. We successfully measure halo masses for massive triaxial halos; the mass determination accuracy is 5 percent for halos with $M = 10^{14.6}~M_\odot$ at $0.0625\leq z \leq 0.2425$, and 5 percent for those with $10^{14.8}~M_\odot$ at $0.0625\leq z \leq 0.4675$, and 20 percent for $M= 10^{15.0} ~M_\odot$ and $M=10^{15.2}~M_\odot$ in the redshift range $0.0625\leq z \leq 0.4675$. The redshift estimate accuracy is consistently below $\Delta z /z \leq 0.05$ for the above halo masses in the range $0.1525\leq z \leq 0.4675$. We also demonstrate that the orientation of triaxial halos and systematic error in our halo model do not affect reconstruction result significantly. Finally, we present results from reconstruction of mass distribution using shear catalogs produced by multiple halos, to show \splinv{}'s capability using realistic shear maps from ongoing and future galaxy surveys.
The probability that the primordial black hole (PBH) binaries formed in the early Universe can be affected by the Hubble expansion of background, which is non-negligible when the number density of PBHs is very low (it is actually this case for supermassive PBHs). In this paper, taking into account the effect of cosmic expansion on the comoving distance of PBH pairs, we worked out the merger rate of PBHs with any extended mass function. The torques by all PBHs and linear density perturbations are also considered. It is found that the merger rate of PBH, $M\gtrsim 10^6M_\odot$, binaries is significantly lower for $f_\text{pbh}\lesssim 0.01$ than expected.
In contrast to our understanding of density field tracers, the modelling of direct statistics pertaining to the cosmic velocity field remains open to significant opportunities for improvement. The lack of accurate modelling for the non-linear domain of pairwise velocities restricts our capacity to fully exploit the information encoded in this observable. We present a robust approach for modelling the mean infall velocities, $v_{12}(r,a)$, with broad applicability spanning sub-megaparsec scales and cosmologies extending beyond the standard LCDM paradigm. Our approach involves solving the full pair-conservation equation using accurate non-linear power spectrum descriptions. To assess the robustness of our model, we extend it to cosmologies beyond the standard LCDM, in particular, the Hu-Sawicki $f(R)$-gravity and Dvali-Gabadadze-Porrati (DGP) modified gravity models. Remarkably, our predictions for pairwise velocities of dark matter particles at kilo-parsec scales exhibit excellent agreement with N-body simulations throughout the entire dynamical range ($0.1 \lesssim \xi \lesssim 1000$, or $r\geq0.4$Mpc/h). Furthermore, we show that different gravity models leave distinct signatures in the shape and dynamics of the mean pairwise velocities, providing a potent test of cosmological gravity laws.
The background of primordial gravitational waves (PGWs) predicted by the inflationary scenario induce incoherence of the electromagnetic field propagating over cosmological distances. We propose a new schema to constrain the underlying inflationary parameters, in particular the tensor-to-scalar ratio r, based on angular size-redshift \theta-z measurement accomplished by very long baseline interferometry (VLBI) surveys. VLBI observations rely on the van-Citter Zernike theorem, which expresses the coherence length in terms of its wavelength and angular size, i.e., d_coh < \lambda/\theta. In this study, we show that the interaction of the radio signal (involved in VLBI) with PGWs, along the propagation from a source located at redshift z_* to the Earth, leads to the blurring of the visibility. The blurring effect is evaluated for the highly-squeezed PGW, where it turns out that in order not to ruin the visibility, the projected baseline of the interferometer must be smaller than \xi_sq(z), which is inversely proportional to the the tensor-to-scalar ratio through \xi_sq \propto r^{-1/2}. Hence, VLBI observations based on interference pattern lead to a constraint on r imposed by the fact that \xi_sq(z) is greater than d_coh of the emitted radiation. In order to evaluate the constraint, we use a sample of compact radio quasars observed in VLBI and located at redshift range 0.46<z<2.73. We obtain a stringent upper-bound on the tensor-to-scalar ratio, r<2.10^{-6}, far beyond present constraints on this parameter. Further issues and caveats that potentially affect the results are reviewed. In particular, the possible effect of quantum-to-classical transition of PGWs is discussed. Ultimately, the background of primordial tensor perturbations may be more constrained with the help of the high-precision VLBI measurement of angular size-redshift of more distant sources.
We investigate the stochastic gravitational waves background arising from the first-order QCD chiral phase transition, considering three distinct sources: bubble collisions, sound waves, and fluid turbulence. Within the framework of the Polyakov-Nambu-Jona-Lasinio (PNJL) model, we calculate the parameters governing the intensity of the gravitational wave phase transition and determine their magnitudes along the adiabatic evolutionary path. We introduce the effective bag constant $B_{\mathrm{eff}}$ related to the dynamical evolution of quarks to evaluate the intensity of the phase transition. By calculating expanded potential at the point of false vacuum, we find that all the bubbles are in the mode of runaway, leading the velocity of the bubble wall to the speed of light. The resulting gravitational wave energy spectrum is estimated, revealing a characteristic amplitude of the generated gravitational waves within the centihertz frequency range. We present the gravitational wave spectrum and compare it with the sensitivity range of detectors, and find that the gravitational wave spectra generated by these sources have the potential to be detected by future detectors such as BBO and $\mu$ARES.
In this contribution to the proceedings of the 182nd Nobel Symposium, I reflect on the concept of "discovery" as it is used by physicists and astronomers. In particular, I comment on how the scientific community distinguishes discoveries from propositions that are supported only by lesser forms of evidence, emphasizing the social nature of this process and remarking on the subjective factors that go into making such judgements. I advocate for an approach that is intentionally Bayesian in nature, in which individuals are encouraged to evaluate and publicly state their priors and to update them systematically. I close by applying these practices to the case example of the Galactic Center Gamma-Ray Excess.
We investigate the accuracy of the perturbative galaxy bias expansion in view of the forthcoming analysis of the Euclid spectroscopic galaxy samples. We compare the performance of an Eulerian galaxy bias expansion, using state-of-art prescriptions from the effective field theory of large-scale structure (EFTofLSS), against a hybrid approach based on Lagrangian perturbation theory and high-resolution simulations. These models are benchmarked against comoving snapshots of the Flagship I N-body simulation at $z=(0.9,1.2,1.5,1.8)$, which have been populated with H$\alpha$ galaxies leading to catalogues of millions of objects within a volume of about $58\,h^{-3}\,{\rm Gpc}^3$. Our analysis suggests that both models can be used to provide a robust inference of the parameters $(h, \omega_{\rm c})$ in the redshift range under consideration, with comparable constraining power. We additionally determine the range of validity of the EFTofLSS model in terms of scale cuts and model degrees of freedom. From these tests, it emerges that the standard third-order Eulerian bias expansion can accurately describe the full shape of the real-space galaxy power spectrum up to the maximum wavenumber $k_{\rm max}=0.45\,h\,{\rm Mpc}^{-1}$, even with a measurement precision well below the percent level. In particular, this is true for a configuration with six free nuisance parameters, including local and non-local bias parameters, a matter counterterm, and a correction to the shot-noise contribution. Fixing either tidal bias parameters to physically-motivated relations still leads to unbiased cosmological constraints. We finally repeat our analysis assuming a volume that matches the expected footprint of Euclid, but without considering observational effects, as purity and completeness, showing that we can get consistent cosmological constraints over this range of scales and redshifts.
We study the possibility of using the $LiteBIRD$ satellite $B$-mode survey to constrain models of inflation producing specific features in CMB angular power spectra. We explore a particular model example, i.e. spectator axion-SU(2) gauge field inflation. This model can source parity-violating gravitational waves from the amplification of gauge field fluctuations driven by a pseudoscalar "axionlike" field, rolling for a few e-folds during inflation. The sourced gravitational waves can exceed the vacuum contribution at reionization bump scales by about an order of magnitude and can be comparable to the vacuum contribution at recombination bump scales. We argue that a satellite mission with full sky coverage and access to the reionization bump scales is necessary to understand the origin of the primordial gravitational wave signal and distinguish among two production mechanisms: quantum vacuum fluctuations of spacetime and matter sources during inflation. We present the expected constraints on model parameters from $LiteBIRD$ satellite simulations, which complement and expand previous studies in the literature. We find that $LiteBIRD$ will be able to exclude with high significance standard single-field slow-roll models, such as the Starobinsky model, if the true model is the axion-SU(2) model with a feature at CMB scales. We further investigate the possibility of using the parity-violating signature of the model, such as the $TB$ and $EB$ angular power spectra, to disentangle it from the standard single-field slow-roll scenario. We find that most of the discriminating power of $LiteBIRD$ will reside in $BB$ angular power spectra rather than in $TB$ and $EB$ correlations.
While the simplest inflationary models predict the primordial perturbations to be near scale-invariant, the primordial power spectrum (PPS) can exhibit oscillatory features in many physically well-motivated models. We search for hints of such features via free-form reconstructions of the PPS based on \textit{Planck} 2018 CMB temperature and polarization anisotropies. In order to robustly invert the oscillatory integrals and handle noisy unbinned data, we draw inspiration from image analysis techniques. In previous works, the Richardson-Lucy deconvolution algorithm for deblurring images has been modified for reconstructing PPS from the CMB temperature angular power spectrum. We extensively develop the methodology by including CMB polarization and introducing two new regularization techniques, also inspired by image analysis and adapted for our cosmological context. Regularization is essential for improving the fit to the temperature and polarization channels (TT, TE and EE) simultaneously without sacrificing one for another. The reconstructions we obtain are consistent with previous findings from temperature-only analyses. We evaluate the statistical significance of the oscillatory features in our reconstructions using mock data and find the observations to be consistent with having a featureless PPS. The machinery developed here will be a complimentary tool in the search for features with upcoming CMB surveys. Our methodology also shows competitive performance in image deconvolution tasks, which have various applications from microscopy to medical imaging.
Cosmological analyses of second-order weak lensing statistics require precise and accurate covariance estimates. These covariances are impacted by two sometimes neglected terms: A negative contribution to the Gaussian covariance due to finite survey area and the super-sample covariance (SSC) which for the power spectrum contains the impact by Fourier modes larger than the survey window. We show here that these two effects are connected and can be seen as correction terms to the "large-field-approximation", the asymptotic case of an infinitely large survey area. We describe the two terms collectively as "Finite-Field-Terms". We derive the covariance of second-order shear statistics from first principles. For this, we use an estimator in real space without relying on an estimator for the power spectrum. The resulting covariance does not scale inversely with the survey area, as naively assumed. This scaling is only correct under the large-field approximation when the contribution of the finite-field terms tends to zero. Furthermore, all parts of the covariance, not only the SSC, depend on the power- and trispectrum at all modes, including those larger than the survey. We also show that it is generally impossible to transform an estimate for the power spectrum covariance into the covariance of a real-space statistic. Such a transformation is only possible in the asymptotic case of the "large-field approximation". Additionally, we find that the total covariance of a real-space statistic can be calculated using correlation functions estimates on spatial scales smaller than the survey window. Consequently, estimating covariances of real-space statistics, in principle, does not require information on spatial scales larger than the survey area. We demonstrate that this covariance estimation method is equivalent to the standard sample covariance method.
Lorentz invariance is one of the fundamental tenets of Special Relativity, and has been extensively tested with laboratory and astrophysical observations. However, many quantum gravity models and theories beyond the Standard Model of Particle Physics predict a violation of Lorentz invariance at energies close to the Planck scale. This article reviews observational and experimental tests of Lorentz invariance violation (LIV) with photons, neutrinos and gravitational waves. Most astrophysical tests of LIV using photons are based on searching for a correlation of the spectral lag data with redshift and energy. These have been primarily carried out using compact objects such as pulsars, Active Galactic Nuclei (AGN), and Gamma-ray bursts (GRB). There have also been some claims for LIV from some of these spectral lag observations with GRBs, which however are in conflict with the most stringent limits obtained from other LIV searches. Searches have also been carried out using polarization measurements from GRBs and AGNs. For neutrinos, tests have been made using both astrophysical observations at MeV energies (from SN 1987A) as well as in the TeV-PeV energy range based on IceCube observations, atmospheric neutrinos, and long-baseline neutrino oscillation experiments. Cosmological tests of LIV entail looking for a constancy of the speed of light as a function of redshift using multiple observational probes, as well as looking for birefringence in Cosmic Microwave Background observations. This article will review all of these aforementioned observational tests of LIV, including results which are in conflict with each other.
We perform a blinded cosmology analysis with cosmic shear two-point correlation functions (2PCFs) measured from more than 25 million galaxies in the Hyper Suprime-Cam three-year shear catalog in four tomographic redshift bins ranging from 0.3 to 1.5. After conservative masking and galaxy selection, the survey covers 416 deg$^2$ of the northern sky with an effective galaxy number density of 15 arcmin$^{-2}$ over the four redshift bins. The 2PCFs adopted for cosmology analysis are measured in the angular range: $7.1 < \theta/{\rm arcmin} < 56.6$ for $\xi_+$ and $31.2 <\theta/{\rm arcmin} < 248$ for $\xi_-$, with a total signal-to-noise ratio of 26.6. We apply a conservative, wide, flat prior on the photometric redshift errors on the last two tomographic bins, and the relative magnitudes of the cosmic shear amplitude across four redshift bins allow us to calibrate the photometric redshift errors. With this flat prior on redshift errors, we find $\Omega_{\rm m}=0.256_{-0.044}^{+0.056}$ and $S_8\equiv \sigma_8 \sqrt{\Omega_{\rm m}/0.3}=0.769_{-0.034}^{+0.031}$ (both 68\% CI) for a flat $\Lambda$ cold dark matter cosmology. We find, after unblinding, that our constraint on $S_8$ is consistent with the Fourier space cosmic shear and the 3$\times$2pt analyses on the same HSC dataset. We carefully study the potential systematics from astrophysical and systematic model uncertainties in our fiducial analysis using synthetic data, and report no biases (including projection bias in the posterior space) greater than $0.5\sigma$ in the estimation of $S_8$. Our analysis hints that the mean redshifts of the two highest tomographic bins are higher than initially estimated. In addition, a number of consistency tests are conducted to assess the robustness of our analysis. Comparing our result with Planck-2018 cosmic microwave background observations, we find a ~$2\sigma$ tension for the $\Lambda$CDM model.
The rapid development of gravitational wave astronomy provides the unique opportunity of exploring the dynamics of the Universe using clustering properties of coalescing binary black hole mergers. Gravitational wave data, along with information coming from future galaxy surveys, have the potential of shedding light about many open questions in Cosmology, including those regarding the nature of dark matter and dark energy. In this work we explore which combination of gravitational wave and galaxy survey datasets are able to provide the best constraints both on modified gravity theories and on the nature of the very same binary black hole events. In particular, by using the public Boltzmann code \texttt{Multi\_CLASS}, we compare cosmological constraints on popular $\Lambda$CDM extensions coming from gravitational waves alone and in conjunction with either deep and localized or wide and shallow galaxy surveys. We show that constraints on extensions of General Relativity will be at the same level of existing limits from gravitational waves alone or one order of magnitude better when galaxy surveys are included. Furthermore, cross-correlating both kind of galaxy survey with gravitational waves datasets will allow to confidently rule in or out primordial black holes as dark matter candidate in the majority of the allowed parameter space.
We show that -- in the framework of general relativity (GR) -- if black holes (BHs) are singularity-free objects, they couple to the large-scale cosmological dynamics. We find that the leading contribution to the resulting growth of the BH mass ($M_{\rm BH}$) as a function of the scale factor $a$ stems from the curvature term, yielding $M_{\rm BH} \propto a^k$, with $k=1$. We demonstrate that such a linear scaling is universal for spherically-symmetric objects, and it is the only contribution in the case of regular BHs. For nonsingular horizonless compact objects we instead obtain an additional subleading model-dependent term. We conclude that GR nonsingular BHs/horizonless compact objects, although cosmologically coupled, are unlikely to be the source of dark energy. We test our prediction with astrophysical data by analysing the redshift dependence of the mass growth of supermassive BHs in a sample of elliptical galaxies at redshift $z=0.8 -0.9$. We also compare our theoretical prediction with higher redshift BH mass measurements obtained with the James Webb Space Telescope (JWST). We find that, while $k=1$ is compatible within $1 \sigma$ with JWST results, the data from elliptical galaxies at $z=0.8 -0.9$ favour values of $k>1$. New samples of BHs covering larger mass and redshift ranges and more precise BH mass measurements are required to settle the issue.
Weak gravitational lensing induces flux dependent fluctuations in the observed galaxy number density distribution. This cosmic magnification (magnification bias) effect in principle enables lensing reconstruction alternative to cosmic shear and CMB lensing. However, the intrinsic galaxy clustering, which otherwise overwhelms the signal, has hindered its application. Through a scaling relation found by principal component analysis of the galaxy clustering in multi-band photometry space, we design a minimum variance linear estimator to suppress the intrinsic galaxy clustering and to reconstruct the lensing convergence map. In combination of the CosmoDC2 galaxy mock and the CosmicGrowth simulation, we test this proposal for a LSST-like galaxy survey with $ugrizY$ photometry bands. The scaling relation holds excellently at multipole $\ell<10^3$, and remains reasonably well to $\ell\sim 3000$. The linear estimator efficiently suppresses the galaxy intrinsic clustering, by a factor of $\sim 10^2$. For galaxies in the photo-z range $0.8<z_\kappa<1.2$, the reconstructed convergence map is cosmic variance limited per $\ell$ mode at $\ell<10^2$, and shot noise limited at $\ell>= 200$. Its cross-correlation with cosmic shear of galaxies can achieve $S/N >= 200$. When the source redshift of cosmic shear galaxies $z_\gamma<z_\kappa$, the systematic error is negligible at all investigated scales ($\ell<3000$). When $z_\gamma\geq z_\kappa$, the systematic error caused by the residual intrinsic galaxy clustering becomes non-negligible. We discuss possible mitigation of the residual intrinsic galaxy clustering required for accurate measurement at $\ell>10^3$. This work further demonstrates the potential of lensing measurement through cosmic magnification to enhance the weak lensing cosmology.
Together with larger spectroscopic surveys such as the Dark Energy Spectroscopic Instrument (DESI), the precision of large scale structure studies and thus the constraints on the cosmological parameters are rapidly improving. Therefore, one must build realistic simulations and robust covariance matrices. We build galaxy catalogues by applying a Halo Occupation Distribution (HOD) model upon the \textsc{FastPM} simulations, such that the resulting galaxy clustering reproduces high resolution $N$-body simulations. While the resolution and halo finder are different from the reference simulations, we reproduce the reference galaxy two-point clustering measurements -- monopole and quadrupole -- to a precision required by the DESI Year 1 Emission Line Galaxy sample down to non-linear scales, i.e. $k<0.5\,h\mathrm{Mpc}$ or $s>10\,\mathrm{Mpc}/h$. Furthermore, we compute covariance matrices based on the resulting \textsc{FastPM} galaxy clustering -- monopole and quadrupole. We study for the first time the effect of fitting on Fourier conjugate [e.g. power spectrum] on the covariance matrix of the Fourier counterpart [e.g. correlation function]. We estimate the uncertainties of the two parameters of a simple clustering model and observe a maximum variation of 20 per cent for the different covariance matrices. Nevertheless, for most studied scales the scatter is between two to ten per cent Consequently, using the current pipeline we can precisely reproduce the clustering of $N$-body simulations and the resulting covariance matrices provide robust uncertainty estimations against HOD fitting scenarios. We expect our methodology will be useful for the coming DESI data analyses and their extension for other studies.
Pulsar Timing Array collaborations have recently announced the discovery of a stochastic gravitational wave background (SGWB) at nanohertz frequencies. We analyze the GW signals from the domination of ultra-low mass primordial black holes (PBHs) in the early universe and show that they can explain this recent discovery. This scenario requires a relatively broad peak in the power spectrum of scalar perturbations from inflation with a spectral index in a narrow range of $1.45$ to $1.6$. The resulting PBH population would have mass around $10^{8}$g, and the initial abundance $\beta_f$ lies between $10^{-10}$ and $10^{-9}$. We find that this explanation is preferred by the data over the generic model, assuming supermassive BHs as the source. These very light PBHs would decay before Big Bang Nucleosynthesis (BBN); however, upcoming third-generation terrestrial laser interferometers would be able to test the model by observing the GW spectrum produced during the formation of the PBHs. Also, the scalar power spectra associated with our scenario will be within the reach of PIXIE probing CMB spectral distortions.
We present two galaxy samples, selected from DESI Legacy Imaging Surveys (LS) DR9, with approximately 20,000 square degrees of coverage and spectroscopic redshift distributions designed for cross-correlations such as with CMB lensing, galaxy lensing, and the Sunyaev-Zel'dovich effect. The first sample is identical to the DESI Luminous Red Galaxy (LRG) sample, and the second sample is an extended LRG sample with 2-3 times the DESI LRG density. We present the improved photometric redshifts, tomographic binning and their spectroscopic redshift distributions and imaging systematics weights, and magnification bias coefficients. The catalogs and related data products will be made publicly available. The cosmological constraints using this sample and Planck lensing maps are presented in a companion paper. We also make public the new set of general-purpose photometric redshifts trained using DESI spectroscopic redshifts, which are used in this work, for all galaxies in LS DR9.
The Square Kilometer array is expected to measure the 21cm signal from the Epoch of Reionization (EoR) in the coming decade, and its pathfinders may provide a statistical detection even earlier. The currently reported upper limits provide tentative constraints on the astrophysical parameters of the models of the EoR. In order to interpret such data with 3D radiative hydrodynamics simulations using Bayesian inference, we present the latest developments of the \textsc{Licorice} code. Relying on an implementation of the halo conditional mass function to account for unresolved star formation, this code now allows accurate simulations of the EoR at $256^3$ resolution. We use this version of \textsc{Licorice} to produce the first iteration of \textsc{LoReLi}, a public dataset now containing hundreds of 21cm signals computed from radiative hydrodynamics simulations. We train a neural network on \textsc{LoReLi} to provide a fast emulator of the \textsc{Licorice} power spectra, \textsc{LorEMU}, which has $\sim 5\%$ rms error relative to the simulated signals. \textsc{LorEMU} is used in a Markov Chain Monte Carlo framework to perform Bayesian inference, first on a mock observation composed of a simulated signal and thermal noise corresponding to 100h observations with the SKA. We then apply our inference pipeline to the latest measurements from the HERA interferometer. We report constraints on the X-ray emissivity, and confirm that cold reionization scenarios are unlikely to accurately represent our Universe.
In this work, we analytically calculate the spectra of primordial perturbations at the end of Type III hilltop inflation models under the slow-roll approximation. We examine the one-loop corrections of the spectra and find that those from the inflaton self-interaction are negligible. On the contrary, the loop effects from the interaction between the inflaton field and the waterfall field can be significant when the vacuum expectation value of the waterfall field is small. The implications are discussed.
According to a standard initial mass function, stars in the range 7-12 Msun constitute ~50% (by number) of the stars more massive than ~7 Msun, but, in spite of this, their evolutionary properties, and in particular their final fate, are still scarcely studied. In this paper we present a detailed study of the evolutionary properties of solar metallicity, non rotating stars in the range 7-15 Msun, from the pre main sequence phase up to the presupernova stage or up to an advanced stage of the thermally pulsing phase, depending on the initial mass. We find that (1) the 7.00 Msun develops a degenerate CO core and evolves as a classical AGB star in the sense that it does not ignite the C burning reactions; (2) stars with the initial mass M >= 9.22 Msun end their life as core collapse supernovae; (3) stars in the range 7.50 <= M/Msun <= 9.20 develop a degenerate ONeMg core and evolve through the thermally pulsing SAGB phase; 4) stars in the mass range 7.50 <= M/Msun <= 8.00 end their life as hybrid CO/ONeMg- or ONeMg- WD; (5) stars with the initial mass in the range 8.50 <= M/Msun <= 9.20 may potentially explode as electron capture supernovae.
Among the binary sources of interest for LISA some are quasimonochromatic, in the sense that the change in the gravitational wave frequency $\Delta f\lesssim 1\;\mbox{yr}^{-1}$ during the observation time. We study these sources in the frequency domain (FD) by taking into account the Doppler shift induced by LISA's motion and the LISA pattern functions, and we compare our results with previous work in the time domain. We also discuss the transition from the quasimonochromatic case to the stationary phase approximation commonly used in Fisher matrix calculations, which applies when $\Delta f\gtrsim 1\;\mbox{yr}^{-1}$.
Narrow Fe K$\alpha$ fluorescent emission lines arising at $\sim$kpc-scale separations from the nucleus have only been detected in a few AGN. The detections require that the extended line emission be spatially resolved and sufficiently bright. Compared to narrow Fe K$\alpha$ lines arising closer to the nucleus, they have much lower fluxes but show substantially larger equivalent widths, EW$_{\rm Fe K\alpha}$. We show that, in the optically-thin limit, a purely analytical argument naturally predicts large, EW$_{\rm FeK\alpha}\sim$1 keV, values for such lines, regardless of the details of equivalent hydrogen column density, $N_H$, or reprocessor geometry. Monte Carlo simulations corroborate this result and show that the simple analytic EW$_{\rm FeK\alpha}$ prescription holds up to higher $N_H$ approaching the Compton-thick regime. We compare to $Chandra$ observations from the literature and discuss that our results are consistent with the large EW$_{\rm FeK\alpha}$ values reported for local AGN, for which the line is detected in extended, up to $\sim$kpc-scale, regions. We argue that large EW$_{\rm FeK\alpha}$ from kpc-scale regions in AGN should be ubiquitous, because they do not depend on the absolute luminosity of the central X-ray source, and are measured only against the scattered continuum. We predict values to be of the order of $\sim$1 keV or larger, even for covering factors $\ll$1, and for arbitrarily small column densities. We propose that the large-scale molecular material that is now routinely being detected with the Atacama Large Millimeter/Submillimeter Array (ALMA) may act as an extended X-ray scattering reprocessor giving rise to $\sim$kpc-scale Fe K$\alpha$ emission.
We present here broadband radio observations of the CSM interacting SN2019oys. SN2019oys was first detected in the optical and was classified as a Type Ib SN. Then, about $\sim 100$ days after discovery, it showed an optical rebrightening and a spectral transition to a spectrum dominated by strong narrow emission lines, which suggests strong interaction with a distant, dense, CSM shell. We modeled the broadband, multi-epoch, radio spectra, covering 2.2 to 36 GHz and spanning from 22 to 1425 days after optical discovery, as a synchrotron emitting source. Using this modeling we characterized the shockwave and the mass-loss rate of the progenitor. Our broadband radio observations show strong synchrotron emission. This emission, as observed 201 and 221 days after optical discovery, exhibits signs of free-free absorption from the material in front of the shock traveling in the CSM. In addition, the steep power law of the optically thin regime points towards synchrotron cooling of the radiating electrons. Analyzing these spectra in the context of the SN-CSM interaction model gives a shock velocity of 14,000 $\rm km \, s^{-1}$, and an electron number density of $2.6 \times 10^5 \, \rm cm^{-3}$ at a distance of $2.6 \times 10^{16}$ cm. This translates to a high mass-loss rate from the progenitor massive star of $6.7 \times 10^{-4} \, \rm M_{\odot} yr^{-1}$ for an assumed wind of 100 $\rm km s^{-1}$ (assuming constant mass-loss rate in steady winds). The late-time radio spectra, 392 and 557 days after optical discovery, are showing broad spectral peaks. We show that this can be explained by introducing a non-homogeneous CSM structure.
Optical surveys have become increasingly adept at identifying candidate Tidal Disruption Events (TDEs) in large numbers, but classifying these generally requires extensive spectroscopic resources. We here present $\texttt{tdescore}$, a simple photometric classifier that is trained using a systematic census of $\sim$3000 nuclear transients from the Zwicky Transient Facility (ZTF). The sample is highly imbalanced, with TDEs representing $<$2% of the total. $\texttt{tdescore}$ is nonetheless able to reject non-TDEs with 99.6% accuracy, yielding a sample of probable TDEs with completeness of 77.0% and a purity of 80.3%. $\texttt{tdescore}$ is thus substantially better than any available TDE photometric classifier scheme in the literature, and performs comparably well to the single-epoch spectroscopy as a method for classifying ZTF nuclear transients, despite relying solely on ZTF data and multi-wavelength catalogue crossmatching. In a novel extension, we use 'SHapley Additive exPlanations' (SHAP) to provide a human-readable justification for each individual $\texttt{tdescore}$ classification, enabling users to understand and form opinions about the underlying classifier reasoning. $\texttt{tdescore}$ serves as a model for photometric identification of TDEs with time-domain surveys, such as the upcoming Rubin observatory.
We present photometry and spectroscopy of the slowly evolving superluminous Type IIn SN2015da. SN2015da is extraordinary for its very high peak luminosity, and also for sustaining a high luminosity for several years. Even at 8\,yr after explosion, SN2015da remains as luminous as the peak of a normal SNII-P. The total radiated energy integrated over this time period (with no bolometric correction) is at least 1.6 FOE. Including a mild bolometric correction, adding kinetic energy of the expanding cold dense shell of swept-up circumstellar material (CSM), and accounting for asymmetry, the total explosion kinetic energy was likely 5-10 FOE. Powering the light curve with CSM interaction requires an energetic explosion and 20 Msun of H-rich CSM, which in turn implies a massive progenitor system above 30 Msun. Narrow P Cyg features show steady CSM expansion at 90 km/s, requiring a high average mass-loss rate of roughly 0.1 Msun/yr sustained for 2 centuries before explosion (although ramping up toward explosion time). No current theoretical model for single-star pre-SN mass loss can account for this. The slow CSM, combined with broad wings of H$\alpha$ indicating H-rich material in the unshocked ejecta, disfavor a pulsational pair instability model for the pre-SN mass loss. Instead, violent pre-SN binary interaction is a likely cuprit. Finally, SN2015da exhibits the characteristic asymmetric blueshift in its emission lines from shortly after peak until the present epoch, adding another well-studied superluminous SNeIIn with unambiguous evidence of post-shock dust formation.
Large widefield surveys make possible the serendipitous discovery of rare sub-classes of pulsars. One such class are "spider"-type pulsar binaries, comprised of a pulsar in a compact orbit with a low-mass (sub)stellar companion. In a search for circularly-polarized radio sources in ASKAP Pilot Survey observations, we discovered highly variable and circularly polarized emission from a radio source within the error region of the $\gamma$-ray source {4FGL}~J1646.5$-$4406. The variability is consistent with the eclipse of a compact, steep-spectrum source behind ablated material from a companion in a $\sim 5.3\,$h binary orbit. Based on the eclipse properties and spatial coincidence with {4FGL} J1646.5$-$4406, we argue that the source is likely a recycled pulsar in a "redback" binary system. Using properties of the eclipses from ASKAP and Murchison Widefield Array observations, we provide broad constraints on the properties of the eclipse medium. We identified a potential optical/infra-red counterpart in archival data consistent with a variable low-mass star. Using the Parkes Radio Telescope "Murriyang" and MeerKAT, we searched extensively for radio pulsations but yielded no viable detections of pulsed emission. We suggest that the non-detection of pulses is due to scattering in the intra-binary material, but scattering from the ISM can also plausibly explain the pulse non-detections if the interstellar dispersion measure exceeds $\sim$600$\,$pc$\,$cm$^{-3}$. Orbital constraints derived from optical observations of the counterpart would be highly valuable for future $\gamma$-ray pulsation searches, which may confirm the source nature as a pulsar.
LAMOST J024048.51+195226.9 (J0240+1952) was recently identified as the second AE Aquarii (AE Aqr)-type cataclysmic variable, possessing the fastest known rotating white dwarf. We performed a Very Long Baseline Interferometry (VLBI) observation of J0240+1952 utilizing the European VLBI Network at 1.7\,GHz, to obtain the first view of the radio morphology on mas scale. Our high-resolution VLBI image clearly shows that the radio emission is compact on mas scale ($\lesssim2$\,AU), with no evidence for a radio jet or extended emission. The compact radio source has an average flux density of $\sim0.37$\,mJy, and its brightness temperature is given at $\gtrsim2.3\times10^{7}$\,K, confirming a non-thermal origin. The emission exhibits irregular variations on a time-scale of tens of minutes, similar to the radio flares seen in AE Aqr. The measured VLBI position of J0240+1952 is consistent with that derived from \textit{Gaia}. Our results favour the model in which the radio emission is attributed to a superposition of synchrotron radiation from expanding magnetized blobs of this system.
In this work, we made an extensive study about the possible presence of anisotropies in strange stars. To accomplish this task, we use three different configurations for the strange matter: the unpaired matter, a two-flavor super-conducting (2SC) strange matter, and a fully three-flavor super-conducting strange matter (CFL). For each configuration, we calculate the relevant quantities for the strange stars, such as the mass-radius relation, the dimensionless tidal parameter, the moment of inertia, and the surface curvature for different degrees of anisotropies. Whenever possible, we compare our results with constraints found in the literature, especially focusing on the existence of very massive pulsars (PSR J0952-0607), as well as very light compact objects (HESS J1731-347).
The unification of quantum mechanics and general relativity has long been elusive. Only recently have empirical predictions of various possible theories of quantum gravity been put to test. The dawn of multi-messenger high-energy astrophysics has been tremendously beneficial, as it allows us to study particles with much higher energies and travelling much longer distances than possible in terrestrial experiments, but more progress is needed on several fronts. A thorough appraisal of current strategies and experimental frameworks, regarding quantum gravity phenomenology, is provided here. Our aim is twofold: a description of tentative multimessenger explorations, plus a focus on future detection experiments. As the outlook of the network of researchers that formed through the COST Action CA18108 "Quantum gravity phenomenology in the multi-messenger approach (QG-MM)", in this work we give an overview of the desiderata that future theoretical frameworks, observational facilities, and data-sharing policies should satisfy in order to advance the cause of quantum gravity phenomenology.
The tidal acceleration experienced by an object at the event horizon of a black hole decreases as one over the square of the black hole's mass. As such there is a maximum mass at which a black hole can tidally disrupt an object outside of its event horizon and potentially produce observable emission. This maximum mass is known as the ``Hills mass'', and in full general relativity is a function of both the black hole's spin $a_\bullet$ and the inclination angle of the incoming object's orbit with respect to the black hole's spin axis $\psi$. In this paper we demonstrate that the Hills mass can be represented by a simple analytical function of $a_\bullet$ and $\psi$, the first general solution of this problem. This general solution is found by utilising the symmetries of a class of critical Kerr metric orbits known as the innermost bound spherical orbits. Interestingly, at fixed black hole spin the maximum Hills mass can lie at incoming orbital inclinations outside of the black hole's equatorial plane $\psi \neq \pi/2$. When compared to previous results in the literature this effect can lead to an increase in the maximum Hills mass (at fixed spin) by as much as a factor of $\sqrt{11/5} \simeq 1.48$ for a maximally rotating black hole. We then demonstrate how Bayesian inference, coupled with an estimate of the mass of a black hole in a tidal disruption event, can be used to place conservative constraints on that black hole's spin. We provide a publicly available code tidalspin which computes these spin distributions.
Thanks to the mid-infrared capability offered by Spitzer, systematic searches of dust in SNe have been carried out over the past decade. Studies have revealed the presence of a substantial amount of dust over a broad range of SN subtypes. How normal SNe present mid-IR excess at later time and turn out to be dusty SNe can be affected by several factors, such as mass-loss history and envelope structure of progenitors and their explosion environment. All these can be combined and related to their environmental properties. A systematic analysis of SNe that exploded under a dusty environment could be of critical importance to measure the properties of the dust-veiled exploding stars, and whether such an intense dust production process is associated with the local environment. In this work, we firstly use the IFS data to study the environmental properties of dusty SNe compared to those of normal ones, and analyze correlations between the environmental properties and their dust parameters. We find that dusty SNe have a larger proportion located at higher SFR regions compared to the normal types. The occurrence of dusty SNe is less dependent on metallicity, with the oxygen abundance spanning from subsolar to oversolar metallicity. We also find the host extinction of dusty SNe scatters a lot, with about 40% of dusty SN located at extremely low extinction environments, and another 30% of them with considerably high host extinction of E(B-V)>0.6 mag.
In this study, we explore the dynamical stability of magnetar bursts within the context of the chaos-randomness phase space for the first time, aiming to uncover unique behaviors compared to various astrophysical transients, including fast radio bursts (FRBs). We analyze burst energy time series data from active magnetar sources SGR J1550-5418 and SGR J1935+2154, focusing on burst arrival time and energy differences between consecutive events. We find a distinct separation in the time domain, where magnetar bursts exhibit significantly lower randomness compared to FRBs, solar flares, and earthquakes, with a slightly higher degree of chaos. In the energy domain, magnetar bursts exhibit a broad consistency with other phenomena, primarily due to the wide distribution of chaos-randomness observed across different bursts and sources. Intriguingly, contrary to expectations from the FRB-magnetar connection, the arrival time patterns of magnetar bursts in our analysis do not exhibit significant proximity to repeating FRBs in the chaos-randomness plane. This finding may challenge the hypothesis that FRBs are associated with typical magnetar bursts but indirectly supports the evidence that FRBs may primarily be linked to special magnetar bursts like peculiar X-ray bursts from SGR J1935+2154 observed simultaneously with Galactic FRB 200428.
We present the results of a joint {\it XMM-Newton} and {\it NuSTAR} campaign on the active galactic nucleus (AGN) IC 4329A, consisting of 9 $\times$ 20 ks {\it XMM-Newton} observations, and 5 $\times$ 20 ks {\it NuSTAR} observations within nine days, performed in August 2021. Within each observation, the AGN is not very variable, with the fractional variability never exceeding 5%. Flux variations are observed between the different observations, on timescales of days, with a 30% ratio between the minimum and the maximum 2-10 keV flux. These variations follow the softer-when-brighter behavior typically observed in AGN. In all observations, a soft excess is clearly present. Consistently with previous observations, the X-ray spectra of the source exhibit a cut-off energy between 140 and 250 keV, constant within the error in the different observations. We detected a narrow iron \ka line consistent with being constant during the monitoring, and likely originating in a distant neutral medium. We do not find evidence of a broad component of the iron line, suggesting that the inner disk does not produce strong reflection. We find that the reflection component is weak ($R_{\rm max}=0.009\pm0.002$). We also found the presence of a warm absorber component together with an ultra-fast outflow. Looking at their energetic, these outflows have enough mechanical power to exercise a significant feedback impact on the AGN surrounding environment.
Tidal disruption events (TDEs) provide a clue to the properties of a central supermassive black hole (SMBH) and an accretion disk around it, and to the stellar density and velocity distributions in the nuclear star cluster surrounding the SMBH. Deviations of TDE light curves from the standard occurring at a parabolic encounter with the SMBH depends on whether the stellar orbit is hyperbolic or eccentric (Hayasaki et al. 2018) and the penetration factor ($\beta$, tidal disruption radius to orbital pericenter ratio). We study the orbital parameters of bound and unbound stars being tidally disrupted by comparison of direct $N$-body simulation data with an analytical model. Starting from the classical steady-state Fokker-Planck model of Cohn & Kulsrud (1978), we develop an analytical model of the number density distribution of those stars as a function of orbital eccentricity ($e$) and $\beta$. To do so fittings of the density and velocity distribution of the nuclear star cluster and of the energy distribution of tidally disrupted stars are required and obtained from $N$-body data. We confirm that most of the stars causing TDEs in a spherical nuclear star cluster originate from the full loss-cone region of phase space, derive analytical boundaries in eccentricity-$\beta$ space, and find them confirmed by $N$-body data. Since our limiting eccentricities are much smaller than critical eccentricities for full accretion or full escape of stellar debris, we conclude that those stars are only very marginally eccentric or hyperbolic, close to parabolic.
Emergent bulk properties of matter governed by the strong nuclear force give rise to physical phenomena across vastly different scales, ranging from the shape of atomic nuclei to the masses and radii of neutron stars. They can be accessed on Earth by measuring the spatial extent of the outer skin made of neutrons that characterises the surface of heavy nuclei. The isotope $^{208}$Pb, owing to its simple structure and neutron excess, has been in this context the target of many dedicated efforts. Here, we determine the neutron skin from measurements of particle distributions and their collective flow in $^{208}$Pb+$^{208}$Pb collisions at ultrarelativistic energy performed at the Large Hadron Collider, which are sensitive to the overall size of the colliding $^{208}$Pb ions. By means of state-of-the-art global analysis tools within the hydrodynamic model of heavy-ion collisions, we infer a neutron skin $\Delta r_{np}=0.217\pm0.058$ fm, consistent with nuclear theory predictions, and competitive in accuracy with a recent determination from parity-violating asymmetries in polarised electron scattering. We establish thus a new experimental method to systematically measure neutron distributions in the ground state of atomic nuclei.
TianQin and LISA are space-based laser interferometer gravitational wave (GW) detectors planned to be launched in the mid-2030s. Both detectors will detect low-frequency GWs around $10^{-2}\,{\rm Hz}$, however, TianQin is more sensitive to frequencies above this common sweet-spot while LISA is more sensitive to frequencies below $10^{-2}\,{\rm Hz}$. Therefore, TianQin and LISA will be able to detect the same sources but with different accuracy depending on the source and its parameters. We consider some of the most important astrophysical sources -- massive black hole binaries, stellar-mass black hole binaries, double white dwarfs, extreme mass ratio inspirals, light and heavy intermediate mass ratio inspirals, as well as the stochastic gravitational background of astrophysical origin -- that TianQin and LISA will be able to detect. For each of these sources, we analyze how far they can be detected (detection distance) and how well their parameters can be measured (detection accuracy) using a Fisher Matrix analysis. We compare the results obtained by the three detection scenarios (TianQin alone, LISA alone, and joint detection by LISA and TianQin) highlighting the gains from joint detection as well as the contribution of TianQin and LISA to a combined study of astrophysical sources. In particular, we consider the different orientations, lifetimes, and duty cycles of the two detectors to explore how they can give a more complete picture when working together.
We present a phenomenological model to investigate the chiral phase transition characterized by parity doubling in dense, beta equilibrated, cold matter. Our model incorporates effective interactions constrained by SU(3) relations and considers baryonic degrees of freedom. By constraining the model with astrophysical data and nuclear matter properties, we find a first-order phase transition within realistic values of the slope parameter L. The inclusion of the baryon octet and negative parity partners, along with a chiral-invariant mass $m_{0}$, allows for a non-massless chiral symmetric phase. Through exploration of parameter space, we identify parameter sets satisfying mass and radius constraints without requiring a partonic phase. The appearance of the parity partner of the nucleon, the N(1535) resonance, suppresses strangeness, pushing hyperonization to higher densities. We observe a mild first-order phase transition to the chirally restored phase, governed by $m_{0}$. Our calculations of surface tension highlight its strong dependence on $m_{0}$. The existence of mixed phases is ruled out since they become energetically too costly. We compare stars with metastable and stable cores using both branches of the equation of state. Despite limited lifespans due to low surface tension values, phase conversion and star contraction could impact neutron stars with masses around 1.3 solar masses or more. We discuss some applications of this model in its non-zero temperatures generalization and scenarios beyond beta equilibrium that can provide insights into core-collapse supernovae, proto-neutron star evolution, and neutron star mergers. Core-collapse supernovae dynamics, influenced by chiral symmetry restoration and exotic hadronic states, affect explosion mechanisms and nucleosynthesis.
We introduce \texttt{Weakhub}, a novel neutrino microphysics library that provides opacities and kernels beyond conventional interactions used in the literature. This library includes neutrino-matter, neutrino-photon, and neutrino-neutrino interactions, along with corresponding weak and strong corrections. A full kinematics approach is adopted for the calculations of $\beta$-processes, incorporating various weak corrections and medium modifications due to the nuclear equation of state. Calculations of plasma processes, electron neutrino-antineutrino annihilation, and nuclear de-excitation are included. We also present the detailed derivations of weak interactions and the coupling of them to the two-moment based general-relativistic multi-group radiation transport in the \texttt{G}eneral-relativistic \texttt{mu}ltigrid \texttt{nu}merical (\texttt{Gmunu}) code. We compare the neutrino opacity spectra for all interactions and estimate their contributions at hydrodynamical points in core-collapse supernova and binary neutron star postmerger remnant, and predict the effects of improved opacities in comparison to conventional ones for a binary neutron star postmerger at a specific hydrodynamical point. We test the implementation of the conventional set of interactions by comparing it to an open-source neutrino library \texttt{NuLib} in a core-collapse supernova simulation. We demonstrate good agreement with discrepancies of less than $\sim 10\%$ in luminosity for all neutrino species, while also highlighting the reasons contributing to the differences. To compare the advanced interactions to the conventional set in core-collapse supernova modelling, we perform simulations to analyze their impacts on neutrino signatures, hydrodynamical behaviors, and shock dynamics, showing significant deviations.
We conduct one-dimensional stellar evolutionary numerical simulations under the assumption that an efficient dynamo operates in the core of massive stars years to months before core collapse and find that the magnetic activity enhances mass loss rate and might trigger binary interaction that leads to outbursts. We assume that the magnetic flux tubes that the dynamo forms in the inner core buoy out to the outer core where there is a steep entropy rise and a molecular weight drop. There the magnetic fields turn to thermal energy, i.e., by reconnection. We simulate this energy deposition where the entropy steeply rises and find that for our simulated cases the envelope radius increases by a factor of 1.2-2 and luminosity by about an order of magnitude. These changes enhance the mass loss rate. The envelope expansion can trigger a binary interaction that powers an outburst. Because magnetic field amplification depends positively on the core rotation rate and operates in cycles, not in all cases the magnetic activity will be powerful enough to change envelope properties. Namely, only a fraction of core-collapse supernovae experiences pre-explosion outbursts.
Fueling and feedback couple supermassive black holes (SMBHs) to their host galaxies across many orders of magnitude in spatial and temporal scales, making this problem notoriously challenging to simulate. We use a multi-zone computational method based on the general relativistic magneto-hydrodynamic (GRMHD) code KHARMA that allows us to span $7$ orders of magnitude in spatial scale, to simulate accretion onto a non-spinning SMBH from an external medium with Bondi radius $R_B\approx 2\times 10^5 \,GM_\bullet/c^2$, where $M_\bullet$ is the SMBH mass. For the classic idealized Bondi problem, spherical gas accretion without magnetic fields, our simulation results agree very well with the general relativistic analytic solution. Meanwhile, when the accreting gas is magnetized, the SMBH magnetosphere becomes saturated with a strong magnetic field. The density profile varies as $\sim r^{-1}$ rather than $r^{-3/2}$ and the accretion rate $\dot{M}$ is consequently suppressed by over 2 orders of magnitude below the Bondi rate $\dot{M}_B$. We find continuous energy feedback from the accretion flow to the external medium at a level of $\sim10^{-2}\dot{M}c^2 \sim 5\times 10^{-5} \dot{M}_B c^2$. Energy transport across these widely disparate scales occurs via turbulent convection triggered by magnetic field reconnection near the SMBH. Thus, strong magnetic fields that accumulate on horizon scales transform the flow dynamics far from the SMBH and naturally explain observed extremely low accretion rates compared to the Bondi rate, as well as at least part of the energy feedback.
We study the statistics of dwarf galaxy populations as a function of environment in the cold dark matter (CDM) and warm dark matter (WDM) cosmogonies, using hydrodynamical simulations starting from initial conditions with matched phases but differing power spectra, and evolved with the EAGLE galaxy formation model. We measure the abundance of dwarf galaxies within 3~Mpc of DM haloes with a present-day halo mass similar to that of the Milky Way (MW), and find that the radial distribution of galaxies $M_{*}>10^7$\Msun is nearly identical for WDM and CDM. However, the cumulative mass function becomes shallower for WDM at lower masses, yielding 50~per~cent fewer dwarf galaxies of $M_{*}\gtrsim10^{5}$~\Msun than CDM. The suppression of low-mass halo counts in WDM relative to CDM increases significantly from high-density regions to low-density regions for haloes in the region of the half-mode mass, $M_\rm{hm}$. The luminous fraction in the two models also diverges from the overdense to the underdense regions for $M>2M_\rm{hm}$, as the increased collapse delay at small densities pushes the collapse to after the reionization threshold. However, the stellar mass--halo mass relation of WDM haloes relative to CDM increases towards lower-density regions. Finally, we conclude that the suppression of galaxies with $M_{*}\gtrsim10^5$\Msun between WDM and CDM is independent of density: the suppression of halo counts and the luminous fraction is balanced by an enhancement in stellar mass--halo mass relation.
Blue supergiants (BSGs) are key objects for understanding the evolution of massive stars. However, discrepancies between theoretical predictions and empirical observations have opened up important questions yet to be answered. Studying statistically significant and unbiased samples of BSGs can help to improve the situation. We aim to perform a homogeneous and comprehensive quantitative spectroscopic analysis of a large sample of Galactic BSGs from the IACOB spectroscopic database. We derive the projected rotational velocity ($v\sin i$) and macroturbulent broadening ($v_{\rm mac}$) using IACOB-BROAD. We used FASTWIND computations to derive effective temperatures ($T_{\rm eff}$), surface gravities, microturbulences ($\xi$), Si and He surface abundances, and the wind-strength parameter. We provide estimates of the above-mentioned quantities for the largest sample of Galactic BSGs spectroscopically analyzed to date, comprising 538 O9-B5 type stars. We find a drop in the relative number of BSGs at ~21 kK, coinciding with a scarcity of fast rotating stars below that temperature. We speculate that this feature might be roughly delineating the location of the empirical Terminal-Age-Main-Sequence in the 15-85Msol range. By investigating the main characteristics of the $v\sin i$ distribution of O stars and BSGs as a function of $T_{\rm eff}$, we propose that an efficient mechanism transporting angular momentum from the stellar core to the surface might be operating along the main sequence. We find correlations between $\xi$, $v_{\rm mac}$ and the spectroscopic luminosity. We also find that no more than 20% of the BSGs have atmospheres enriched in He, and suggest that the origin of this specific sub-sample of BSGs might be in binary evolution. We do not find clear empirical evidence of an increase in the wind-strength over the wind bi-stability region towards lower $T_{\rm eff}$.
We investigate the effect of radial anisotropy on the rate of orbit decay in parabolic encounters of identical spherical galaxies. Our galaxy models have Hernquist density profiles and Osipkov--Merritt velocity distributions. We find that radially anisotropic models merge in as little as half the time of their isotropic counterparts. Anisotropic models are more susceptible to tidal deformation; this accelerates the transfer of orbital angular momentum to internal degrees of freedom. Even during the initial approach, the anisotropic models become more distorted, and arrive at pericentre already having lost substantial amounts of angular momentum. Our results may have implications for estimates of merger rates and persistence of tidal tails.
Extremely metal-poor galaxies (XMPGs) at relatively low redshift are excellent laboratories for studying galaxy formation and evolution in the early universe. Much effort has been spent on identifying them from large-scale spectroscopic surveys or spectroscopic follow-up observations. Previous work has identified a few hundred XMPGs. In this work, we obtain a large sample of 223 XMPGs at $z<1$ from the early data of the Dark Energy Spectroscopic Instrument (DESI). The oxygen abundance is determined using the direct $T_{\rm e}$ method based on the detection of the [O III]$\lambda$4363 line. The sample includes 95 confirmed XMPGs based on the oxygen abundance uncertainty; remaining 128 galaxies are regarded as XMPG candidates. These XMPGs are only 0.01% of the total DESI observed galaxies. Their coordinates and other proprieties are provided in the paper. The most XMPG has an oxygen abundance of $\sim 1/34 Z_{\odot}$, stellar mass of about $1.5\times10^7 M_{\odot}$ and star formation rate of 0.22 $M_{\odot}$ yr$^{-1}$. The two most XMPGs present distinct morphologies suggesting different formation mechanisms. The local environmental investigation shows that XMPGs preferentially reside in relatively low-density regions. Many of them fall below the stellar mass-metallicity relations (MZRs) of normal star-forming galaxies. From a comparison of the MZR with theoretical simulations, it appears that XMPGs are good analogs to high-redshift star-forming galaxies. The nature of these XMPG populations will be further investigated in detail with larger and more complete samples from the on-going DESI survey.
We present a new total intensity image of M31 at 1.248 GHz, observed with the Five-hundred-meter Aperture Spherical radio telescope (FAST) with an angular resolution of 4 arcmin and a sensitivity of about 16 mK. The new FAST image clearly reveals weak emission outside the ring due to its high sensitivity on large-scale structures. We derive a scale length of 2.7 kpc for the cosmic ray electrons and find that the cosmic ray electrons propagate mainly through diffusion by comparing the scale length at 4.8 GHz. The spectral index of the total intensity varies along the ring, which can be attributed to the variation of the spectra of synchrotron emission. This variation is likely caused by the change of star formation rates along the ring. We find that the azimuthal profile of the non-thermal emission can be interpreted by an axisymmetric large-scale magnetic field with varying pitch angle along the ring, indicating a complicated magnetic field configuration in M31.
Low Surface Brightness Galaxies (LSBGs) are excellent probes of quenching and other environmental processes near massive galaxies. We study an extensive sample of LSBGs near massive hosts in the local universe that are distributed across a diverse range of environments. The LSBGs with surface-brightness $\mu_{\rm eff,g}> $24.2 mag arcsec$^{-2}$ are drawn from the Dark Energy Survey Year 3 catalog while the hosts with masses $9.0< log(M_{\star}/M_{\odot})< 11.0$ comparable to the Milky Way and the Large Magellanic Cloud are selected from the z0MGS sample. We study the projected radial density profiles of LSBGs as a function of their color and surface brightness around hosts in both the rich Fornax-Eridanus cluster environment and the low-density field. We detect an overdensity with respect to the background density, out to 2.5 times the virial radius for both hosts in the cluster environment and the isolated field galaxies. When the LSBG sample is split by $g-i$ color or surface brightness $\mu_{\rm eff,g}$, we find the LSBGs closer to their hosts are significantly redder and brighter, like their high surface-brightness counterparts. The LSBGs form a clear 'red sequence' in both the cluster and isolated environments that is visible beyond the virial radius of the hosts. This suggests a pre-processing of infalling LSBGs and a quenched backsplash population around both host samples. However, the relative prominence of the 'blue cloud' feature implies that pre-processing is ongoing near the isolated hosts compared to the cluster hosts.
Massive, starbursting galaxies in the early Universe represent some of the most extreme objects in the study of galaxy evolution. One such source is HFLS3 (z~6.34), which was originally identified as an extreme starburst galaxy with mild gravitational magnification ($\mu$~2.2). Here, we present new observations of HFLS3 with the JWST/NIRSpec IFU in both low (PRISM/CLEAR; R~100) and high spectral resolution (G395H/290LP; R~2700), with high spatial resolution (~0.1") and sensitivity. Thanks to the combination of the NIRSpec data and a new lensing model with accurate spectroscopic redshifts, we find that the 3"x3" field is crowded, with a lensed arc (C, $z=6.3425\pm0.0002$), two galaxies to the south (S1 and S2, $z=6.3592\pm0.0001$), two galaxies to the west (W1, $z=6.3550\pm0.0001$; W2, $z=6.3628\pm0.0001$), and two low-redshift interlopers (G1, $z=3.4806\pm0.0001$; G2, $z=2.00\pm0.01$). We present spectral fits and morpho-kinematic maps for each bright emission line from the R2700 data for all sources except G2. From a line ratio analysis, the galaxies in component C are likely powered by star formation, while we cannot rule out or confirm the presence of AGN in the other high-redshift sources. We perform gravitational lens modelling, finding evidence for a two-source composition of the lensed central object and a comparable magnification factor ($\mu$=2.1-2.4) to previous work. The projected distances and velocity offsets of each galaxy suggest that they will merge within the next ~1Gyr. Finally, we examine the dust extinction-corrected SFR(Ha) of each z>6 source, finding that the total star formation ($510\pm140$Msol/yr, magnification-corrected) is distributed across the six z~6.34-6.36 objects over a region of diameter ~11kpc. Altogether, this suggests that HFLS3 is not a single starburst galaxy, but instead is a merging system of star-forming galaxies in the Epoch of Reionisation.
We developed and trained a pipeline of three machine learning (ML) models than can predict which sources are more likely to be an AGN and to be detected in specific radio surveys. Also, it can estimate redshift values for predicted radio-detectable AGNs. These models, which combine predictions from tree-based and gradient-boosting algorithms, have been trained with multi-wavelength data from near-infrared-selected sources in the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) Spring field. Training, testing, calibration, and validation were carried out in the HETDEX field. Further validation was performed on near-infrared-selected sources in the Stripe 82 field. In the HETDEX validation subset, our pipeline recovers 96% of the initially labelled AGNs and, from AGNs candidates, we recover 50% of previously detected radio sources. For Stripe 82, these numbers are 94% and 55%. Compared to random selection, these rates are two and four times better for HETDEX, and 1.2 and 12 times better for Stripe 82. The pipeline can also recover the redshift distribution of these sources with $\sigma_{\mathrm{NMAD}}$ = 0.07 for HETDEX ($\sigma_{\mathrm{NMAD}}$ = 0.09 for Stripe 82) and an outlier fraction of 19% (25% for Stripe 82), compatible with previous results based on broad-band photometry. Feature importance analysis stresses the relevance of near- and mid-infrared colours to select AGNs and identify their radio and redshift nature. Combining different algorithms in ML models shows an improvement in the prediction power of our pipeline over a random selection of sources. Tree-based ML models (in contrast to deep learning techniques) facilitate the analysis of the impact that features have on the predictions. This prediction can give insight into the potential physical interplay between the properties of radio AGNs (e.g. mass of black hole and accretion rate).
The solar physics community is entering a golden era that is ripe with next-generation ground- and space-based facilities. With ever-increasing resolving power stemming from the newest observational telescopes, it becomes more challenging to obtain (near-)simultaneous measurements at high spatial, temporal and spectral resolutions, while operating at the diffraction limit of these new facilities. Hence, in recent years there has been increased interest in the capabilities integral field units (IFUs) offer towards obtaining the trifecta of high spatial, temporal and spectral resolutions contemporaneously. To date, IFUs developed for solar physics research have focused on mid-optical and infrared measurements. Here, we present an IFU prototype that has been designed for operation within the near-ultraviolet to mid-optical wavelength range, hence providing additional spectral coverage to the instrument suites developed to date. The IFU was constructed as a low-budget proof-of-concept for the upcoming 2m class Indian National Large Solar Telescope and employs circular cross-section fibres to guide light into a Czerny-Turner configuration spectrograph, with the resulting spectra captured using a high quantum efficiency scientific CMOS camera. Mapping of each input fibre allows for the reconstruction of two-dimensional spectral images, with frame rates exceeding 20 per second possible while operating in a non-polarimetric configuration. The science verification data presented here highlights the suitability of fibre-fed IFUs operating at near-ultraviolet wavelengths for solar physics research. Importantly, the successful demonstration of this type of instrument paves the way for further technological developments to make a future variant suitable for upcoming ground-based and space-borne telescope facilities.
Next-generation radio interferometers like the Square Kilometer Array have the potential to unlock scientific discoveries thanks to their unprecedented angular resolution and sensitivity. One key to unlocking their potential resides in handling the deluge and complexity of incoming data. This challenge requires building radio interferometric imaging methods that can cope with the massive data sizes and provide high-quality image reconstructions with uncertainty quantification (UQ). This work proposes a method coined QuantifAI to address UQ in radio-interferometric imaging with data-driven (learned) priors for high-dimensional settings. Our model, rooted in the Bayesian framework, uses a physically motivated model for the likelihood. The model exploits a data-driven convex prior, which can encode complex information learned implicitly from simulations and guarantee the log-concavity of the posterior. We leverage probability concentration phenomena of high-dimensional log-concave posteriors that let us obtain information about the posterior, avoiding MCMC sampling techniques. We rely on convex optimisation methods to compute the MAP estimation, which is known to be faster and better scale with dimension than MCMC sampling strategies. Our method allows us to compute local credible intervals, i.e., Bayesian error bars, and perform hypothesis testing of structure on the reconstructed image. In addition, we propose a novel blazing-fast method to compute pixel-wise uncertainties at different scales. We demonstrate our method by reconstructing radio-interferometric images in a simulated setting and carrying out fast and scalable UQ, which we validate with MCMC sampling. Our method shows an improved image quality and more meaningful uncertainties than the benchmark method based on a sparsity-promoting prior. QuantifAI's source code: https://github.com/astro-informatics/QuantifAI.
The classification of variable objects provides insight into a wide variety of astrophysics ranging from stellar interiors to galactic nuclei. The Zwicky Transient Facility (ZTF) provides time series observations that record the variability of more than a billion sources. The scale of these data necessitates automated approaches to make a thorough analysis. Building on previous work, this paper reports the results of the ZTF Source Classification Project (SCoPe), which trains neural network and XGBoost machine learning (ML) algorithms to perform dichotomous classification of variable ZTF sources using a manually constructed training set containing 170,632 light curves. We find that several classifiers achieve high precision and recall scores, suggesting the reliability of their predictions for 112,476,749 light curves across 40 ZTF fields. We also identify the most important features for XGB classification and compare the performance of the two ML algorithms, finding a pattern of higher precision among XGB classifiers. The resulting classification catalog is available to the public, and the software developed for SCoPe is open-source and adaptable to future time-domain surveys.
Photonic lanterns (PLs) are tapered waveguides that gradually transition from a multi-mode fiber geometry to a bundle of single-mode fibers. In astronomical applications, PLs can efficiently couple multi-mode telescope light into a multi-mode fiber entrance and convert it into multiple single-mode beams. The output beams are highly stable and suitable for feeding into high-resolution spectrographs or photonic chip beam combiners. For instance, by using relative intensities in the output cores as a function of wavelength, PLs can enable spectroastrometry. In addition, by interfering beams in the output cores with a beam combiner in the backend, PLs can be used for high-throughput interferometric imaging. When used on an Extremely Large Telescope (ELT), with its increased sensitivity and angular resolution, the imaging and spectroastrometric capabilities of PLs will be extended to higher contrast and smaller angular scales. We study the potential spectroastrometry and imaging science cases of PLs on ELTs, including study of exomoons, broad-line regions of quasars, and inner circumstellar disks.
Some ideas for onboard processing and photometry with an astrometry satellite are presented, especially designed for GaiaNIR which may be launched about 2045 as a successor of Gaia. - Increased sensitivity, reduced image overlap, and simpler PSF calibration in GaiaNIR will result if the proposed initial processing of data from the detectors is implemented, because the across-scan smearing will become insignificant. - Filter photometry is required for high angular resolution as needed for astrometric and astrophysical reasons. Low-dispersion spectra are questioned because they fail at high star density. This will be a much greater problem with GaiaNIR than it is with Gaia because of the larger number of stars expected. It was the aim to collect in this note all arguments about GaiaNIR photometry which can be stated with words only, in correspondence with readers. The remaining work to be done for the definition of photometric equipment on the satellite requires further quantitative assessments and comparison of various options. Finally, 1) an advantage of filters is that the photometric observations can also be used for astrometry, 2) the XP spectra in Gaia will give very good astrophysical data for about 400 million single stars with G <~ 18.5 mag, but filters would have been better for all fainter and for all multiple stars, and 3) it is presently not clear which advantages for astrophysics low-dispersion spectra in the NIR might have over filters.
In July 2021, faculty from the UAlbany Department of Physics participated in a week-long field expedition with the organization UAPx to collect data on UAPs in Avalon, California, located on Catalina Island, and nearby. This paper reviews both the hardware and software techniques which this collaboration employed, and contains a frank discussion of the successes and failures, with a section about how to apply lessons learned to future expeditions. Both observable-light and infrared cameras were deployed, as well as sensors for other (non-EM) emissions. A pixel-subtraction method was augmented with other similarly simple methods to provide initial identification of objects in the sky and/or the sea crossing the cameras' fields of view. The first results will be presented based upon approximately one hour in total of triggered visible/night-vision-mode video and over 600 hours of untriggered (far) IR video recorded, as well as 55 hours of (background) radiation measurements. Following multiple explanatory resolutions of several ambiguities that were potentially anomalous at first, we focus on the primary remaining ambiguity captured at approximately 4am Pacific Time on Friday, July 16: a dark spot in the visible/near-IR camera possibly coincident with ionizing radiation that has thus far resisted a prosaic explanation. We conclude with quantitative suggestions for serious researchers in this still-nascent field of hard-science-based UAP studies, with an ultimate goal of identifying UAPs without confirmation bias toward either mundane or speculative conclusions.
The extreme sensitivity required for direct observation of gravitational waves by the Advanced LIGO detectors means that environmental noise can potentially contaminate gravitational wave signals. Consequently, environmental monitoring efforts have been undertaken and novel noise mitigation techniques have been developed which have helped keep environmental artifacts from influencing gravitational wave detections for the $90$ gravitational wave events detected from 2015--2020 by the aLIGO detectors. The increasing rate of gravitational wave detections due to detector sensitivity improvements requires sophisticated, reliable and automated ways to monitor and assess the degree of environmental coupling between gravitational wave detectors and their surroundings. We introduce a computational tool, PEMcheck, for quantifying the degree of environmental coupling present in gravitational wave signals using data from the network of environmental monitoring sensors. We study its performance when applied to the $79$ gravitational waves confidently detected in LIGO's third observing run and test its performance in the case of extreme environmental contamination of gravitational wave data. We find that PEMcheck's automated analysis identifies only a small number of gravitational waves that merit further study by environmental noise experts due to possible contamination, a substantial improvement over the manual vetting that occurred for every gravitational wave candidate in previous observing runs. Overall, PEMcheck works as intended. Consequently, PEMcheck will play a critical role in event validation during LIGO's fourth observing run.
The future Laser Interferometer Space Antenna (LISA) mission, which has successfully passed Mission Formulation phase, is in planning to be launched in 2030s. One of the ubiquitous LISA sources are the white-dwarf binaries (WDB) of which $\sim$40 are guaranteed sources as of now, making LISA unique in comparison to its ground-based counterpart. The current hardware design in planning necessitates a thorough check to determine whether the various locking schemes influence the guaranteed sources' signals significantly in order to re-consider that what is hard-coded in the phasemeter before launch for pre-science operations phase. Comparison of the phasemeter output of a face-on (V407Vul) binary and an edge-on (ZTFJ1539) binary indicates that the non-swap locking scheme, N2a, is optimal for instrument calibration. Additionally, the influence of the $\sim 7$ min orbital period edge-on source in two of the locking schemes yields a difference of maximum $\leq 10\%$ at the Time Delay Interferometry (TDI) output for data stream of one day. Simplified analyses show that neither of the locking schemes is favoured in the post-processing level. We find similar amplitudes in the TDI output stream for the face-on system V407Vul and the edge-on system ZTFJ1539 which leads to a significantly smaller inclination bias for the non-swap locking scheme. Additionally, a larger amplitude for edge-on systems will benefit most verification systems as the population of verification systems is biased towards edge-on systems as they are easier to detect in electromagnetic data.
Bubbles of nothing (BoNs) describe the decay of spacetimes with compact dimensions and are thus of fundamental importance for many higher dimensional theories proposed beyond the Standard Model. BoNs admit a 4-dimensional description in terms of a singular Coleman-de Luccia (CdL) instanton involving the size modulus field, stabilized by some potential $V(\phi)$. Using the so-called tunneling potential ($V_t$) approach, we study which types of BoNs are possible and for which potentials $V(\phi)$ can they be present. We identify four different types of BoN, characterized by different asymptotic behaviours at the BoN core and corresponding to different classes of higher dimensional theories, which we also classify. Combining numerous analytical and numerical examples, we study the interplay of BoN decays with other standard decay channels, identify the possible types of quenching of BoN decays and show how BoNs for flux compactifications can also be described in 4 dimensions by a multifield $V_t$. The use of the $V_t$ approach greatly aids our analyses and offers a very simple picture of BoNs which are treated in the same language as any other standard vacuum decays.
This paper investigates charged particle deflection in a Kerr spacetime background with a dipole magnetic field, focusing on the equatorial plane and employing the weak field approximation. We employ the Jacobi-Randers metric to unify the treatment of the gravitational and electromagnetic effects on charged particles. Furthermore, we utilize the Gauss-Bonnet theorem to calculate the deflection angle through curvature integrals. The difference between the prograde and retrograde deflection angles is linked to the non-reversibility of metrics and geodesics in Finsler geometry, revealing that this difference can be considered a Finslerian effect. We analyze the impact of both gravitomagnetic field and dipole magnetic field on particle motion and deflection using the Jacobi-Randers magnetic field. The model considered in this paper exhibits interesting features in the second-order approximation of ($M/b$). When $q\mu=2MaE$, the Jacobi-Randers metric possesses reversible geodesics, leading to equal prograde and retrograde deflection angles. In this case, the gravitomagnetic field and dipole magnetic field cancel each other out, distinguishing it from scenarios involving only the gravitomagnetic field or the dipole magnetic field. We also explore the magnetic field's impact on gravitational lensing of charged particles.
We construct and analyze a model of neutron star in the $\kappa$-deformed space-time. This is done by first deriving $\kappa$-deformed generalization of the Einstein tensor, starting from the non-commutative generalization of the metric tensor. By generalising the energy-momentum tensor to the non-commutative space-time and exploiting the $\kappa$-deformed dispersion relation, we then set up Einstein's field equations in the $\kappa$-deformed space-time. As we adopt a realisation of the non-commutative coordinates in terms of the commutative coordinates and their derivatives, our model is constructed in terms of commutative variables. Now by treating the interior of the star to be a perfect fluid as in the commutative space-time, we investigate the modification to neutron star's mass due to non-commutativity of the space-time. We show that the non-commutativity of the space-time enhances the mass limit of the neutron star. Using recent observational limit on the upper bound on the mass of neutron stars, we find the deformation parameter to be $|a|\sim 10^{-44}m$.
The group $SO(d+1,1)$ makes an appearance both as the conformal group of Euclidean space in $d$ dimensions and as the isometry group of de Sitter spacetime in $d+1$ dimensions. While this common feature can be taken as a hint towards holography on de Sitter space, understanding the representation theory has importance for cosmological applications where de Sitter spacetime is relevant. Among the categories of $SO(d+1,1)$ unitary irreducible representations, discrete series is important in physical applications because they are expected to capture gauge fields. However, they are also the most difficult ones to recognize in field theoretical examples compared to representations from the other categories. Here we point towards some examples where we are able to recognize discrete series representations from fields on de Sitter and highlight some of the properties of these representations.
This article is devoted to investigate the effects of $f(R)$ theory in the dynamics of a massless particle near the horizon of a static spherically symmetric (SSS) black hole. Deriving the equations of motion within $f(R)$ gravitational theories, novel solutions for charged and neutral black holes are obtained, introducing a dimensional parameter $a$ in $f(R)=R-2a\sqrt{R}$. Departing from General Relativity, these solutions showcase unique properties reliant on the dynamics of Ricci scalar. Analysis shows that chaos manifests within a specific energy range, with $a$ playing a crucial role. The study underscores the general applicability of the spherically symmetric metric, revealing insights into particle dynamics near black hole horizons. Despite an initially integrable nature, the introduction of harmonic perturbation leads to chaos, aligning with the Kolmogorov-Arnold-Moser theory. This research contributes to a nuanced understanding of black hole dynamics, emphasizing the importance of alternative theories of gravity.
The Schwarzschild geometry is investigated within the context of effective-field-theory models of gravity. Starting from its harmonic-coordinate expression, we derive the metric in standard coordinates by keeping the leading one-loop quantum contributions in their most general form. We examine the metric horizons and the nature of the hypersurfaces having constant radius; furthermore, a possible energy-extraction process which violates the null energy condition is described, and both timelike and null geodesics are studied. Our analysis shows that there is no choice of the sign of the constant parameter embodying the quantum correction to the metric which leaves all the features of the classical Schwarzschild solution almost unaffected.
The dependence of the parameters of the evolution of scalarly charged Black Holes (BHs) in the early Universe on the parameters of field-theoretic theories of interaction, the influence of the geometric factor of the structure of the relative position of BHs on the limitation of their maximum mass are studied, the problem of the metric of a scalarly charged BH in a medium of expanding scalarly charged matter is discussed, the expression is obtained for the macroscopic cosmological constant at late stages of expansion, generated by quadratic fluctuations of the metric, connecting the value of the cosmological constant with the BH masses and their concentration. Keywords: scalarly charged plasma, cosmological model, Higgs scalar field, gravitational stability, spherical perturbations, black hole formation, effective cosmological constant.
We show under the assumption of small data that solutions to the Einstein-Vlasov system with a pure magnetic field and Bianchi I symmetry isotropise and tend to dust solutions. We also obtain the decay rates for the main variables. This generalises part of the work [V.~G.~LeBlanc, Classical Quantum Gravity 14, 2281-2301 (1997)] concerning the future behaviour of orthogonal perfect fluids with a linear equation of state in the presence of a magnetic field to the Vlasov case.
We show that vacuum type N Kundt spacetimes in an arbitrary dimension admit a Kerr-Schild (KS) double copy. This is mostly done in a coordinate-independent way using the higher-dimensional Newman-Penrose formalism. We also discuss two kinds of non-uniqueness of an electromagnetic field corresponding to a given KS metric (i.e., its single copy) - these originate, respectively, from the rescaling freedom in the KS vector and from the non-uniqueness of the splitting of the KS metric in the flat part and the KS part. In connection to this, we show that the subset of KS pp-waves admits both null and non-null electromagnetic single copies. Since vacuum type N Kundt spacetimes are universal solutions of virtually any higher-order gravities and null fields in such backgrounds are immune to higher-order electromagnetic corrections, the KS-Kundt double copy demonstrated in the present paper also applies to large classes of modified theories.
We determine the adiabatic tidal contributions to the radiation reacted momentum impulse $\Delta p_i^\mu$ and scattering angle $\theta$ between two scattered massive bodies (neutron stars) at next-to-next-to-leading post-Minkowskian (PM) order. The state-of-the-art three-loop (4PM) worldline quantum field theory toolkit using dimensional regularization is employed to establish the classical observables. We encounter divergent terms in the gravito-electric and gravito-magnetic quadrupolar sectors necessitating the addition of post-adiabatic counterterms in this classical theory. This leads us to include also the leading post-adiabatic tidal contributions to the observables. The resulting renormalization group flow of the associated post-adiabatic Love numbers is established and shown to agree with a recent gravito-electric third post-Newtonian analysis in the non-relativistic limit.
During a binary black hole merger, multiple intermediary marginally outer trapped tubes connect the initial pair of apparent horizons with the final (single) apparent horizon. The marginally outer trapped surfaces (MOTSs) that foliate these tubes can have complicated geometries as well as non-spherical topologies. In particular, toroidal MOTSs form inside both of the original black holes during the early stages of a head-on merger that starts from time-symmetric initial data [1]. We show that toroidal MOTSs also form in the maximal analytic extension of the Schwarzschild spacetime as Kruskal time advances from the $T=0$ moment of time symmetry. As for the merger simulations, they cross the Einstein-Rosen bridge and are tightly sandwiched between the apparent horizons in the two asymptotic regions at early times. This strongly suggests that their formation is a consequence of the initial conditions rather than merger physics. Finally, we consider MOTSs of spherical topology in the Kruskal-Szekeres slicing and study their properties. All of these are contained within the apparent horizon but some do not enclose the wormhole.
We consider gauge theories on Poisson manifolds emerging as semiclassical approximations of noncommutative spacetime with Lie algebra type noncommutativity. We prove an important identity, which allows to obtain simple and manifestly gauge-covariant expressions for the Euler-Lagrange equations of motion, the Bianchi and the Noether identities. We discuss the non-Lagrangian equations of motion, and apply our findings to the $\kappa$-Minkowski case. We construct a family of exact solutions of the deformed Maxwell equations in the vacuum. In the classical limit, these solutions recover plane waves with left-handed and right-handed circular polarization, being classical counterparts of photons. The deformed dispersion relation appears to be nontrivial.
In this work, we consider a gravitational wave interacting with a quantum harmonic oscillator in the transverse-traceless gauge. We take the gravitational wave to be carrying the signatures of both plus and cross polarization at first. We then try to obtain a suitable form of the Lewis invariant using the most general form possible while considering only quadratic order contributions from both position and momentum variables. In order to progress further, we then drop the cross terms obtaining a separable Hamiltonian in terms of the first and the second spatial coordinates. We then obtain two Lewis invariants corresponding to each separable parts of the entire Hamiltonian of the system. Using both Lewis invariants, one can obtain two Ermakov-Pinney equations, from which we finally obtain the corresponding Lewis phase and eventually the Berry phase for the entire system. Finally, we obtain some explicit expressions of the Berry phase for a plane polarized gravitational wave with different choices of the harmonic oscillator frequency.
Most of the existing proposals for laboratory tests of a quantum nature of gravity are based on the use of two delocalized masses or harmonically bound masses prepared in pure quantum states with large enough spatial extent. Here a setup is proposed that is based on a single delocalized mass coupled to a harmonically trapped test mass (undergoing first expansion and then compression) that moves under the action of gravity. We investigate the in-principle feasibility of such an experiment, which turns out to crucially depend on the ability to tame Casimir-Polder forces. We thus proceed with a design aimed at achieving this, trying at the same time to take advantage of these forces rather than only fighting them.
I derive a smooth and global Kruskal-Szekeres coordinate chart for the maximal extension of the non-extremal Reissner-Nordstr\"om geometry that provides a generalization to the standard inner and outer Kruskal-Szekeres coordinates. The Kruskal-Szekeres diagram associated to this coordinate chart, whose existence is an interesting fact in and on itself, provides a simple alternative with a transparent physical interpretation to the conformal diagram of the spacetime.
We present a new covariant method of construction of the (position space) propagators in the $N$-dimensional (Euclidean) anti-de Sitter background for any gravitational theory with the Lagrangian that is an analytic expression in the metric, curvature, and covariant derivative. We show that the propagators (in Landau gauge) for all such theories can be expressed using the heat kernels for scalars and symmetric transverse-traceless rank-2 tensors on the hyperbolic $N$-space. The latter heat kernels are constructed explicitly and shown to be directly related to the former if an improved bi-scalar representation is used. Our heat kernel approach is first tested on general relativity, where we find equivalent forms of the propagators. Then it is used to obtain explicit expressions for propagators for various higher-derivative as well as infinite-derivative/nonlocal theories of gravity. As a by-product, we also provide a new derivation of the equivalent action (correcting a mistake in the original derivation) and an extension of the quadratic action to arbitrary ${N\geq 3}$ dimensions.
In AdS/CFT, observables on the boundary are invariant under renormalization group (RG) flow in the bulk. In this paper, we study holographic entanglement entropy under bulk RG flow and find that it is indeed invariant. We focus on tree-level RG flow, where massive fields in a UV theory are integrated out to give the IR theory. We explicitly show that in several simple examples, holographic entanglement entropy calculated in the UV theory agrees with that calculated in the IR theory. Moreover, we give an argument for this agreement to hold for general tree-level RG flow. Along the way, we generalize the replica method of calculating holographic entanglement entropy to bulk theories that include matter fields with nonzero spin.
By considering the concept of the modified Chaplygin gas (MCG) as a single fluid model unifying dark energy and dark matter, we construct a static, spherically charged black hole (BH) solution in the framework of General Relativity. The $P-V$ criticality of the charged anti-de Sitter (AdS) BH with a surrounding MCG is explored in the context of the extended phase space, where the negative cosmological constant operates as a thermodynamical pressure. This critical behavior shows that the small/large BH phase transition is analogous to the van der Waals liquid/gas phase transition. Accordingly, along the $P-V$ phase spaces, we derive the BH equations of state and then numerically evaluate the corresponding critical quantities. Similarly, critical exponents are identified, along with outcomes demonstrating the scaling behavior of thermodynamic quantities near criticality into a universal class. The use of \emph{geometrothermodynamic} (GT) tools finally offers a new perspective on discovering the critical phase transition point. At this stage, we apply a class of GT tools, such as Weinhold, Ruppeiner, HPEM, and Quevedo classes I and II. The findings are therefore non-trivial, as each GT class metric captures at least either the physical limitation point or the phase transition critical point. Overall, this paper provides a detailed study of the critical behavior of the charged AdS BH with surrounding MCG.
We consider the bending angle of the trajectory of a photon incident from and deflected to infinity around a Reissner-Nordstr\"om black hole. We treat the bending angle as a function of the squared reciprocal of the impact parameter and the squared electric charge of the background normalized by the mass of the black hole. It is shown that the bending angle satisfies a system of two inhomogeneous linear partial differential equations with polynomial coefficients. This system can be understood as an isomonodromic deformation of the inhomogeneous Picard-Fuchs equation satisfied by the bending angle in the Schwarzschild spacetime, where the deformation parameter is identified as the background electric charge. Furthermore, the integrability condition for these equations is found to be a specific type of the Painlev\'e VI equation that allows an algebraic solution. We solve the differential equations both at the weak and strong deflection limits. In the weak deflection limit, the bending angle is expressed as a power series expansion in terms of the squared reciprocal of the impact parameter and we obtain the explicit full-order expression for the coefficients. In the strong deflection limit, we obtain the asymptotic form of the bending angle that consists of the divergent logarithmic term and the finite O(1) term supplemented by linear recurrence relations which enable us to straightforwardly derive higher order coefficients. In deriving these results, the isomonodromic property of the differential equations plays an important role. Lastly, we briefly discuss the applicability of our method to other types of spacetimes such as a spinning black hole.
We study the optical properties of spacetime around a novel regular black hole (BH) in general relativity (GR) coupled to nonlinear electrodynamics (NED), which is asymptotically flat. First, we study the angular velocity and Lyapunov exponent in unstable photon circular orbits in the novel spherically symmetric BH spacetime. Later, the rotating regular BH solution is obtained using the Newmann-Janis algorithm, and the event horizon properties of the BH are determined. We analyze the effective potential for the circular motion of photons in the spacetime of the novel rotating BH. Also, we analyze the photon sphere around the novel BH and its shadow using celestial coordinates. We obtain that an increase of the BH spin and charge as well as NED field nonlinearity parameters causes an increase in the distortion parameter of the BH shadow, while, the area of the shadow and its oblateness decrease. Moreover, we also obtain the constraint values for the BH charge and the nonlinearity parameters using Event Horizon Telescope data from shadow sizes of supermassive BHs Sgr A* and M87*. Finally, the emission rate of BH evaporation through Hawking radiation is also studied.
We find a new solution to calculate the relativistic orbital periastron advance of a test-body subject to a post-Newtonian (PN) central force field, for relativistic models and theories beyond Einstein. This analitycal formula has general validity that includes all the PN contributions to the dynamics and is useful for high-precision gravitational tests. The solution is directly applicable to corrective potentials of various forms, without the need for numerical integration. Later, we apply it to the Scalar Tensor Fourth Order Gravity (STFOG) and NonCommutative Geometry, providing corrections to the Newtonian potential of Yukawa-like form $V(r)=\alpha \frac{e^{-\beta r}}{r}$, and we conduct the first analysis involving all the PN terms for these theories. The same work is performed with a Schwarzschild geometry perturbed by a Quintessence Field, leading to a power-law potential $V(r)=\alpha_q {r}^q$. Finally, by using astrometric data of the Solar System planetary precessions and those of S2 Star around Sgr A*, we infer new theoretical constraints and improvements in the bounds for $\beta$. The resulting simulated orbits turn out to be compatible with General Relativity.
To maximize the discovery potential of high-energy colliders, experimental searches should be sensitive to unforeseen new physics scenarios. This goal has motivated the use of machine learning for unsupervised anomaly detection. In this paper, we introduce a new anomaly detection strategy called FORCE: factorized observables for regressing conditional expectations. Our approach is based on the inductive bias of factorization, which is the idea that the physics governing different energy scales can be treated as approximately independent. Assuming factorization holds separately for signal and background processes, the appearance of non-trivial correlations between low- and high-energy observables is a robust indicator of new physics. Under the most restrictive form of factorization, a machine-learned model trained to identify such correlations will in fact converge to the optimal new physics classifier. We test FORCE on a benchmark anomaly detection task for the Large Hadron Collider involving collimated sprays of particles called jets. By teasing out correlations between the kinematics and substructure of jets, our method can reliably extract percent-level signal fractions. This strategy for uncovering new physics adds to the growing toolbox of anomaly detection methods for collider physics with a complementary set of assumptions.
We introduce the first generative model trained on the JetClass dataset. Our model generates jets at the constituent level, and it is a permutation-equivariant continuous normalizing flow (CNF) trained with the flow matching technique. It is conditioned on the jet type, so that a single model can be used to generate the ten different jet types of JetClass. For the first time, we also introduce a generative model that goes beyond the kinematic features of jet constituents. The JetClass dataset includes more features, such as particle-ID and track impact parameter, and we demonstrate that our CNF can accurately model all of these additional features as well. Our generative model for JetClass expands on the versatility of existing jet generation techniques, enhancing their potential utility in high-energy physics research, and offering a more comprehensive understanding of the generated jets.
We set up a numerical S-matrix bootstrap problem to rigorously constrain bound state couplings given by the residues of poles in elastic amplitudes. We extract upper bounds on these couplings that follow purely from unitarity, crossing symmetry, and the Roy equations within their proven domain of validity. First we consider amplitudes with a single spin 0 or spin 2 bound state, both with or without a self-coupling. Subsequently we investigate amplitudes with the spectrum of bound states corresponding to the estimated glueball masses of pure SU(3) Yang-Mills. In the latter case the 'glue-hedron', the space of allowed couplings, provides a first-principles constraint for future lattice estimates.
We study the phenomenology of superheavy decaying dark matter with mass around $10^{10}$ GeV which can arise in the low-energy limit of string compactifications. Generic features of string theory setups (such as high scale supersymmetry breaking and epochs of early matter domination driven by string moduli) can accommodate superheavy dark matter with the correct relic abundance. In addition, stringy instantons induce tiny $R$-parity violating couplings which make dark matter unstable with a lifetime well above the age of the Universe. Adopting a model-independent approach, we compute the flux and spectrum of high-energy gamma rays and neutrinos from three-body decays of superheavy dark matter and constrain its mass-lifetime plane with current observations and future experiments. We show that these bounds have only a mild dependence on the exact nature of neutralino dark matter and its decay channels. Applying these constraints to an explicit string model sets an upper bound of ${\cal O}(0.1)$ on the string coupling, ensuring that the effective field theory is in the perturbative regime.
Kinetic mixing of the dark photon, the gauge boson of a hidden $U(1)_D$, with the Standard Model (SM) gauge fields to induce an interaction between ordinary matter and dark matter (DM) at 1-loop requires the existence of portal matter (PM) fields having both dark and SM charges. As discussed in earlier work, these same PM fields can also lead to other loop-level mechanisms besides kinetic mixing that can generate significant interactions between SM fermions and the dark photon in a manner analogous to those that can be generated between a Dirac neutrino and a SM photon, \ie, dark moments. In either case, there are reasons to believe, \eg, due to the RGE running of the $U(1)_D$ gauge coupling, that PM fields may have $\sim$ TeV-scale masses that lie at or above those directly accessible to the HL-LHC. If they lie above the reach of the HL-LHC, then the only way to possibly explore the physics at this high scale is via indirect measurements made at lower energies, \eg, at lepton colliders operating in the $m_Z$ to 1 TeV range. In particular, processes such as $e^+e^- \to \gamma+$DM or $e^+e^-\to \bar ff$, where $f$ is a SM fermion, may be most useful in this regard. Here we explore these possibilities within the framework of a simple toy PM model, introduced in earlier work, based on a non-abelian dark gauge group completion operating at the PM scale. In the KM setup, we show these efforts fail due to the inherently tiny cross sections in the face of substantial SM backgrounds. However, in the case of interactions via induced dark moments, since they necessarily take the form of higher dimensional operators whose influence grows with energy, we show that access to PM-scale information may become possible for certain ranges of the toy model parameters for both of these $e^+e^-$ processes at a 1 TeV collider.
Analytical expressions for the Transverse Momentum Dependent (TMD, or unintegrated) gluon and sea quark densities in nuclei are derived at leading order of QCD running coupling. The calculations are performed in the framework of the rescaling model and Kimber-Martin-Ryskin (KMR) prescription, where the Bessel-inspired behavior of parton densities at small Bjorken $x$ values, obtained in the case of flat initial conditions in the double scaling QCD approximation, is applied. The derived expressions are used to evaluate the inclusive heavy flavor production in proton-lead collisions at the LHC. We find a good agreement of our results with latest experimental data collected by the CMS and ALICE Collaborations at $\sqrt s = 5.02$ GeV.
We propose a method to compute the effective potential of QCD from gap equations by introducing the homotopy transformation between solutions of the equation of motion. Via this method, the effective potential can be obtained beyond the bare vertex approximation, which then generalizes the Cornwall, Jackiw and Tomboulis (CJT) effective potential for bilocal composite operators. Moreover, the extended effective potential is set to be a function of self energy instead of the composite operator, which is the key point to make the potential bounded from below as for the auxiliary field (AF) Potential. We then apply the effective potential in the coexistence region where there exists at least two solutions, for instance, in vacuum with small current quark mass, and the first order phase transition region in finite temperature and chemical potential, which provides the in-medium behavior of the latent heat and false vacuum energy.
In the past twenty years, many new hadrons that are difficult to be explained within the conventional quark model have been discovered in the quarkonium region, which are called exotic hadrons. Belle II experiment, as the next-generation $B$ factory, provides a good platform to explore them. The charmonium-like states can be produced at Belle II in several ways, such as $B$ meson decays, initial-state radiation processes, two-photon collisions, and double charmonium productions. The bottomonium-like states can be produced directly in $e^+e^-$ colliding energies at Belle II with low continuum backgrounds. Belle II plans to perform a high-statistics energy scan from the $B\bar B$ threshold up to the highest possible energy of 11.24 GeV to search for new $Y_b$ states with $J^{PC}$ = $1^{--}$, $X_b$ (the bottom counterpart of $\chi_{c1}(3872)$ (also known as $X(3872)$)), and partners of $Z_b$ states. In this paper, we give a mini-review on the status and prospects of exotic hadrons at Belle II.
A qualitative calculation and discussion of two vortex states collisions are given in the scalar $\phi ^4$ model. Three kinds of vortex states -- Bessel, general monochromatic, and Laguerre-Gaussian vortex states -- are considered. It is found that the total final momentum distribution in collision of physical vortex states displays general topological structures, which depend on the initial vortex states' topological charges, which are proportional to the orbital angular momenta. This peculiar matching provides a novel observable, the topological number of momentum distribution, and it may represent a new fascinating research direction in particle physics. We also find that the situation when the angular momenta of the two colliding Laguerre-Gaussian states combine to zero can be recognized by the total final momentum distribution close to the collision axis. Both features can be used to measure the orbital angular momentum of vortex states.
We study the jet quenching parameter $\hat q$ in a non-equilibrium plasma using the QCD effective kinetic theory. We discuss subleading terms at large jet momentum p, show that our expression for $\hat q$ reproduces thermal results at small and large transverse momentum cutoffs for infinite p and construct an interpolation between these limits to be used in phenomenological applications. Using simple non-equilibrium distributions that model pertinent features of the bottom-up thermalization scenario, we analytically assess how anisotropy, under- or overoccupation affect the jet quenching parameter. Our work provides more details on the $\hat q$ formula used in our preceding work [arXiv:2303.12595] and sets the stage for further numerical studies of jet momentum broadening in the initial stages of heavy-ion collisions from QCD kinetic theory.
The upper limit of the branching ratio of the rare $\Xi_b^-\rightarrow \Xi^- \gamma$ decay is obtained as $BR(\Xi_b^-\rightarrow \Xi^- \gamma)<1.3\times10^{-4}$ by the LHCb. In the present work we study this decay within the light cone QCD sum rules employing the $\Xi_b$ distribution amplitudes. At first stage, the form factors entering the $\Xi_b^-\rightarrow \Xi^- \gamma$ decay are obtained. Next, using the results for the form factors the corresponding branching ratio for this decay is estimated to be $BR(\Xi_b^-\rightarrow \Xi^- \gamma)=(4.8\pm 1.3)\times 10^{-5}$. This value lies below the upper limit established by the LHCb collaboration. Our finding for the branching ratio is also compared with the results of the other theoretical approaches existing in the literature.
We present an evaluation of the glueball spectrum for configurations produced with $N_f=1$ dynamical fermions as a function of the $m_{\rm PCAC}$ mass. We obtained masses of states that fall into the irreducible representations of the octahedral group of rotations in combination with the quantum numbers of charge conjugation $C$ and parity $P$. Due to the low signal to noise ratio, practically, we can only extract masses for the irreducible representations $R^{PC}=$ $A_1^{++}$, $E^{++}$, $T_2^{++}$ as well as $A_1^{-+}$. We make use of the Generalized Eigenvalue Problem (GEVP) with an operator basis consisting only of gluonic operators. Throughout this work we are aiming towards the identification of the effects of light dynamical quarks on the glueball spectrum and how this compares to the statistically more precise spectrum of SU(3) pure gauge theory. We used large gauge ensembles which consist of ${\sim {~\cal O}}(10 {\rm K})$ configurations. Our findings demonstrate that the low-lying spectrum of the scalar, tensor as well as pseudo-scalar glueballs receive negligible contributions from the inclusion of $N_f=1$ dynamical fermions.
Though the Standard Model (SM) provides a very elegant description of the interactions among fundamental particles, there are ample evidences suggesting that new physics is needed. In particular, extending the scalar sector has enough motivation from vacuum stability, electroweak phase transition and various other sectors. Among different such extensions, the two-Higgs-doublet model (THDM) is the simplest one that preserves the electroweak $\rho$ parameter. Flavour-changing neutral currents (FCNC) are usually avoided by implementing additional discrete symmetries, but this type of models are subject to severe phenomenological constraints. In the more general framework of the aligned THDM (ATHDM) tree-level FCNCs are avoided by choosing the same flavour structure for the Yukawa couplings of the two scalar doublets, which results in weaker phenomenological constraints. Here, we present a global fit of the ATHDM, using the package HEPfit that performs a bayesian analysis on the parameter-space of this model with the help of stability and perturbativity bounds, experimental data for various flavour and electroweak precision observables, and constraints from Higgs searches at the LHC. This global fit has been performed assuming that all additional scalars are heavier than the SM Higgs and that there are no extra sources of CP violation beyond the CKM phase.
Dilepton decays of pseudoscalar mesons have been drawing particular interest, thanks to their sensitivity to both the QCD dynamics at low energy and also signals beyond the Standard Model. In this context, we shortly review our recent study on an improved Standard-Model prediction for the rare decay $\pi^0\to e^+e^-$, and compare it with the first determination on the lattice that predicted also the $\pi^0\to \gamma\gamma $ decay width as a byproduct. In addition, we discuss our recent work on $K_L\to\ell^+\ell^-$ decays and its connection to lattice QCD. We comment on the current uncertainty estimates and discuss how they could be improved profiting from future experiments and progress in lattice QCD.
Recently, several models have been suggested to reduce the tension between Gallium and reactor antineutrino spectral ratio data which is found in the framework of 3+1 active-sterile neutrino mixing. Among these models, we consider the extensions of 3+1 mixing with a finite wavepacket size, or the decay of the heaviest neutrino $\nu_4$, or the possibility to have a broad $\nu_4$ mass distribution. We consider the reactor antineutrino rate data and we show that these models cannot liminate the tension between Gallium and reactor rate data that is found in the 3+1 neutrino mixing framework. Indeed, we show that the parameter goodness of fit remains small. We consider also a model which explains the Gallium Anomaly with non-standard decoherence in the framework of three-neutrino mixing. We find that it is compatible with the reactor rate data.
We study the inclusive and exclusive cross sections of $e^+e^-\to$ hadrons for center-of-mass energies between 3.70 GeV and 3.83 GeV to infer the mass, width, and couplings of the $\psi(3770)$ resonance. By using a coupled-channel $K$-matrix approach, we setup our analysis to respect unitarity and the analyticity properties of the underlying scattering amplitudes. We fit three models to the full dataset and identify our nominal results through a statistical model comparison. We find that, accounting for the interplay between the $\psi(2S)$ and the $\psi(3770)$, no further pole is required to describe the $\psi(3770)$ line shape. In particular we derive from the pole location $M_{\psi(3770)} = 3777.1 \pm 0.3$ MeV and $\Gamma_{\psi(3770)} = 26.3 \pm 0.8$ MeV. Moreover, we find the decay to $D^+D^-$ and $D^0\bar{D}^0$ to be consistent with isospin symmetry and obtain a sizeable branching ratio $\mathcal{B}(\psi(3770) \to \textrm{non-}D\bar{D}) = (22 \pm 6)\%$.
As the study of three-hadron physics from lattice QCD matures, it is necessary to develop proper analysis tools in order to reliably study a variety of phenomena, including resonance spectroscopy and nuclear structure. Reconstructing the three-particle scattering amplitude requires solving integral equations, which can be written in terms of data-constrained dynamical functions and physical on-shell quantities. The driving term in these equations is the so-called one-particle exchange, which leads to a kinematic divergence for particles on-mass-shell. A vital component in defining three-particle amplitudes with definite parity and total angular momentum, which are used in spectroscopic studies, is to project the one-particle exchange into definite partial waves. We present a general procedure to construct exact analytic partial wave projections of the one-particle exchange contribution for any system composed of three spinless hadrons. Our result allows one full control over the analytic structure of the projection, which we explore for some low-lying partial waves with applications to three pions.
Two-loop multi-leg form factors in off-shell kinematics require knowledge of planar and nonplanar double box Feynman diagrams with massless internal propagators. These are complicated functions of Mandelstam variables and external particle virtualities. The latter serve as regulators of infrared divergences, thus making these observables finite in four space-time dimensions. In this paper, we use the method of canonical differential equations for calculation of (non)planar double box integrals in the near mass-shell kinematical regime, i.e., where virtualities of external particles are much smaller than the Mandelstam variables involved. We deduce a basis of master integrals with uniform transcendental weight based on the analysis of leading singularities by means of the Baikov representation as well as an array of complementary techniques. We dub the former asymptotically canonical since it is valid in the near mass-shell limit of interest. We iteratively solve resulting differential equations up to weight four in terms of multiple polylogarithms.
Intensity Frontier experiments are often evaluated by the smallest coupling it can probe, irrespective of what particle can be found or the scientific significance of its detection. In this work, we propose a new framework that determines the number of events required to characterize new particle properties. For example, we show that Heavy Neutral Leptons require 100 events to establish the neutrino mass hierarchy, and 1000 events to reveal the Majorana phase of active neutrinos. Ultimately, this framework presents a more objective way to connect experiments to their scientific outcomes.
We investigate the thermodynamic geometry of the quark-meson model at finite temperature, $T$, and quark number chemical potential, $\mu$. We extend previous works by the inclusion of fluctuations exploiting the functional renormalization group approach. We use recent developments to recast the flow equation into the form of an advection-diffusion equation. We adopt the local potential approximation for the effective average action. We focus on the thermodynamic curvature, $R$, in the $(\mu,T)$ plane, in proximity of the chiral crossover, up to the critical point of the phase diagram. We find that the inclusion of fluctuations results in a smoother behavior of $R$ near the chiral crossover. Moreover, for small $\mu$, $R$ remains negative, signaling the fact that bosonic fluctuations reduce the capability of the system to completely overcome the fermionic statistical repulsion of the quarks. We investigate in more detail the small $\mu$ region by analyzing a system in which we artificially lower the pion mass, thus approaching the chiral limit in which the crossover is actually a second order phase transition. On the other hand, as $\mu$ is increased and the critical point is approached, we find that $R$ is enhanced and a sign change occurs, in agreement with mean field studies. Hence, we completely support the picture that $R$ is sensitive to a crossover and a phase transition, and provides information about the effective behavior of the system at the phase transition.
We present a new way to interpret Top Standard Model measurements going beyond the SMEFT framework. Instead of the usual paradigm in Top EFT, where the main effects come from tails in momenta distributions, we propose an interpretation in terms of new physics which only shows up at loop-level. The effects of these new states, which can be lighter than required within the SMEFT, appear as distinctive structures at high momenta, but may be suppressed at the tails of distributions. As an illustration of this phenomena, we present the explicit case of a UV model with a $\mathcal{Z}_2$ symmetry, including a Dark Matter candidate and a top-partner. This simple UV model reproduces the main features of this class of signatures, particularly a momentum-dependent form factor with more structure than the SMEFT. As the new states can be lighter than in SMEFT, we explore the interplay between the reinterpretation of direct searches for colored states and Dark Matter, and Top measurements made by ATLAS and CMS in the differential $t \bar t$ final state. We also compare our method with what one would expect using the SMEFT reinterpretation, finding that using the full loop information provides a better discriminating power.
Phase transitions are part of everyday life, yet are also believed to be part of the history of our universe, where the nature of particle interactions change as the universe settles into its vacuum state. The discovery of the Higgs, and measurement of its mass suggests that our vacuum may not be entirely stable, and that a further phase transition could take place. This article is based on a talk in the Oldenberg Series, and reviews how we find the probability of these phase transitions, discussing past work on how black holes can dramatically change the result! Apart from a brief update at the end, this article mostly follows the content of the talk.
Double beta decays are rare nuclear processes that can occur in two modes: two-neutrino double beta decay, observed in the Standard Model, and neutrinoless double beta decay, a hypothetical process with profound implications for Particle Physics. To draw reliable conclusions from their experimental constraints, it is necessary to have accurate predictions of the underlying hadronic interactions described by quantum chromodynamics (QCD), a non-Abelian gauge theory with the symmetry group SU(3). QCD predictions require non-perturbative methods for calculating observables, and lattice QCD (LQCD), a numerical method based on QCD formulated on a finite space-time grid, is the only reliable first-principles technique for obtaining quantitative results. However, LQCD needs formal prescriptions to match numerical results with observables. This thesis provides such prescriptions for double beta decays using the finite volume effects in the LQCD framework. Matching relations that connect two-nucleon double beta decay amplitudes to quantities accessible via LQCD calculations, namely the nuclear matrix elements and two-nucleon energy spectra in a finite volume are provided. The impact of uncertainties is examined on the precision with which low-energy constants of the corresponding effective field theories can be determined from future LQCD calculations. Hamiltonian simulation of QCD is another non-perturbative method of solving QCD which can be more suitable in some cases than the conventional LQCD. The rise of tensor network methods and quantum simulation has made Hamiltonian simulation of lattice gauge theories (LGTs) a reality. Towards the goal of simulating QCD, a loop-string-hadron (LSH) formulation of an SU(3) LGT with matter in 1+1 dimensions is developed in this thesis, motivated by recent studies that showed the LSH formulation of an SU(2) LGT to be advantageous over other formulations.
We examine the dynamics of quarks and gauge fields in the lowest energy states in the QED and QCD interactions by combining Schwinger's longitudinal confinement in (1+1)D with Polyakov's transverse confinement in (2+1)D in a ``stretch (2+1)D'' flux tube model in (3+1)D. For such QED and QED systems in the flux tube configuration with cylindrical symmetry, we separate out the transverse and longitudinal degrees of freedom, approximate the non-Abelian QCD in the quasi-Abelian approximation, and solve the derived equations to study the collective excitations. We find stable collective QED and QCD excitations showing up as confined QED and QCD mesons, in support of previous theoretical studies and recent observations of the anomalous hypothetical X17 and E38 particles. Future theoretical lattice gauge calculations of QED in (3+1)D with the inclusion of the Schwinger longitudinal confinement mechanism and experimental confirmation of the hypothetical X17 and E38 particles will shed definitive light on quark confinement in the QED interaction in (3+1)D.
We apply the joint threshold and transverse momentum dependent (TMD) factorization theorem to introduce new threshold-TMD distribution functions, including threshold-TMD parton distribution functions (PDFs) and fragmentation functions (FFs). We apply Soft-Collinear Effective Theory and renormalization group methods to carry out QCD evolution for both threshold-TMD PDFs and FFs. We show the universality of threshold-TMD functions among three standard processes, i.e. the Drell-Yan production in $pp$ collisions, semi-inclusive deep-inelastic scattering and back-to-back two hadron production in $e^+e^-$ collisions. In the end, we present the numerical predictions for different threshold-TMD functions and also transverse momentum distributions at $pp$, $ep$, and $e^+e^-$ collisions.
Quantum computers are expected to give major speed-ups for the simulation of quantum systems. In this work, we present quantum gates that simulate the colour part of the interactions of quarks and gluons in perturbative quantum chromodynamics (QCD). As a first application, we implement these circuits on a simulated noiseless quantum computer and use them to calculate colour factors for various examples of Feynman diagrams. This work constitutes a first key step towards a quantum simulation of generic scattering processes in perturbative QCD.
A number of extensions of the Standard Model predict the existence of new light, weakly-coupled particles that couple to the visible sector through higher-dimensional operators containing one or two photons, suppressed by a high new physics scale, and thus have long lifetimes. In this work, we study the prospects for detecting three $\sim\,$sub-GeV such long-lived particles (LLP) at intensity frontier experiments: a massive spin-2 mediator ($G$), a dark axion portal, and a light neutralino coupled to ALPino or gravitino. We consider the production and visible decays of these particles in several current and proposed beam dump experiments (CHARM, NuCal, SeaQuest, NA62, SHiP) as well as in the LHC detectors (FASER, FASER$\nu$, FLArE, MATHUSLA). In addition to the usual displaced vertex signature, we also examine the impact of electron scattering signature and the Primakoff-like process which leads to conversion of $G$ into a photon or to the secondary LLP production via upscattering of the lighter dark sector state on dense material put in front of the detector. In all cases, we find the SHiP experiment could provide the strongest constraints for the displaced vertex search, while FASER2/FPF could provide complementary coverage of the $\gamma c\tau \sim 1$ m decay length region of the parameter space.
The flavor-changing neutral current (FCNC) decays of charmed hadrons with missing energy $(\not\!\!E)$ can serve as potentially promising hunting grounds for hints of new physics, as the standard-model backgrounds are very suppressed. A few of such processes have been searched for in recent experiments, specifically $D^0\to\,\not\!\!E$ by Belle and $D^0\to\pi^0$$\not\!\!E$ and $\Lambda_c^+\to p\!\not\!\!E$ by BESIII, resulting in upper bounds on their branching fractions. We consider them to illuminate the possible contributions of the quark transition $c\to u\!\not\!\!E$ with a couple of invisible spinless bosons carrying away the missing energy, assuming that they are not charge conjugates of each other and hence can have unequal masses. We find that these data are complementary in that they constrain different sets of the underlying operators and do not cover the same ranges of the bosons' masses, but there are regions not yet accessible. From the allowed parameter space, we show that other $D$-meson decays, such as $D\to\rho$$\not\!\!E$, and the charmed-baryon ones $\Xi_c\to(\Sigma,\Lambda)$$\not\!\!E$ can have sizable branching fractions and therefore may offer further probes of the new-physics interactions. We point out the importance of $D^0\to\gamma\!\not\!\!E$ which are not yet searched for but could access parts of the parameter space beyond the reach of the other modes. In addition, we look at a scenario where the invisibles are instead fermionic, namely sterile neutrinos, and a scalar leptoquark mediates $c\to u\!\not\!\!E$. We discuss the implications of the aforesaid bounds for this model. The predictions we make for the various charmed-hadron decays in the different scenarios may be testable in the near future by BESIII and Belle II.
Statistical models serve as the cornerstone for hypothesis testing in empirical studies. This paper introduces a new cross-platform Python-based package designed to utilise different likelihood prescriptions via a flexible plug-in system. This framework empowers users to propose, examine, and publish new likelihood prescriptions without developing software infrastructure, ultimately unifying and generalising different ways of constructing likelihoods and employing them for hypothesis testing within a unified platform. We propose a new simplified likelihood prescription, surpassing previous approximation accuracies by incorporating asymmetric uncertainties. Moreover, our package facilitates the integration of various likelihood combination routines, thereby broadening the scope of independent studies through a meta-analysis. By remaining agnostic to the source of the likelihood prescription and the signal hypothesis generator, our platform allows for the seamless implementation of packages with different likelihood prescriptions, fostering compatibility and interoperability.
Hints for an additional Higgs boson with a mass of about 95 GeV originate from LEP and searches in the diphoton channel by CMS and ATLAS. A search for resonant production of SM plus BSM Higgs bosons in the diphoton plus bb channel by CMS showed some excess for a 650 GeV resonance decaying into the SM Higgs plus a 95 GeV Higgs boson. We investigate whether these phenomena can be interpreted simultaneously within the NMSSM subject to the latest constraints on couplings of the SM Higgs boson, on extra Higgs bosons from the LHC, and on dark matter direct detection cross sections. We find that the hints for a 95 GeV Higgs boson in the diphoton channel by CMS and ATLAS and in the diphoton plus bb channel by CMS can be fitted simultaneously within the 2 sigma level.
We consider magnetic Weyl semimetals. First of all we review relation of intrinsic anomalous Hall conductivity, band contribution to intrinsic magnetic moment, and the conductivity of chiral separation effect (CSE) to the topological invariants written in terms of the Wigner transformed Green functions (with effects of interaction and disorder taken into account). Next, we concentrate on the CSE. The corresponding bulk axial current would result in accumulation of particles and holes of opposite chiralities at the surface of the sample. However, this accumulation is compensated by the flow of the states in momentum space along the Fermi arcs. Together with the bulk CSE current this flow forms closed Weyl orbits. Their detection can be considered as experimental discovery of chiral separation effect. Previously it was proposed to detect Weyl orbits through the observation of quantum oscillations \cite{Potter_2014} . We propose the alternative way to detect existence of Weyl orbits through the observation of their contributions to Hall conductance.
I work out in full detail the derivation of the Cross-Free Family (CFF) representation for the box diagram, and highlight the differences with its Time Ordered Perturbation Theory (TOPT) representation. I briefly discuss implications for the threshold singularity structure of the diagram.
This work is based on the author's PhD thesis. The main result of the thesis is the use of the boost operator to develop a systematic method to construct new integrable spin chains with nearest-neighbour interaction and characterized by an R-matrix of non-difference form. This method has the advantage of being more feasible than directly solving the Yang-Baxter equation. We applied this approach to various contexts, in particular, in the realm of open quantum systems, we achieved the first classification of integrable Lindbladians. These operators describe the dynamics of physical systems in contact with a Markovian environment. Within this classification, we discovered a novel deformation of the Hubbard model spanning three sites of the spin chain. Additionally, we applied our method to classify models with $\mathfrak{su}(2)\oplus \mathfrak{su}(2)$ symmetry and we recovered the matrix part of the S-matrix of $AdS_5 \times S^5$ derived by requiring centrally extended $\mathfrak{su}(2|2)$ symmetry. Furthermore, we focus on spin 1/2 chain on models of 8-Vertex type and we showed that the models of this class satisfy the free fermion condition. This enables us to express the transfer matrix associated to some of the models in a diagonal form, simplifying the computation of the eigenvalues and eigenvectors. The thesis is based on the works: 2003.04332, 2010.11231, 2011.08217, 2101.08279, 2207.14193, 2301.01612, 2305.01922.
Castellano, Ruiz, and Valenzuela recently observed a remarkable "pattern" in infinite-distance limits of moduli spaces in quantum gravity, which relates the field space variation of the mass of the lightest tower of particles to the field space variation of the species scale. In this work, we show how a version of this pattern can be proven to hold for BPS particles and strings throughout the vector multiplet moduli space of a 5d supergravity theory, even in regions where the particle masses and string tensions are substantially modified relative to their asymptotic behavior in the infinite-distance limits. This suggests that a suitably defined version of the pattern may hold not merely in the asymptotic limits of moduli space, but in the interior as well.
We show that a class of $L$-loop conformal ladder graphs correspond to twisted partition functions of free massive complex scalars in $d=2L+1$ dimensions. The graphs arise as four-point functions in certain two- and four-dimensional conformal fishnet models. The twisted thermal two-point function of the scalars is a generator of such conformal graphs for all loops. We argue that this correspondence is seeded by a system of two decoupled harmonic oscillators twisted by an imaginary chemical potential. We find a number of algebraic and differential relations among the conformal graphs which mirror the underlying free dynamics.
This work aims to investigate the classical-level duality between the $SIM(1)$-Maxwell-Chern-Simons (MCS) model and its self-dual counterpart. Initially, our focus is on free-field cases to establish equivalence through two distinct approaches: comparing the equations of motion and utilizing the master Lagrangian method. In both instances, the classical correspondence between the self-dual field and the MCS dual field undergoes modifications due to very special relativity (VSR). Specifically, duality is established only when the associated VSR-mass parameters are the same. Furthermore, we analyze the duality when the self-dual model is minimally coupled to fermions. As a result, we show that Thirring-like interactions, corrected for non-local VSR contributions, are included in the MCS model. Additionally, we demonstrate the equivalence of the fermion sectors in both models, thereby concluding the proof of classical-level duality.
We study to what extent, and in what form, the notion of gauge-string duality may persist at finite $N$. It is shown, in the half-BPS sector, that the states of D3 giant graviton branes in $\mathrm{AdS}_5 \times S^5$ are holographically dual to certain auxiliary ghosts that compensate for finite $N$ trace relations in $U(N)$ $\mathcal{N}=4$ super Yang-Mills. The complex formed from spaces of states of bulk D3 giants is observed to furnish a BRST-like resolution of the half-BPS Hilbert space of $U(N)$ $\mathcal{N}=4$ SYM at finite $N$. We argue that the identification between the states of certain bulk D-branes and the auxiliary ghosts in the boundary holds rather generally at vanishing 't Hooft coupling $\lambda = 0$. We propose that a complex, which furnishes a BRST-like resolution of the finite $N$ Hilbert space of a boundary $U(N)$ gauge theory at $\lambda = 0$, should be identified as the space of states of the dual string theory in the $\alpha' \to \infty$ limit. The Lefschetz trace formula provides the holographic map in this regime, where bulk observables are computed by taking the alternating sum of the expectation values in an ensemble of states built on each open string vacuum. The giant graviton expansion is recovered and generalized in a limit of the resolution.
The Schur index of the Higgs branch of four-dimensional $\mathcal{N}=2$ SCFTs is related to the spectrum of non-unitary two-dimensional CFTs. The rank one case has been shown to lead to the non-unitary CFTs with Deligne-Cvitanovic (DC) exceptional sequence of Lie groups. We show that a subsequence $(A_0, A_{\frac{1}{2}}, A_1, A_2, D_4)$ within the non-unitary sequence is related to a subsequence in the Mathur-Mukhi-Sen (MMS) sequence of unitary theories. We show that 2D non-unitary $G_2$ theory is related to unitary $E_6$ theory, and using this result along with the Galois conjugation, we propose that the $G_2$ Higgs branch is a sub-branch of the $E_6$ Higgs branch.
Two-dimensional conformal field theory is a powerful tool to understand the geometry of surfaces. Here, we study Liouville conformal field theory in the classical (large central charge) limit, where it encodes the geometry of the moduli space of Riemann surfaces. Generalizing previous work, we employ this to study moduli spaces of higher genus surfaces, surfaces with boundaries, and surfaces with cone points. In each case, the knowledge of classical conformal blocks provides an extremely efficient approximation to the Weil-Petersson metric on moduli space. We find detailed agreement with analytic results for volumes and geodesic lengths on moduli space.
In the perturbative treatment of interacting quantum field theories, if the interaction Lagrangian changes adiabatically in a finite interval of time, secular growths may appear in the truncated perturbative series also when the interaction Lagrangian density is returned to be constant. If this happens, the perturbative approach does not furnish reliable results in the evaluation of scattering amplitudes or expectation values. In this paper we show that these effects can be avoided for adiabatically switched-on interactions, if the spatial support of the interaction is compact and if the background state is suitably chosen. We start considering equilibrium background states and show that, when thermalisation occurs (interaction Lagrangian of spatial compact support), secular effects are avoided. Furthermore, no secular effects pop up if the limit where the Lagrangian is supported everywhere in space is taken after thermalisation (large time limit), in contrast to the reversed order. This result is generalized showing that if the interaction Lagrangian is spatially compact, secular growths are avoided for generic background states which are only invariant under time translation and to states whose explicit dependence of time is not too strong. Finally, as an example, we apply the presented theorems to study a complex scalar and a Dirac field in a classical external electromagnetic potential, on a background KMS state, to manifest that a spatially compact supported interaction does not give rise to secular growths.
This is the author's PhD thesis. Two main sections address various aspects of mirror symmetry for compact Calabi-Yau threefolds and the roles that classically modular varieties play in string theory compactifications. The main results include a study, and finding an application to the higher genus problem, of infinite Coxeter symmetries in the sets of Gopakumar-Vafa invariants; provision of a new class of solutions to the supersymmetric flux vacuum equations which have elsewhere been conjectured to give weight-two modular manifolds; provision of two new conjectural examples of weight-four modular varieties (rank-two attractors); and discussion of a set of numerical relations between infinite sums of Gromov-Witten invariants and critical L-values.
We find all polynomial tau-functions of the n-th reduced BKP hierarchy (=n-th Sawada-Kotera hierarchy). The name comes from the fact that for n=3 the simplest equation of the hierarchy is the famous Sawada-Kotera equation.
We study the on-shell scattering amplitudes in quantum gravity for high-energy collisions in the eikonal approximation. We first evaluate the $n$-loop 2-particle scattering amplitude in the high energy and low momentum transfer limit. We do so in a symmetrized manner by finding the contributions of each of the particle worldines to the scattering amplitude and gluing them together via the $n$ intermediate particle exchanges. In this limit on applying the eikonal approximation and summing over all $n$-loop Feynman diagrams we obtain a closed form for the 2 particle scattering amplitude. Finally, we extend this approach to obtain a generalized eikonal approximation for $N$-particle scattering at high energies and small momentum transfers. The generalised form of the scattering amplitude can then be used to evaluate the bound states of the system.
Quantum teleportation can be used to define a notion of parallel transport which characterizes the entanglement structure of a quantum state \cite{Czech:2018kvg}. This suggests one can formulate a gauge theory of entanglement. In \cite{Wong:2022mnv}, it was explained that measurement based quantum computation in one dimension can be understood in term of such a gauge theory (MBQC). In this work, we give an alternative formulation of this "entanglement gauge theory" as an extended topological field theory. This formulation gives a alternative perspective on the relation between the circuit model and MBQC. In addition, it provides an interpretation of MBQC in terms of the extended Hilbert space construction in gauge theories, in which the entanglement edge modes play the role of the logical qubit.
We study the structure of wave functions in complex Chern-Simons theory on the complement of a hyperbolic knot, emphasizing the similarities with the topological string/spectral theory correspondence. We first conjecture a hidden integrality structure in the holomorphic blocks and show that this structure guarantees the cancellation of potential singularities in the full non-perturbative wave function at rational values of the coupling constant. We then develop various techniques to determine the wave function at such rational points. Finally, we illustrate our conjectures and obtain explicit results in the examples of the figure-eight and the three-twist knots. In the case of the figure-eight knot, we also perform a direct evaluation of the state integral in the rational case and observe that the resulting wave function has the features of the ground state for a quantum mirror curve.
Dissipation is a ubiquitous phenomenon in nature that affects the fate of chaotic quantum dynamics. To characterize the interplay between quantum chaos and dissipation in generic quantum many-body systems, we consider a minimal dissipative Floquet many-body system. We study the dissipative form factor (DFF), an extension of the spectral form factor to open quantum systems, of the random phase model in the presence of arbitrary one-site nonunitary gates (quantum channels). In the limit of large local Hilbert space dimension, we obtain an exact expression for the DFF averaged over the random unitary gates, with simple, closed-form expressions in the limit of large times. We find that, for long enough times, the system always relaxes (i.e., the DFF decays) with two distinctive regimes characterized by the presence or absence of gap closing. While the system can sustain a robust ramp for a long (but finite) time interval in the gap-closing regime, relaxation is ``assisted'' by quantum chaos in the regime where the gap remains nonzero. In the latter regime, we find that, if the thermodynamic limit is taken first, the gap does not close even in the dissipationless limit.
We investigate the twist-3 generalized parton distributions (GPDs) for the valence quarks of the proton within the basis light-front quantization (BLFQ) framework. We first solve for the mass spectra and light-front waved functions (LFWFs) in the leading Fock sector using an effective Hamiltonian. Using the LFWFs we then calculate the twist-3 GPDs via the overlap representation. By taking the forward limit, we also get the twist-3 parton distribution functions (PDFs), and discuss their properties. Our prediction for the twist-3 scalar PDF agrees well with the CLAS experimental extractions.
The following is a master thesis centered around the concept of localisation and the Third Way Theory. This thesis discusses various aspects of supersymmetric localisation in one and three dimensions, and contains original results with regards to the Third Way Theory. It starts off with the Witten index for a one-dimensional supersymmetric system and derives various aspects through localisation. After this, the thesis moves on to the Third Way Theory. First, it offers a review of the Third Way Theory, a deformation of topologically massive Yang-Mills theory in three dimensions. Then it moves on to original results. These include a supersymmetrisation of the Third Way Theory and consequently a localisation of the Third Way Theory, which is to say, a method of deriving non-perturbative results.
Matrix elements in different representations are connected by quadratic relations. If matrix elements are those of a $\textit{group element}$, i.e. satisfying the property $\Delta(X) = X\otimes X$, then their generating functions obey bilinear Hirota equations and hence are named $\tau$-functions. However, dealing with group elements is not always easy, especially for non-commutative algebras of functions, and this slows down the development of $\tau$-function theory and the study of integrability properties of non-perturbative functional integrals. A simple way out is to use arbitrary elements of the universal enveloping algebra, and not just the group elements. Then the Hirota equations appear to interrelate a whole system of generating functions, which one may call $\textit{generalized}$ $\tau$-functions. It was recently demonstrated that this idea can be applicable even to a somewhat sophisticated case of the quantum toroidal algebra. We consider a number of simpler examples, including ordinary and quantum groups, to explain how the method works and what kind of solutions one can obtain.
We introduce a new Lie-algebraic approach to explicitly construct the motivic coaction and single-valued map of multiple polylogarithms in any number of variables. In both cases, the appearance of multiple zeta values is controlled by conjugating generating series of polylogarithms with Lie-algebra generators associated with odd zeta values. Our reformulation of earlier constructions of coactions and single-valued polylogarithms preserves choices of fibration bases, exposes the correlation between multiple zeta values of different depths and paves the way for generalizations beyond genus zero.
Loop Vertex Expansion (LVE) was developed for the construction of QFT models with local and non-local interactions. Using LVE, one can prove the analyticity in the finite cardioid-like domain in the complex plain of the coupling constant of the free energies and cumulants of various vector, matrix, or tensor-type models. Here, applying the idea of choosing the initial approximation depending on the coupling constant, we construct the analytic continuation of the free energy of the quartic matrix model beyond the standard LVE cardioid over the branch cut and for arbitrary large couplings.
Using carefully chosen projections, we consider different Carroll limits of relativistic Dirac fermions in any spacetime dimensions. These limits define Carroll fermions of two types: electric and magnetic. The latter type transforms as a reducible but indecomposable representation of the Carroll group. We also build action principles for all Carroll fermions we introduce; in particular, in even dimensions we provide an action principle for a minimal magnetic Carroll fermion, having the same number of components as a Dirac spinor. We then explore the coupling of these fermions to magnetic Carroll gravity in both its first-order and second-order formulations.
We explore connections between three structures associated with the cohomology of the moduli of 1-dimensional stable sheaves on $\mathbb{P}^2$: perverse filtrations, tautological classes, and refined BPS invariants for local $\mathbb{P}^2$. We formulate the $P=C$ conjecture identifying the perverse filtration with the Chern filtration for the free part of the cohomology. This can be viewed as an analog of de Cataldo--Hausel--Migliorini's $P=W$ conjecture for Hitchin systems. Our conjecture is compatible with the enumerative invariants of local $\mathbb{P}^2$ calculated by refined Pandharipande--Thomas theory or Nekrasov partition functions. It provides a cohomological lift of a conjectural product formula of the asymptotic refined BPS invariants. We prove the $P=C$ conjecture for degrees $\leq 4$.
The tensor functor called $\alpha$-induction arises from a Frobenius algebra object, or a Q-system, in a braided unitary fusion category. In the operator algebraic language, it gives extensions of endomorphism of $N$ to $M$ arising from a subfactor $N\subset M$ of finite index and finite depth giving a braided fusion category of endomorpshisms of $N$. It is also understood in terms of Ocneanu's graphical calculus. We study this $\alpha$-induction for bi-unitary connections, which give a characterization of finite-dimensional nondegenerate commuting squares and gives certain 4-tensors appearing in recent studies of 2-dimensional topological order. We show that the resulting $\alpha$-induced bi-unitary connections are flat if we have a commutative Frobenius algebra, or a local Q-system. Examples related to chiral conformal field theory and the Dynkin diagrams are presented.
We extend the study of superinstantons presented in 1905.01513 to include orthosymplectic supergroup gauge theories, $B_{n_0|n_1}$, $C_n$, and $D_{n_0|n_1}$. We utilize equivariant localization to obtain the LMNS contour integral formula for the instanton partition function, and we investigate the Seiberg--Witten geometries associated with these theories. We also explore the brane configurations involving positive and negative branes together with O-planes that realize the orthosymplectic supergroup theories.
We study the renormalized action and the renormalized presymplectic potential for Maxwell fields on Anti de Sitter backgrounds of any dimensions. We then use these results to explicitly derive finite boundary charges for angle-dependent asymptotic symmetries. We consider both Poincar\'e and Bondi coordinates, the former allowing us to control the systematics for arbitrary $D$, the latter being better suited for a smooth flat limit.
We discuss the generalization of the local renormalization group approach to theories in which Weyl symmetry is gauged. These theories naturally correspond to scale invariant - rather than conformal invariant - models in the flat space limit. We argue that this generalization can be of use when discussing the issue of scale vs conformal invariance in quantum and statistical field theories. The application of Wess-Zumino consistency conditions constrains the form of the Weyl anomaly and the beta functions in a nonperturbative way. In this work we concentrate on two dimensional models including also the contributions of the boundary. Our findings suggest that the renormalization group flow between scale invariant theories differs from the one between conformal theories because of the presence of a new charge that appears in the anomaly. It does not seem to be possible to find a general scheme for which the new charge is zero, unless the theory is conformal in flat space. Two illustrative examples involving flat space's conformal and scale invariant models that do not allow for a naive application of the standard local treatment are given.
By incorporating two gauge connections, transgression forms provide a generalization of Chern-Simons actions that are genuinely gauge-invariant on bounded manifolds. In this work, we show that, when defined on a manifold with a boundary, the Hamiltonian formulation of a transgression field theory can be consistently carried out without the need to implement regularizing boundary terms at the level of first-class constraints. By considering boundary variations of the relevant functionals in the Poisson brackets, the surface integral in the very definition of a transgression action can be translated into boundary contributions in the generators of gauge transformations and diffeomorphisms. This prescription systematically leads to the corresponding surface charges of the theory, reducing to the general expression for conserved charges in (higher-dimensional) Chern-Simons theories when one of the gauge connections in the transgression form is set to zero.
Self-dual Yang-Mills theory admits an underlying infinite dimensional symmetry algebra, which has been obtained from mode expansion of Mellin transformed 4d scattering amplitudes and separately, Koszul duality on twistor space. In this paper, we propose to derive an explicit 2d realization of the algebra by performing a particular gauge transformation on the twistor action for self-dual Yang-Mills. The gauge parameter used in the transformation generates pure gauge connections corresponding to large gauge transformations on 4d Minkowski space, which localises part of the twistor action to a CP1 after scaling reduction of twistor space. Under a projection, it can be mapped to the celestial sphere at the light-cone cut of the origin on null infinity. Geometrically, this is the common boundary celestial sphere shared by Euclidean AdS3 or Lorentzian dS3 slices of Minkowski space. We comment on the geometric meaning of the derivation from the perspective of minitwistor spaces of the 3d slices embedded in 4d Minkowski space. Using the action functional of this 2d CFT, we compute its stress-energy tensor and central charge. By a further marginal deformation, we calculate correlation functions of current algebra generators purely from the 2d side which reproduce 4d MHV form factors.
We discuss two-dimensional conformal field theories (CFTs) which are invariant under gauging a non-invertible global symmetry. At every point on the orbifold branch of $c=1$ CFTs, it is known that the theory is self-dual under gauging a $\mathbb{Z}_2\times \mathbb{Z}_2$ symmetry, and has $\mathsf{Rep}(H_8)$ and $\mathsf{Rep}(D_8)$ fusion category symmetries as a result. We find that gauging the entire $\mathsf{Rep}(H_8)$ fusion category symmetry maps the orbifold theory at radius $R$ to that at radius $2/R$. At $R=\sqrt{2}$, which corresponds to two decoupled Ising CFTs (Ising$^2$ in short), the theory is self-dual under gauging the $\mathsf{Rep}(H_8)$ symmetry. This implies the existence of a topological defect line in the Ising$^2$ CFT obtained from half-space gauging of the $\mathsf{Rep}(H_8)$ symmetry, which commutes with the $c=1$ Virasoro algebra but does not preserve the fully extended chiral algebra. We bootstrap its action on the $c=1$ Virasoro primary operators, and find that there are no relevant or marginal operators preserving it. Mathematically, the new topological line combines with the $\mathsf{Rep}(H_8)$ symmetry to form a bigger fusion category which is a $\mathbb{Z}_2$-extension of $\mathsf{Rep}(H_8)$. We solve the pentagon equations including the additional topological line and find 8 solutions, where two of them are realized in the Ising$^2$ CFT. Finally, we show that the torus partition functions of the Monster$^2$ CFT and Ising$\times$Monster CFT are also invariant under gauging the $\mathsf{Rep}(H_8)$ symmetry.
We combine perturbation theory with analytic and numerical bootstrap techniques to study the critical point of the long-range Ising (LRI) model in two and three dimensions. This model interpolates between short-range Ising (SRI) and mean-field behaviour. We use the Lorentzian inversion formula to compute infinitely many three-loop corrections in the two-dimensional LRI near the mean-field end. We further exploit the exact OPE relations that follow from bulk locality of the LRI to compute infinitely many two-loop corrections near the mean-field end, as well as some one-loop corrections near SRI. By including such exact OPE relations in the crossing equations for LRI we set up a very constrained bootstrap problem, which we solve numerically using SDPB. We find a family of sharp kinks for two- and three-dimensional theories which compare favourably to perturbative predictions, as well as some Monte Carlo simulations for the two-dimensional LRI.
We present an overview of both older and recent developments concerning scale separation in string theory. We focus on parametric scale separation obtained at the classical level in flux compactifications down to AdS vacua. We review the scenarios that have been proposed to achieve a hierarchy of scales between spacetime and the internal space, built from a low-dimensional perspective. We then discuss how they have been understood to arise from proper higher-dimensional descriptions. Eventually, limitations of these constructions as well as Swampland and holographic arguments addressing the question of scale separation in string theory are discussed. The purpose of the review is to draw an accurate picture of the state of the art of the subject at the moment.
Topological entanglement entropy (TEE), the sub-leading term in the entanglement entropy of topological order, is the direct evidence of the long-range entanglement. While effective in characterizing topological orders on closed manifolds, TEE is model-dependent when entanglement cuts intersect with physical gapped boundaries. In this paper, we study the origin of this model-dependence by introducing a model-independent picture of partitioning the topological orders with gapped boundaries. In our picture, the entanglement boundaries (EBs), i.e. the virtual boundaries of each subsystem induced by the entanglement cuts, are assumed to be gapped boundaries with boundary defects. At this model-independent stage, there are two choices one has to make manually in defining the bi-partition: the boundary condition on the EBs, and the coherence between certain boundary states. We show that TEE appears because of a constraint on the defect configurations on the EBs, which is choice-dependent in the cases where the EBs touch gapped boundaries. This choice-dependence is known as the ambiguity in entanglement entropy. Different models intrinsically employ different choices, rendering TEE model-dependent. For Z2 toric code, the ambiguity can be fully characterized by two parameters that respectively quantifies the EB condition and the coherence. In particular, calculations compatible with the folding trick naturally choose EB conditions that respect electric-magnetic duality and set specific parameter values.
Classically, the ground states of $\mathcal{N}=4$ supersymmetric Yang-Mills theory on $\mathbb{R}\times S^3/\Gamma$ where $\Gamma$ is a discrete ADE subgroup of $SU(2)$ are represented by flat Wilson lines winding around the ADE singularity. By a duality relating such ground states to WZW conformal blocks, the ground state degeneracy cannot be lifted by quantum corrections. Using the superconformal index, we compute the supersymmetric Casimir energy of each flat Wilson line for $SU(2)$ SYM on different ADE singularities and find that the flat Wilson lines all have the same supersymmetric Casimir energy. We argue that this exact degeneracy is peculiar to $\mathcal{N}=4$ supersymmetry and show that the degeneracy is lifted when the number of supersymmetry is reduced. In particular, we uncover a surprising result for the ground state structure of the conformal $\mathcal{N}=2$ $SU(2)$ four-flavor theory on $S^3/\Gamma$. For $\mathcal{N}=4$ SYM, S-duality maps the ground state Wilson lines to ground state t' Hooft lines taking values in the Langlands dual group. We show that the supersymmetric Casimir energy of the t' Hooft line ground states is the same as the Wilson line ground states. This can be viewed as a ground state test of S-duality.
In High Energy Physics, detailed and time-consuming simulations are used for particle interactions with detectors. To bypass these simulations with a generative model, the generation of large point clouds in a short time is required, while the complex dependencies between the particles must be correctly modelled. Particle showers are inherently tree-based processes, as each particle is produced by the decay or detector interaction of a particle of the previous generation. In this work, we present a significant extension to DeepTreeGAN, featuring a critic, that is able to aggregate such point clouds iteratively in a tree-based manner. We show that this model can reproduce complex distributions, and we evaluate its performance on the public JetNet 150 dataset.
This paper presents the measurement of fiducial and differential cross sections for both the inclusive and electroweak production of a same-sign $W$-boson pair in association with two jets ($W^\pm W^\pm jj$) using 139 fb$^{-1}$ of proton-proton collision data recorded at a centre-of-mass energy of $\sqrt{s}=13$ TeV by the ATLAS detector at the Large Hadron Collider. The analysis is performed by selecting two same-charge leptons, electron or muon, and at least two jets with large invariant mass and a large rapidity difference. The measured fiducial cross sections for electroweak and inclusive $W^\pm W^\pm jj$ production are $2.92 \pm 0.22\, \text{(stat.)} \pm 0.19\, \text{(syst.)}$ fb and $3.38 \pm 0.22\, \text{(stat.)} \pm 0.19\, \text{(syst.)}$ fb, respectively, in agreement with Standard Model predictions. The measurements are used to constrain anomalous quartic gauge couplings by extracting 95% confidence level intervals on dimension-8 operators. A search for doubly charged Higgs bosons $H^{\pm\pm}$ that are produced in vector-boson fusion processes and decay into a same-sign $W$ boson pair is performed. The largest deviation from the Standard Model occurs for an $H^{\pm\pm}$ mass near 450 GeV, with a global significance of 2.5 standard deviations.
In particle physics, Monte Carlo (MC) event generators are needed to compare theory to the measured data. Many MC samples have to be generated to account for theoretical systematic uncertainties, at a significant computational cost. Therefore, the MC statistic becomes a limiting factor for most measurements and the significant computational cost of these programs a bottleneck in most physics analyses. In this contribution, the Deep neural network using Classification for Tuning and Reweighting (DCTR) approach is evaluated for the reweighting of two systematic uncertainties in MC simulations of top quark pair production within the CMS experiment. DCTR is a method, based on a Deep Neural Network (DNN) technique, to reweight simulations to different model parameters by using the full kinematic information in the event. This methodology avoids the need for simulating the detector response multiple times by incorporating the relevant variations in a single sample.
The HL-LHC run is anticipated to start at the end of this decade and will pose a significant challenge for the scale of the HEP software and computing infrastructure. The mission of the U.S. CMS Software & Computing Operations Program is to develop and operate the software and computing resources necessary to process CMS data expeditiously and to enable U.S. physicists to fully participate in the physics of CMS. We have developed a strategic plan to prioritize R&D efforts to reach this goal for the HL-LHC. This plan includes four grand challenges: modernizing physics software and improving algorithms, building infrastructure for exabyte-scale datasets, transforming the scientific data analysis process and transitioning from R&D to operations. We are involved in a variety of R&D projects that fall within these grand challenges. In this talk, we will introduce our four grand challenges and outline the R&D program of the U.S. CMS Software & Computing Operations Program.
A search for a long-lived, heavy neutral lepton ($\mathcal{N}$) in 139 fb$^{-1}$ of $\sqrt{s}=13$ TeV $pp$ collision data collected by the ATLAS detector at the Large Hadron Collider is reported. The $\mathcal{N}$ is produced via $W \rightarrow \mathcal{N} \mu$ or $W \rightarrow \mathcal{N} e$ and decays into two charged leptons and a neutrino, forming a displaced vertex. The $\mathcal{N}$ mass is used to discriminate between signal and background. No signal is observed, and limits are set on the squared mixing parameters of the $\mathcal{N}$ with the left-handed neutrino states for the $\mathcal{N}$ mass range $3$ GeV $< m_{\mathcal{N}} < 15$ GeV. For the first time, limits are given for both single-flavor and multiflavor mixing scenarios motivated by neutrino flavor oscillation results for both the normal and inverted neutrino-mass hierarchies.
The ATLAS detector at the Large Hadron Collider has been used to measure jet substructure modification and suppression in Pb+Pb collisions at a nucleon-nucleon center-of-mass energy $\sqrt{s_{_\mathrm{NN}}}=5.02~\mathrm{TeV}$ in comparison with $pp$ collisions at $\sqrt{s}=5.02~\mathrm{TeV}$. The Pb+Pb data, collected in 2018, have an integrated luminosity of $1.72~\mathrm{nb^{-1}}$, while the $pp$ data, collected in 2017, have an integrated luminosity of $260~\mathrm{pb}^{-1}$. Jets used in this analysis are clustered using the anti-$k_{t}$ algorithm with a radius parameter $R=0.4$. The jet constituents, defined by both tracking and calorimeter information, are used to determine the angular scale $r_\mathrm{g}$ of the first hard splitting inside the jet by reclustering them using the Cambridge-Aachen algorithm and employing the soft-drop grooming technique. The nuclear modification factor, $R_\mathrm{AA}$, used to characterize jet suppression in Pb+Pb collisions, is presented differentially in $r_\mathrm{g}$, jet transverse momentum, and in intervals of collision centrality. The $R_\mathrm{AA}$ value is observed to depend significantly on jet $r_\mathrm{g}$. Jets produced with the largest measured $r_\mathrm{g}$ are found to be twice as suppressed as those with the smallest $r_\mathrm{g}$ in central Pb+Pb collisions. The $R_\mathrm{AA}$ values do not exhibit a strong variation with jet $p_\mathrm{T}$ in any of the $r_\mathrm{g}$ intervals. The $r_\mathrm{g}$ and $p_\mathrm{T}$ dependence of jet $R_\mathrm{AA}$ is qualitatively consistent with a picture of jet quenching arising from coherence and provides the most direct evidence in support of this approach.
A search for pair-produced scalar or vector leptoquarks decaying into a $b$-quark and a $\tau$-lepton is presented using the full LHC Run 2 (2015-2018) data sample of 139 fb$^{-1}$ collected with the ATLAS detector in proton-proton collisions at a centre-of-mass energy of $\sqrt{s} =13$ TeV. Events in which at least one $\tau$-lepton decays hadronically are considered, and multivariate discriminants are used to extract the signals. No significant deviations from the Standard Model expectation are observed and 95% confidence-level upper limits on the production cross-section are derived as a function of leptoquark mass and branching ratio $B$ into a $\tau$-lepton and $b$-quark. For scalar leptoquarks, masses below 1460 GeV are excluded assuming $B=100$%, while for vector leptoquarks the corresponding limit is 1650 GeV (1910 GeV) in the minimal-coupling (Yang-Mills) scenario.
Using $e^+e^-$ collision data corresponding to a total integrated luminosity of 12.9 $fb^{-1}$ collected with the BESIII detector at the BEPCII collider, the exclusive Born cross sections and the effective form factors of the reaction $e^+e^-\rightarrow\Xi^{-}\bar\Xi^{+}$ are measured via the single baryon-tag method at 23 center-of-mass energies between 3.510 and 4.843 GeV. Evidence for the decay $\psi(3770)\rightarrow\Xi^{-}\bar\Xi^{+}$ is observed with a significance of 4.5$\sigma$ by analyzing the measured cross sections together with earlier BESIII results. For the other charmonium(-like) states $\psi(4040)$, $\psi(4160)$, $Y(4230)$, $Y(4360)$, $\psi(4415)$, and $Y(4660)$, no significant signal of their decay to $\Xi^-\bar \Xi^+$ is found. For these states, upper limits of the products of the branching fraction and the electronic partial width at the 90% confidence level are provided.
We compare the well known Rayleigh-Ritz variational method (RRVM) with a recently proposed approach based on supersymmetric quantum mechanics and the Gram-Schmidt orthogonalization method (SSQMGS). We apply both procedures to a particular class of double-well harmonic oscillators that had been conveniently chosen for the application of the latter approach. The RRVM eigenvalues converge smoothly from above providing much more accurate results with less computational effort. Present results show that the unproved SSQMGS upper bounds do not hold.
Ion trap systems built upon microfabricated chips have emerged as a promising platform for quantum computing to achieve reproducible and scalable structures. However, photo-induced charging of materials in such chips can generate undesired stray electric fields that disrupt the quantum state of the ion, limiting high-fidelity quantum control essential for practical quantum computing. While crude understanding of the phenomena has been gained heuristically over the past years, explanations for the microscopic mechanism of photo-generated charge carrier dynamics remains largely elusive. Here, we present a photo-induced charging model for semiconductors, whose verification is enabled by a systematic interaction between trapped ions and photo-induced stray fields from exposed silicon surfaces in our chip. We use motion-sensitive qubit transitions to directly characterize the stray field and analyze its effect on the quantum dynamics of the trapped ion. In contrast to incoherent errors arising from the thermal motion of the ion, coherent errors are induced by the stray field, whose effect is significantly imprinted during the quantum control of the ion. These errors are investigated in depth and methods to mitigate them are discussed. Finally, we extend the implications of our study to other photo-induced charging mechanisms prevalent in ion traps.
The quantum approximate optimization algorithm (QAOA) shows great computational potential on combinatorial optimization problems. It is a promising algorithm on near-term quantum devices, but one of the difficulty in application of QAOA is the complexity of parameter setting. In this paper, an adiabatic-passage based parameter setting method is proposed and applied to 3-SAT. And in simulation, the optimization cost is significantly reduced, approximately between sublinear to logarithmic on the depth $p$ of QAOA. The efficiency of this method mainly stems from two aspects, one is the problem-oriented preprocessing of Hamiltonian, and the other is the parameter space adjustment based on the continuity of adiabatic passage. Firstly, a random model for 3-SAT is provided and the problem Hamiltonian of this model is designed as a random matrix. Based on the statistical property of randomized Hamiltonian, the Hamiltonian of QAOA is preprocessed and the parameter setting is seperated from the overal property of the problem. As a result, a good initialization can be obtained. Secondly, the optimal adiabatic passage is introduced and actrually, the QAOA can be regarded as the parameterization of adiabatic passage and the optimization as the search of the optimal adiabatic passage. Based on this, the adiabatic passage is parameterized as another parameter space with better continuity and the adiabatic-passage based parameter setting method is proposed.
The classical simulation of highly-entangling quantum dynamics is conjectured to be generically hard. Thus, recently discovered measurement-induced transitions between highly entangling and low-entanglement dynamics are phase transitions in classical simulability. Here, we study simulability transitions beyond entanglement: noting that some highly-entangling dynamics (e.g., integrable systems or Clifford circuits) are easy to classically simulate, thus requiring "magic"--a subtle form of quantum resource--to achieve computational hardness, we ask how the dynamics of magic competes with measurements. We study the resulting "dynamical magic transitions" focusing on random monitored Clifford circuits doped by T gates (injecting magic). We identify dynamical "stabilizer-purification"--the collapse of a superposition of stabilizer states by measurements--as the mechanism driving this transition. We find cases where transitions in magic and entanglement coincide, but also others with a magic and simulability transition in a highly (volume-law) entangled phase. In establishing our results, we use Pauli-based computation, a scheme distilling the quantum essence of the dynamics to a magic state register subject to mutually commuting measurements. We link stabilizer-purification to "magic fragmentation" wherein these measurements separate into disjoint, O(1)-weight blocks, and relate this to the spread of magic in the original circuit becoming arrested.
Quantum Markov models are employed ubiquitously in quantum physics and in quantum information theory due to their relative simplicity and analytical tractability. In particular, these models are known to give accurate approximations for a wide range of quantum optical and mesoscopic systems. However, in general, the validity of the Markov approximation entails assumptions regarding properties of the system of interest and its environment, which may not be satisfied or accurate in arbitrary physical systems. Therefore, developing useful modelling tools for general non-Markovian quantum systems for which the Markov approximation is inappropriate or deficient is an undertaking of significant importance. This work considers non-Markovian principal quantum systems that can be embedded in a larger Markovian quantum system with one or more compound baths consisting of an auxiliary quantum system and a quantum white noise field, and derives a set of coupled stochastic and quantum master equations for embedded non-Markovian quantum systems. The case of a purely Hamiltonian coupling between the principal and auxiliary systems as a closed system without coupling to white noises is included as a special case. The results are expected to be of interest for (open-loop and feedback) control of continuous-time non-Markovian systems and studying reduced models for numerical simulation of such systems. They may also shed more light on the general structure of continuous-time non-Markovian quantum systems.
A chiral coordinate Bethe ansatz method is developed to study the periodic XYZ chain. We construct a set of chiral vectors with fixed number of kinks. All vectors are factorized and have simple structures. Under roots of unity conditions, the Hilbert space has an invariant subspace and our vectors form a basis of this subspace. We propose a Bethe ansatz solely based on the action of the Hamiltonian on the chiral vectors, avoiding the use of transfer matrix techniques. This allows to parameterize the expansion coefficients and derive the homogeneous Bethe ansatz equations whose solutions give the exact energies and eigenstates. Our analytic results agree with earlier approaches, notably by Baxter, and are supported by numerical calculations.
We describe a quantum error-detecting and error-correcting code embedded within irreducible representations of SU(3). These logical qutrits inherit the He(3) symmetries induced by the representation, while protecting against small SU(3) displacements. We explore the general methodology for finding codes from structure-inducing representations of groups, together with symmetries inherited from finite subgroups, extending the case of spin representations of SU(2).
We analyse recent contradictory results and conclusions about the spectrum of the screened Coulomb potential $V(r)=-r^{-1}e^{-C/r}$. The well known Hellmann-Feynman theorem shows that all the bound states of the Coulomb potential ($C=0$) remain bounded as $C$ increases. We derive a simple approximate analytical expression for the eigenvalues for sufficiently small values of the screening parameter $C$ and an approximate asymptotic expression for the asymptotic behaviour of the s-state eigenvalues when $C\rightarrow \infty $. Present results are expected to resolve the discrepancy about the spectrum of the quantum-mechanical model just mentioned.
Entangled two-photon absorption (ETPA) could form the basis of nonlinear quantum spectroscopy at very low photon fluxes, since, at sufficiently low photon fluxes, ETPA scales linearly with the photon flux. When different pairs start to overlap temporally, accidental coincidences are thought to give rise to a 'classical' quadratic scaling which dominates the signal at large photon fluxes and thus recovers a supposedly classical regime, where any quantum advantage is thought to be lost. Here we scrutinize this assumption and demonstrate that quantum-enhanced absorption cross sections can persist even to very large photon numbers. To this end, we use a minimal model for quantum light, which can interpolate continuously between the entangled pair and a high-photon-flux limit, to derive analytically ETPA cross sections and the intensity crossover regime. We investigate the interplay between spectral and spatial degrees of freedom, how linewidth broadening of the sample impacts the experimentally achievable enhancement.
Quantum subspace methods (QSMs) are a class of quantum computing algorithms where the time-independent Schrodinger equation for a quantum system is projected onto a subspace of the underlying Hilbert space. This projection transforms the Schrodinger equation into an eigenvalue problem determined by measurements carried out on a quantum device. The eigenvalue problem is then solved on a classical computer, yielding approximations to ground- and excited-state energies and wavefunctions. QSMs are examples of hybrid quantum-classical methods, where a quantum device supported by classical computational resources is employed to tackle a problem. QSMs are rapidly gaining traction as a strategy to simulate electronic wavefunctions on quantum computers, and thus their design, development, and application is a key research field at the interface between quantum computation and electronic structure. In this review, we provide a self-contained introduction to QSMs, with emphasis on their application to the electronic structure of molecules. We present the theoretical foundations and applications of QSMs, and we discuss their implementation on quantum hardware, illustrating the impact of noise on their performance.
The coherent quantum Zeno dynamics (QZD) is a special unitary time evolution in which a quantum population transition gets constrained in a subspace of the entire Hilbert space. We show that coherent QZD can be categorized by orders for the first time, where only the zeroth-order type has been well investigated. In this paper, we focus on the little-known first-order coherent QZD (FC-QZD). We also construct some chain-like systems described by the tight-binding model which establishes FC-QZD in the form of a surprisingly nonlocal end-to-end population transition.
In this paper we consider the scalability of Multi-Angle QAOA with respect to the number of QAOA layers. We found that MA-QAOA is able to significantly reduce the depth of QAOA circuits, by a factor of up to 4 for the considered data sets. However, MA-QAOA is not optimal for minimization of the total QPU time. Different optimization initialization strategies are considered and compared for both QAOA and MA-QAOA. Among them, a new initialization strategy is suggested for MA-QAOA that is able to consistently and significantly outperform random initialization used in the previous studies.
Titanium:Sapphire (Ti:Sa) lasers have been essential for advancing fundamental research and technological applications. Ti:Sa lasers are unmatched in bandwidth and tuning range, yet their use is severely restricted due to their large size, cost, and need for high optical pump powers. Here, we demonstrate a monocrystalline Ti:Sa-on-insulator (Ti:SaOI) photonics platform which enables dramatic miniaturization, cost-reduction, and scalability of Ti:Sa technology. First, through fabrication of low-loss whispering gallery mode resonators, we realize a Ti:Sa laser operating with an ultra-low lasing threshold of 290 $\mu$W. Then, through orders-of-magnitude improvement in mode confinement in Ti:SaOI waveguides, we realize the first integrated solid-state (i.e., non-semiconductor) optical amplifier operating below 1 $\mu$m, with an ultra-wide bandwidth of 700 - 950 nm and peak gain of 64 dB/cm. We demonstrate unprecedented 17 dB distortion-free amplification of picosecond pulses to up to 2.3 nJ pulse energy, corresponding to a peak power of 1.0 kW. Finally, we demonstrate the first tunable integrated Ti:Sa laser, featuring narrow linewidths and a 24.7 THz tuning range, which, for the first time, can be pumped with low-cost, miniature, off-the-shelf green laser diodes. This opens doors to new modalities of Ti:Sa lasers (now occupying a footprint less than 0.15 mm$^2$), such as massively-scalable Ti:Sa laser array systems for a variety of applications. As a proof-of-concept demonstration, we employ a Ti:SaOI laser array as the sole optical control for a cavity quantum electrodynamics experiment with artificial atoms in silicon carbide. This work is a key step towards the democratization of Ti:Sa technology through a three orders-of-magnitude reduction in cost and footprint, as well as the introduction of solid-state broadband amplification of sub-micron wavelength light.
Financial services is a prospect industry where unlocked near-term quantum utility could yield profitable potential, and, in particular, quantum machine learning algorithms could potentially benefit businesses by improving the quality of predictive models. Quantum kernel methods have demonstrated success in financial, binary classification tasks, like fraud detection, and avoid issues found in variational quantum machine learning approaches. However, choosing a suitable quantum kernel for a classical dataset remains a challenge. We propose a hybrid, quantum multiple kernel learning (QMKL) methodology that can improve classification quality over a single kernel approach. We test the robustness of QMKL on several financially relevant datasets using both fidelity and projected quantum kernel approaches. We further demonstrate QMKL on quantum hardware using an error mitigation pipeline and show the benefits of QMKL in the large qubit regime.
Quantum Annealers (QAs) operate as single-instruction machines, lacking a SWAP operation to overcome limited qubit connectivity. Consequently, multiple physical qubits are chained to form a program qubit with higher connectivity, resulting in a drastically diminished effective QA capacity by up to 33x. We observe that in QAs: (a) chain lengths exhibit a power-law distribution, a few dominant chains holding substantially more qubits than others; and (b) about 25% of physical qubits remain unused, getting isolated between these chains. We propose Skipper, a software technique that enhances the capacity and fidelity of QAs by skipping dominant chains and substituting their program qubit with two readout results. Using a 5761-qubit QA, we demonstrate that Skipper can tackle up to 59% (Avg. 28%) larger problems when eleven chains are skipped. Additionally, Skipper can improve QA fidelity by up to 44% (Avg. 33%) when cutting five chains (32 runs). Users can specify up to eleven chain cuts in Skipper, necessitating about 2,000 distinct quantum executable runs. To mitigate this, we introduce Skipper-G, a greedy scheme that skips sub-problems less likely to hold the global optimum, executing a maximum of 23 quantum executables with eleven chain trims. Skipper-G can boost QA fidelity by up to 41% (Avg. 29%) when cutting five chains (11 runs).
BosonSampling is the leading candidate for demonstrating quantum computational advantage in photonic systems. While we have recently seen many impressive experimental demonstrations, there is still a formidable distance between the complexity-theoretic hardness arguments and current experiments. One of the largest gaps involves the ratio of photons to modes: all current hardness evidence assumes a "high-mode" regime in which the number of linear optical modes scales at least quadratically in the number of photons. By contrast, current experiments operate in a "low-mode" regime with a linear number of modes. In this paper we bridge this gap, bringing the hardness evidence for the low-mode experiments to the same level as had been previously established for the high-mode regime. This involves proving a new worst-to-average-case reduction for computing the Permanent that is robust to large numbers of row repetitions and also to distributions over matrices with correlated entries.
Data visualization is important in understanding the characteristics of data that are difficult to see directly. It is used to visualize loss landscapes and optimization trajectories to analyze optimization performance. Popular optimization analysis is performed by visualizing a loss landscape around the reached local or global minimum using principal component analysis. However, this visualization depends on the variational parameters of a quantum circuit rather than quantum states, which makes it difficult to understand the mechanism of optimization process through the property of quantum states. Here, we propose a quantum data visualization method using quantum kernels, which enables us to offer fast and highly accurate visualization of quantum states. In our numerical experiments, we visualize hand-written digits dataset and apply $k$-nearest neighbor algorithm to the low-dimensional data to quantitatively evaluate our proposed method compared with a classical kernel method. As a result, our proposed method achieves comparable accuracy to the state-of-the-art classical kernel method, meaning that the proposed visualization method based on quantum machine learning does not degrade the separability of the input higher dimensional data. Furthermore, we visualize the optimization trajectories of finding the ground states of transverse field Ising model and successfully find the trajectory characteristics. Since quantum states are higher dimensional objects that can only be seen via observables, our visualization method, which inherits the similarity of quantum data, would be useful in understanding the behavior of quantum circuits and algorithms.
In recent years, Classical Convolutional Neural Networks (CNNs) have been applied for image recognition successfully. Quantum Convolutional Neural Networks (QCNNs) are proposed as a novel generalization to CNNs by using quantum mechanisms. The quantum mechanisms lead to an efficient training process in QCNNs by reducing the size of input from $N$ to $log_2N$. This paper implements and compares both CNNs and QCNNs by testing losses and prediction accuracy on three commonly used datasets. The datasets include the MNIST hand-written digits, Fashion MNIST and cat/dog face images. Additionally, data augmentation (DA), a technique commonly used in CNNs to improve the performance of classification by generating similar images based on original inputs, is also implemented in QCNNs. Surprisingly, the results showed that data augmentation didn't improve QCNNs performance. The reasons and logic behind this result are discussed, hoping to expand our understanding of Quantum machine learning theory.
Quantum systems are unavoidably open to their surrounding degrees of freedom. The theory of open quantum systems is thus crucial to understanding the fluctuations, dissipation, and decoherence of a quantum system of interest. Typically, the bath is modeled as an ensemble of harmonic oscillators, which yields Gaussian statistics of the bath influence on the quantum systems. However, there are also phenomena in which the bath consists of two-state systems, spins, or anharmonic oscillators; therefore, the non-Gaussian properties of the bath become important. Nevertheless, a theoretical framework to describe quantum systems under the influence of such non-Gaussian baths is not well established. Here, we develop a theory describing quantum dissipative systems affected by Poisson noise properties of the bath as Poisson noise is fundamental in describing non-Gaussian white noises. In contrast to past studies that modeled the bath as a classical stochastic noise source producing only pure dephasing, we introduce a quantum bath model that allows for the consistent description of dissipative quantum systems. The property of the constructed bath model is consistent with the Poisson noise properties when the bath correlation time is short and the bath interacts with the quantum system strongly but discretely. The obtained results reveal non-Gaussian bath effects in the white noise regime, and they provide an essential step toward describing open quantum dynamics under the influence of generic non-Gaussian baths. Our findings can be used to design baths with non-Gaussian properties for dissipative quantum state engineering in quantum information science, as well as to explore non-Gaussian bath effects in biophysical chemistry and condensed matter physics.
The eigenstate thermalization hypothesis for translation invariant quantum spin systems has been proved recently by using random matrices. In this paper, we study the stronger subsystem version of eigenstate thermalization hypothesis for translation invariant quantum systems without referring to random matrices. By showing the small upper bounds on the quantum variance or the Belavkin-Staszewski relative entropy, we prove the subsystem eigenstate thermalization hypothesis for translation invariant quantum systems in an elementary way.
We report the development of a continuous-wave and pulsed X-band electron spin resonance (ESR) spectrometer for the study of spins on ordered surfaces down to cryogenic temperatures. The spectrometer operates in ultra-high vacuum and utilizes a half-wavelength microstrip line resonator realized using epitaxially grown copper films on single crystal Al$_2$O$_3$ substrates. The one-dimensional microstrip line resonator exhibits a quality factor of more than 200 at room temperature, close to the upper limit determined by radiation losses. The surface characterizations of the copper strip of the resonator by atomic force microscope, low-energy electron diffraction, and scanning tunneling microscope show that the surface is atomically clean, flat, and single crystalline. Measuring the ESR spectrum at 15 K from a few nm thick molecular film of YPc$_2$, we find a continuous-wave ESR sensitivity of $6.5 \cdot 10^{10}~\text{spins}/\text{G} \cdot \text{Hz}^{1/2}$ indicating that a signal-to-noise ratio of $7.7~\text{G} \cdot \text{Hz}^{1/2}$ is expected from a monolayer of YPc$_2$ molecules. Advanced pulsed ESR experimental capabilities including dynamical decoupling and electron-nuclear double resonance are demonstrated using free radicals diluted in a glassy matrix.
We develop and implement an exact conical intersection nonadiabatic wave packet dynamics method that combines the local diabatic representation, Strang splitting for the total molecular propagator, and discrete variable representation with uniform grids. By employing the local diabatic representation, this method captures all non-adiabatic effects, including nonadiabatic transitions, electronic coherences, and geometric phases. Moreover, it is free of singularities in the first and second derivative couplings, and does not require a smooth gauge of electronic wavefunction phase. We further show that in contrast to the adiabatic representation, the split-operator method can be directly applied to the full molecular propagator with the locally diabatic ansatz. The Fourier series, employed as the primitive nuclear basis functions, is universal and can be applied to all types of reactive coordinates. The combination of local diabatic representation, Strang splitting, and Fourier basis allows exact modeling of conical intersection quantum dynamics directly with adiabatic electronic states that can be obtained from standard electronic structure computations.
The persistent separability of certain quantum states, known as symmetric absolutely separable (SAS), under symmetry-preserving global unitary transformations is of key significance in the context of quantum resources for bosonic systems. In this work, we develop criteria for detecting SAS states of any number of qubits. Our approach is based on the Glauber-Sudarshan $P$ representation for finite-dimensional quantum systems. We introduce three families of SAS witnesses, one linear and two non-linear in the state $\rho$, formulated respectively as an algebraic inequality or a quadratic optimization problem. These witnesses are capable of identifying more SAS states than previously known counterparts. We also explore the geometric properties of the subsets of SAS states detected by our witnesses, shedding light on their distinctions.
Controllable single-photon routing plays an important role in quantum networks. We investigate single-photon scattering in two one-dimensional (1D) waveguides by a three-level emitter with a cascade configuration, which is a dipole coupled to an extra cavity. The tunneling path for the transmission of a single photon is switched by whether the extra cavity contains photons. For the setup, the Autler-Townes splitting is modulated by the extra cavity, in which the transmission valley (reflection range) width is tunable in terms of the number of photons in the extra cavity. Our investigation will be beneficial to single-photon routing in quantum networks using quantifiable photon numbers in an extra cavity.
We present two families of quantum speed limits (QSLs) for finite-dimensional quantum systems undergoing a general physical process. These QSLs were obtained using Schatten $\alpha$-norms, firstly exploiting the geometric features of the space of quantum states endowed with some inner product, and secondly employing the Holder's inequality for matrix norms. In particular, for the case of single-qubit states, we find that the geometric QSL is independent of the Schatten norm chosen, thus revealing a universality behavior of such quantifiers. Furthermore, we provide a comparison of these quantum speed limits with existing paradigmatic QSLs in literature, thus showing that the latter results represent particular cases of a general class of QSLs related to Schatten $\alpha$-norms. Noteworthy, we address necessary and sufficient conditions for the tightness of the quantum speed limit that mostly depends on the populations and quantum coherences of the evolved single-qubit state, and also present a geometric interpretation for these set of conditions. Finally, we compare the two QSL obtained for the dynamics of single-qubit states, also presenting an inequality between them that has a clear geometrical meaning.
Improving the measurement precision of low temperature is significant in fundamental science and advanced quantum technology application. However, the measurement precision of temperature $T$ usually diverges as $T$ tends to 0. Here, by utilizing a heavy impurity to measure the temperature of a Bose gas, we obtain the Landau bound to precision $\delta^2 T\propto T^2$ to avoid the divergence. Moreover, when the initial momentum of the heavy impurity is fixed and non-zero, the measurement precision can be $\delta^2 T\propto T^3$ to break the Landau bound. We derive the momentum distribution of the heavy impurity at any moment and obtain the optimal measurement precision of the temperature by calculating the Fisher information. As a result, we find that enhancing the expectation value of the initial momentum can help to improve the measurement precision. In addition, the momentum measurement is the optimal measurement of the temperature in the case of that the initial momentum is fixed and not equal to 0. The kinetic energy measurement is the optimal measurement in the case of that the expectation value of the initial momentum is 0. Finally, we obtain that the temperatures of two Bose gases can be measured simultaneously. The simultaneous measurement precision is proportional to $T^2$ when two temperatures are close to $T$.
Cold ions in traps are well-established, highly controllable quantum systems with a wide variety of applications in quantum information, precision spectroscopy, clocks and chemistry. Nanomechanical oscillators are used in advanced sensing applications and for exploring the border between classical and quantum physics. Here, we report on the implementation of a hybrid system combining a metallic nanowire with laser-cooled ions in a miniaturised ion trap. We demonstrate resonant and off-resonant coupling of the two systems and the coherent motional excitation of the ion by the mechanical drive of the nanowire. The present results open up avenues for mechanically manipulating the quantum motion of trapped ions, for the development of ion-mechanical hybrid quantum systems and for the sympathetic cooling of mechanical systems by trapped ions and vice versa.
In this paper, we generalize to three players the well-known CHSH quantum game. To do so, we consider all possible 3 variables Boolean functions and search among them which ones correspond to a game scenario with a quantum advantage (for a given entangled state). In particular we provide two new three players quantum games where, in one case, the best quantum strategy is obtained when the players share a $GHZ$ state, while in the other one the players have a better advantage when they use a $W$ state as their quantum resource. To illustrate our findings we implement our game scenarios on an online quantum computer and prove experimentally the advantage of the corresponding quantum resource for each game.
$^{171}$Yb$^{3+}$-doped Y$_2$SiO$_5$ crystals are a promising platform for optical quantum memories in long-distance quantum communications. The relevance of this material lies in $^{171}$Yb long optical and spin coherence times, along with a large hyperfine splitting, enabling long quantum storage over large bandwidths. Mechanisms affecting the optical decoherence are however not precisely known, especially since low-temperature measurements have so far focused on the 2 to 4 K range. In this work, we performed two- and three-pulse photon echoes and spectral hole burning to determine optical homogeneous linewidths in two 171 Yb:YSO crystals doped at 2 and 10 ppm. Experiments were performed in the 40 mK to 18 K temperature range, leading to linewidths between 320 Hz, among the narrowest reported for rare-earth ions, and several MHz. Our results show that above 6 K the homogeneous linewidth is mainly due to an elastic two-phonon process which results in a slow broadening with temperature, the homogeneous linewidth reaching only 25 kHz at 10 K. At lower temperatures, interactions with $^{89}$Yb nuclear spin-flips, paramagnetic defects or impurities, and also Yb-Yb interactions for the higher concentrated crystal, are likely the main limiting factor to the homogeneous linewidth. In particular, we conclude that the direct effect of spin and optical excited state lifetime is a minor contribution to optical decoherence in the whole temperature range studied. Our results indicate possible paths and regimes for further decreasing the homogeneous linewidths or maintaining narrow lines at higher $^{171}$Yb concentration.
Two theorems with applications to the quantum theory of measurements are stated and proven. The first one clarifies and amends von Neumann's Measurement Postulate used in the Copenhagen interpretation of quantum mechanics. The second one clarifies the relationship between ``events'' and ``measurements'' and the meaning of measurements in the $ETH$-Approach to quantum mechanics.
Our cavity quantum electrodynamics calculations demonstrate generation of steady-state entanglement between a plasmonically coupled pair of quantum dots by using single-mode squeezed light source. We show that strong coupling of plasmons to the incoming light source and the pairwise nature of squeezed photon generation enable the formation of entanglement between the initially unexcited quantum dots. The entanglement of quantum dots, measured as concurrence, can be improved replacing a pulsed source of light to continuous pumping of squeezed photons. Unlike previously introduced schemes the concurrence is robust against variations in the system parameters. Specifically, the generation of entanglement does not rely on fine tuning of plasmon quantum dot coupling. This work provides a new perspective for robust entangled state preparation in open quantum systems.
The Wehrl entropy of a quantum state is the entropy of the coherent-state distribution function (Husimi function), and is non-zero even for pure states. We investigate the Wehrl entropy for $N$ spin-1/2 particles with respect to SU(2)$^{\otimes N}$ coherent states (i.e., the direct products of spin coherent states of each particle). We focus on: (1) The statistical interpretation of this Wehrl entropy. (2) The relationship between the Wehrl entropy and quantum entanglement. For (1), despite the coherent states not forming a group of orthonormal bases, we prove that the Wehrl entropy can still be interpreted as the entropy of a probability distribution with clear physical meaning. For (2), we numerically calculate the Wehrl entropy of various entangled pure states with particle number $2\leq N\leq 20$. Our results show that for the large-$N$ ($N\gtrsim 10$) systems the Wehrl entropy of the highly chaotic entangled states are much larger than that of the regular ones (e.g., the GHZ state). These results, together with the fact that the Wehrl entropy is invariant under local unitary transformations, indicate that the Wehrl entropy can reflect the complexity of the quantum entanglement (entanglement complexity) of many-body pure states, as A. Sugita proposed directly from the definitions of the Husimi function and Wehrl entropy (Jour. Phys. A 36, 9081 (2003)). Furthermore, the Wehrl entropy per particle can serve as a quantitative description of this complexity. We further show that the many-body pure entangled states can be classified into three types, according to the behaviors of the Wehrl entropy per particle in the limit $N\rightarrow\infty$, with the states of each type having very different entanglement complexity.
We perform a comparative study of the phase noise induced in the lasers used for Bragg diffraction in a Bose-Einstein condensate-based quantum gravimeter where the Bragg beams are generated using two different configurations. In one of the configurations, the Bragg beams that form the moving optical lattice are generated using two different acousto-optic modulators. In the second configuration, the Bragg beams are generated using a single acousto-optic modulator carrying two phase-locked frequencies. The second configuration shows a suppression of phase noise by a factor of 4.7 times in the frequency band upto 10 $kHz$, the primary source of noise, which is the background acoustic noise picked up by optical components and the optical table. We report a sensitivity of 99.7 $\mu Gal/\sqrt Hz$ for an interferometric time of 10 $ms$.
Here we show that the concepts behind such terms as entanglement, qubits, quantum gates, quantum error corrections, unitary time evolution etc., which are usually ascribed to quantum systems, can be adequately realized on a set of coupled classical pendulums.
In a conducting medium held at finite temperature, free carriers are performing Brownian motion and generate fluctuating electromagnetic fields. We compute the averaged Lorentz force density that turns out nonzero in a thin sub-surface layer, pointing towards the surface, while vanishing in the bulk. This is an elementary example of rectified fluctuations, similar to the Casimir force or radiative heat transport. Our results also provide an experimental way to distinguish between the Drude and so-called plasma models.
In radio-frequency trap, the temperature of ion ensembles converges towards a hot equilibrium due to radio-frequency heating. This effect is detrimental to the stability of trapped ensembles and is the justification of cooling. The intensity of this radio-frequency heating increases with the amplitude of the radio-frequency field $q_x$. Using an analytical empirical formula, we show that the lifetime of the ion ensemble $t_0$ under cold condition increases with $q_x$ according to a power law $t_0\propto q_x^A$, and does not vary significantly for the several ion quantities $N$ tested. The temperature of the explosive onset $B$ decreases linearly with $q_x$. We also show that non-linear instabilities due to trapping parameters decreases $t_0$ and $B$, and produce a local increase of heating rate for certain temperature ranges.
We introduce a novel approach for estimating the spectrum of quantum many-body Hamiltonians, and more generally, of Hermitian operators, using quantum time evolution. In our approach we are evolving a maximally mixed state under the Hamiltonian of interest and collecting specific time-series measurements to estimate its spectrum. We demonstrate the advantage of our technique over currently used classical statistical sampling methods. We showcase our approach by experimentally estimating the spectral decomposition of a 2-qubit Heisenberg Hamiltonian on an IBM Quantum backend. For this purpose, we develop a hardware-efficient decomposition that controls $n$-qubit Pauli rotations against the physically closest qubit alongside expressing two-qubit rotations in terms of the native entangling interaction. This substantially reduced the accumulation of errors from noisy two-qubit operations in time evolution simulation protocols. We conclude by discussing the potential impact of our work and the future directions of research it opens.
The Linear Combination of Unitaries (LCU) method is a powerful scheme for the block encoding of operators but suffers from high overheads. In this work, we discuss the parallelisation of LCU and in particular the SELECT subroutine of LCU based on partitioning of observables into groups of commuting operators, as well as the use of adaptive circuits and teleportation that allow us to perform required Clifford circuits in constant depth. We only require an $O(\log n)$ factor increase in the number of qubits in order to produce a significant depth reduction, with evidence suggesting that for practical molecular Hamiltonians, the depth saving is $O(n)$, and calculate a depth saving of $20\times$ for SELECT on a H$_2$O Hamiltonian, even though small problem sizes are the worst case for our scheme. We discuss the implications of our method in the fault-tolerant setting, noting that parallelisation reduces the $T$-depth by the same factor as the logical algorithm, without changing the $T$-count.
The quantum singular value transformation has revolutionised quantum algorithms. By applying a polynomial to an arbitrary matrix, it provides a unifying picture of quantum algorithms. However, polynomials are restricted to definite parity and real coefficients, and finding the circuit (the phase factors) has proven difficult in practice. Recent work has removed these restrictions and enabled faster computation of phase factors, yet only for unitary matrices. Here we propose two generalisations. The generalised quantum singular value transformation allows complex polynomials for arbitrary matrices. For Hermitian matrices, we propose the generalised quantum eigenvalue transformation that even allows polynomials of indefinite parity. While we find that the polynomial might have to be downscaled compared to the quantum singular value transformation, the higher expressivity of polynomials and faster computation of phase factors can sometimes result in advantages. The results are achieved with various block encoding (or projected unitary encoding) techniques, including qubitisation, Hermitianisation, and multiplication. We show how to multiply block-encoded matrices with only one extra qubit, and introduce measure-early multiplication to further avoid the extra qubit and decrease average circuit length.
Communication complexity quantifies how difficult it is for two distant computers to evaluate a function $f(X,Y)$ where the strings $X$ and $Y$ are distributed to the first and second computer, respectively and under the constraint of exchanging as few bits as possible. Surprisingly, some nonlocal boxes, which are resources shared by the two computers, are so powerful that they allow to collapse communication complexity, in the sense that any Boolean function $f$ can be correctly estimated with the exchange of only one bit of communication. The Popescu-Rohrlich (PR) box is an example of such a collapsing resource, but a comprehensive description of the set of collapsing nonlocal boxes remains elusive. In this work, we carry out an algebraic study of the structure of wirings connecting nonlocal boxes, thus defining the notion of the "product of boxes" $\mathtt{P}\boxtimes\mathtt{Q}$, and we show related associativity and commutativity results. This gives rise to the notion of the "orbit of a box", unveiling surprising geometrical properties about the alignment and parallelism of distilled boxes. The power of this new framework is that it allows to prove previously-reported numerical intuitions concerning the best way to wire consecutive boxes, and to numerically and analytically recover recently-identified noisy PR boxes that collapse communication complexity for different types of noise models.
In this paper, we explore the impact of noise on quantum computing, particularly focusing on the challenges when sampling bit strings from noisy quantum computers as well as the implications for optimization and machine learning applications. We formally quantify the sampling overhead to extract good samples from noisy quantum computers and relate it to the layer fidelity, a metric to determine the performance of noisy quantum processors. Further, we show how this allows us to use the Conditional Value at Risk of noisy samples to determine provable bounds on noise-free expectation values. We discuss how to leverage these bounds for different algorithms and demonstrate our findings through experiments on a real quantum computer involving up to 127 qubits. The results show a strong alignment with theoretical predictions.
Energy filter methods in combination with quantum simulation can efficiently access the properties of quantum many-body systems at finite energy densities [Lu et al. PRX Quantum 2, 020321 (2021)]. Classically simulating this algorithm with tensor networks can be used to investigate the microcanonical properties of large spin chains, as recently shown in [Yang et al. Phys. Rev. B 106, 024307 (2022)]. Here we extend this strategy to explore the properties of off-diagonal matrix elements of observables in the energy eigenbasis, fundamentally connected to the thermalization behavior and the eigenstate thermalization hypothesis. We test the method on integrable and non-integrable spin chains of up to 60 sites, much larger than accessible with exact diagonalization. Our results allow us to explore the scaling of the off-diagonal functions with the size and energy difference, and to establish quantitative differences between integrable and non-integrable cases
Parametric amplifiers play a crucial role in modern quantum technology by enabling the enhancement of weak signals with minimal added noise. Traditionally, Josephson junctions have been the primary choice for constructing parametric amplifiers. Nevertheless, high-kinetic inductance thin films have emerged as viable alternatives to engineer the necessary nonlinearity. In this work, we introduce and characterize a Kinetic Inductance Parametric Amplifier (KIPA) built using high-quality NbN superconducting thin films. The KIPA addresses some of the limitations of traditional Josephson-based parametric amplifiers, excelling in dynamic range, operational temperature, and magnetic field resilience. We demonstrate a quantum-limited amplification (> 20 dB) with a 20 MHz gain-bandwidth product, operational at fields up to 6 Tesla and temperatures as high as 850 mK. Harnessing kinetic inductance in NbN thin films, the KIPA emerges as a robust solution for quantum signal amplification, enhancing research possibilities in quantum information processing and low-temperature quantum experiments. Its magnetic field compatibility and quantum-limited performance at high temperatures make it an invaluable tool, promising new advancements in quantum research.
We study superfluid properties of alkali-earth-like Fermi atomic systems in the presence of orbital Feshbach resonance. Using a two-band description of the ground state and excited state and a mean-field approximation of the intra-band atomic pairing, we investigate the phase transitions and crossover between superfluid/normal phases. Defining an effective scattering length by combining both inter-band and intra-band interactions, we derive closed form gap and number density equations for both ground state and excited state atomic bands. We find that our zero-temperature analytical results and finite-temperature numerical results indicate that the system can show smooth crossover between Bardeen, Cooper, and Schreifer (BCS) and Bose-Einstein Condensate (BEC) superfluidity for atoms in each band. In addition, we find that inter-band and intra-band interactions can induce quantum phase transitions between BCS/BEC superfluid states of atoms in one band to that of the other. We anticipate that our closed form analytical results can be used as a bench mark for future experimental and theoretical investigations and will have an impact on the current understanding of two-band superconductors such as MgB$_2$.
While the concept of the entanglement spectrum has heretofore been utilized to address various many-body systems, the models describing an itinerant spinless-fermion excitation coupled to zero-dimensional bosons (e.g. dispersionless phonons) have as yet not received much attention in this regard. To fill this gap, the ground-state entanglement spectrum of a model that includes two of the most common types of short-ranged, nonlocal excitation-phonon interaction -- the Peierls- and breathing-mode couplings -- is numerically evaluated here. This model displays a sharp, level-crossing transition at a critical coupling strength, which signifies the change from a nondegenerate ground state at the quasimomentum $K_{\textrm{gs}}=0$ to a twofold-degenerate one corresponding to a symmetric pair of nonzero quasimomenta. Another peculiarity of this model is that in the special case of equal Peierls- and breathing-mode coupling strengths the bare-excitation Bloch state with the quasimomentum $0$ or $\pi$ is its exact eigenstate. Moreover, below a critical coupling strength this state is the ground state of the model. Thus, the sharp transition between a bare excitation and a heavily phonon-dressed (polaronic) one can be thought of as a transition between vanishing and finite entanglement. It is demonstrated here that the smallest ground-state entanglement-spectrum eigenvalue to a large extent mimics the behavior of the entanglement entropy itself and vanishes in this special case of the model; by contrast, all the remaining eigenvalues diverge in this case. The implications of excitation-phonon entanglement for $W$-state engineering in superconducting and neutral-atom-based qubit arrays serving as analog simulators of this model are also discussed.
We present an efficient three-dimensional dark-focus optical trapping potential for neutral atoms and Bose-Einstein condensates. This "optical bottle" is created by a single blue-detuned light field exploiting the phenomenon of conical refraction occurring in biaxial crystals. The conversion of a Gaussian input beam to the bottle beam has an efficiency of close to 100 % and the optical setup requires the addition of the biaxial crystal and a circular polarizer only. Based on the conical-refraction theory, we derive the general form of the potential, the trapping frequencies, and the potential barrier heights. We present experiments on confining a $^{87}$Rb Bose-Einstein condensate in three dimensions. We determine the trap shape, the vibrational frequencies along the weak axis, as well as the lifetime of ultracold atoms in this type of potential.
The search for topological excitations such as Majorana fermions has spurred interest in the boundaries between distinct quantum states. Here, we explore an interface between two prototypical phases of electrons with conceptually different ground states: the integer quantum Hall insulator and the s-wave superconductor. We find clear signatures of hybridized electron and hole states similar to chiral Majorana fermions, to which we refer as chiral Andreev edge states (CAES). They propagate along the interface in the direction determined by magnetic field and their interference can turn an incoming electron into an outgoing electron or a hole, depending on the phase accumulated by the CAES along their path. Our results demonstrate that these excitations can propagate and interfere over a significant length, opening future possibilities for their coherent manipulation.
We have investigated the heralded generation of two-qubit dual-rail-encoded states by programmable linear optics. Two types of schemes generating the states from four single photons, which is the minimal possible to accomplish the task, have been considered. The schemes have different detection patterns heralding successful generation events, namely, one-mode heralding, in which the two auxiliary photons are detected in one mode, and two-mode heralding, in which single photons are detected in each of the two modes simultaneously. We have shown that the dependence of the schemes' success probabilities on the target state's degree of entanglement are essentially different. In particular, one-mode heralding yields better efficiency for highly-entangled states, if the programmable interferometers can explore the full space of the unitary transfer matrices,. It is reversed in case of weakly-entangled states where two-mode heralding is better. We have found a minimal decomposition of the scheme with two-mode heralding that is programmed by one variable phase shift. We infer that the linear optical schemes designed specifically for generation of two-qubit states are more efficient than schemes implementing gate-based circuits with known two-qubit linear optical gates. Our results yield substantial reduction of physical resources needed to generate two-qubit dual-rail-encoded photonic states.
The vanishing band gap of graphene has long presented challenges for making high-quality quantum point contacts (QPCs) -- the partially transparent p-n interfaces introduced by conventional split-gates tend to short the QPC. This complication has hindered the fabrication of graphene quantum Hall Fabry-P\'erot interferometers, until recent advances have allowed split-gate QPCs to operate utilizing the highly resistive $\nu=0$ state. Here, we present a simple recipe to fabricate QPCs by etching a narrow trench in the graphene sheet to separate the conducting channel from self-aligned graphene side gates. We demonstrate operation of the individual QPCs in the quantum Hall regime, and further utilize these QPCs to create and study a quantum Hall interferometer.
The speed of elementary quantum gates, particularly two-qubit gates, ultimately sets the limit on the speed at which quantum circuits can operate. In this work, we experimentally demonstrate commonly used two-qubit gates at nearly the fastest possible speed allowed by the physical interaction strength between two superconducting transmon qubits. We achieve this quantum speed limit by implementing experimental gates designed using a machine learning inspired optimal control method. Importantly, our method only requires the single-qubit drive strength to be moderately larger than the interaction strength to achieve an arbitrary two-qubit gate close to its analytical speed limit with high fidelity. Thus, the method is applicable to a variety of platforms including those with comparable single-qubit and two-qubit gate speeds, or those with always-on interactions. We expect our method to offer significant speedups for non-native two-qubit gates that are typically achieved with a long sequence of single-qubit and native two-qubit gates.
We propose an experimental method for evaluating the adiabatic condition during quantum annealing (QA), which will be essential for solving practical problems. The adiabatic condition consists of the transition matrix element and the energy gap, and our method simultaneously provides information about these components without diagonalizing the Hamiltonian. The key idea is to measure the power spectrum of a time domain signal by adding an oscillating field during QA, and we can estimate the values of the transition matrix element and energy gap from the measurement output. Our results provides a powerful experimental basis for analyzing the performance of QA.
The information-theoretic unconditional security offered by quantum key distribution has spurred the development of larger quantum communication networks. However, as these networks grow so does the strong need to reduce complexity and overheads. Polarization based entanglement distribution networks are a promising approach due to their scalability and lack of trusted nodes. Nevertheless, they are only viable if the birefringence of all optical distribution fibres in the network is compensated to preserve the polarization based quantum state. The brute force approach would require a few hundred fibre polarization controllers for even a moderately sized network. Instead, we propose and investigate four different methods of polarization compensation. We compare them based on complexity, effort, level of disruption to network operations and performance.
Spectral characterization of noise environments that lead to the decoherence of qubits is critical to developing robust quantum technologies. While dynamical decoupling offers one of the most successful approaches to characterize noise spectra, it necessitates applying large sequences of $\pi$ pulses that increase the complexity and cost of the method. Here, we introduce a noise spectroscopy method that utilizes only the Fourier transform of free induction decay or spin echo measurements, thus removing the need for the application many $\pi$ pulses. We show that our method faithfully recovers the correct noise spectra for a variety of different environments (including $1/f$-type noise) and outperforms previous dynamical decoupling schemes while significantly reducing their experimental overhead. We also discuss the experimental feasibility of our proposal and demonstrate its robustness in the presence of statistical measurement error. Our method is applicable to a wide range of quantum platforms and provides a simpler path toward a more accurate spectral characterization of quantum devices, thus offering possibilities for tailored decoherence mitigation.
Creating a massive spatial quantum superposition, such as the Schr\"odinger cat state, where the mass and the superposition size within the range $10^{-19}-10^{-14}$ kg and $\Delta x \sim 10~{\rm nm}-100~\mu {\rm m}$, is a challenging task. The methods employed so far rely either on wavepacket expansion or on a quantum ancilla, e.g. single spin dependent forces, which scale inversely with mass. In this paper, we present a novel approach that combines gravitational acceleration and diamagnetic repulsion to generate a large spatial superposition in a relatively short time. After first creating a modest initial spatial superposition of $1~\mu {\rm m}$, achieved through techniques such as the Stern-Gerlach (SG) apparatus, we will show that we can achieve an $\sim 10^{2}-10^{3}$ fold improvement to the spatial superposition size ($1~{\rm \mu m}\rightarrow 980~\mu {\rm m}$) between the wave packets in less than $0.02$~s by using the Earth's gravitational acceleration and then the diamagnetic repulsive scattering of the nanocrystal, neither of which depend on the object mass. Finally, the wave packet trajectories can be closed so that spatial interference fringes can be observed. Our findings highlight the potential of combining gravitational acceleration and diamagnetic repulsion to create and manipulate large spatial superpositions, offering new insights into creating macroscopic quantum superpositions.
We analyze theoretically the motional quantum dynamics of a levitated dielectric sphere interacting with the quantum electromagnetic field beyond the point-dipole approximation. To this end, we derive a Hamiltonian describing the fundamental coupling between photons and center-of-mass phonons, including Stokes and anti-Stokes processes, and the coupling rates for a dielectric sphere of arbitrary refractive index and size. We then derive the laser recoil heating rates and the information radiation patterns (the angular distribution of the scattered light that carries information about the center-of-mass motion) and show how to evaluate them efficiently in the presence of a focused laser beam, in either a running- or a standing-wave configuration. This information is crucial to implement active feedback cooling of optically levitated dielectric spheres beyond the point-dipole approximation. Our results predict several experimentally feasible configurations and parameter regimes where optical detection and active feedback can simultaneously cool to the ground state the three-dimensional center-of-mass motion of dielectric spheres in the micrometer regime. Scaling up the mass of the dielectric particles that can be cooled to the center-of-mass ground state is relevant not only for testing quantum mechanics at large scales but also for current experimental efforts that search for new physics (e.g., dark matter) using optically levitated sensors.
Continuous-variable quantum computing architectures based upon the Gottesmann-Kitaev-Preskill (GKP) encoding have emerged as a promising candidate because one can achieve fault-tolerance with a probabilistic supply of GKP states and Gaussian operations. Furthermore, by generalising to rectangular-lattice GKP states, a bias can be introduced and exploited through concatenation with qubit codes that show improved performance under biasing. However, these codes (such as the XZZX surface code) still require weight-four stabiliser measurements and have complex decoding requirements to overcome. In this work, we study the code-capacity behaviour of a rectangular-lattice GKP encoding concatenated with a repetition code under an isotropic Gaussian displacement channel. We find a numerical threshold of $\sigma = 0.599$ for the noise's standard deviation, which outperforms the biased GKP planar surface code with a trade-off of increased biasing at the GKP level. This is all achieved with only weight-two stabiliser operators and simple decoding at the qubit level. Furthermore, with moderate levels of bias (aspect ratio $\leq 2.4$) and nine or fewer data modes, significant reductions in logical error rates can still be achieved for $\sigma \leq 0.3$, opening the possibility of using GKP-biased repetition codes as a simple low-level qubit encoding for further concatenation.
We investigate the full counting statistics of charge transport in $U(1)$-symmetric random unitary circuits. We consider an initial mixed state prepared with a chemical potential imbalance between the left and right halves of the system, and study the fluctuations of the charge transferred across the central bond in typical circuits. Using an effective replica statistical mechanics model and a mapping onto an emergent classical stochastic process valid at large onsite Hilbert space dimension, we show that charge transfer fluctuations approach those of the symmetric exclusion process at long times, with subleading $t^{-1/2}$ quantum corrections. We discuss our results in the context of fluctuating hydrodynamics and macroscopic fluctuation theory of classical non-equilibrium systems, and check our predictions against direct matrix-product state calculations.
Although quantum circuits have been ubiquitous for decades in quantum computing, the first complete equational theory for quantum circuits has only recently been introduced. Completeness guarantees that any true equation on quantum circuits can be derived from the equational theory. We improve this completeness result in two ways: (i) We simplify the equational theory by proving that several rules can be derived from the remaining ones. In particular, two out of the three most intricate rules are removed, the third one being slightly simplified. (ii) The complete equational theory can be extended to quantum circuits with ancillae or qubit discarding, to represent respectively quantum computations using an additional workspace, and hybrid quantum computations. We show that the remaining intricate rule can be greatly simplified in these more expressive settings, leading to equational theories where all equations act on a bounded number of qubits. The development of simple and complete equational theories for expressive quantum circuit models opens new avenues for reasoning about quantum circuits. It provides strong formal foundations for various compiling tasks such as circuit optimisation, hardware constraint satisfaction and verification.
Atom-interferometric quantum sensors could revolutionize navigation, civil engineering, and Earth observation. However, operation in real-world environments is challenging due to external interference, platform noise, and constraints on size, weight, and power. Here we experimentally demonstrate that tailored light pulses designed using robust control techniques mitigate significant error sources in an atom-interferometric accelerometer. To mimic the effect of unpredictable lateral platform motion, we apply laser-intensity noise that varies up to 20$\%$ from pulse-to-pulse. Our robust control solution maintains performant sensing, while the utility of conventional pulses collapses. By measuring local gravity, we show that our robust pulses preserve interferometer scale factor and improve measurement precision by 10$\times$ in the presence of this noise. We further validate these enhancements by measuring applied accelerations over a 200 $\mu g$ range up to 21$\times$ more precisely at the highest applied noise level. Our demonstration provides a pathway to improved atom-interferometric inertial sensing in real-world settings.
Mitigation and calibration schemes are central to maximize the computational reach of today's Noisy Intermediate Scale Quantum (NISQ) hardware, but these schemes are often specialized to exclusively address either coherent or decoherent error sources. Quantifying the two types of errors hence constitutes a desirable feature when it comes to benchmarking error suppression tools. In this paper, we present a scalable and cycle-centric methodology for obtaining a detailed estimate of the coherent contribution to the error profile of a hard computing cycle. The protocol that we suggest is based on Cycle Error Reconstruction (CER), also known as K-body Noise Reconstruction (KNR). This protocol is similar to Cycle Benchmarking (CB) in that it provides a cycle-centric diagnostic based on Pauli fidelity estimation [1]. We introduce an additional hyper-parameter in CER by allowing the hard cycles to be folded multiple times before being subject to Pauli twirling. Performing CER for different values of our added hyper-parameter allows estimating the coherent error contributions through a generalization of the fidelity decay formula. We confirm the accuracy of our method through numerical simulations on a quantum simulator, and perform proof-of-concept experiments on three IBM chips, namely ibmq_guadalupe, ibmq_manila, and ibmq_montreal. In all three experiments, we measure substantial coherent errors biased in Z.
We study the photophysics of molecular aggregates from a quantum optics perspective, with emphasis on deriving scaling laws for the fast non-radiative relaxation of collective electronic excitations, referred to as Kasha's rule. Aggregates exhibit an energetically broad manifold of collective states with delocalized electronic excitations originating from near field dipole-dipole exchanges between neighboring monomers. Photo-excitation at optical wavelengths, much larger than the monomer-monomer average separation, addresses almost exclusively symmetric collective states, which for an arrangement known as H-aggregate, show an upward hypsochromic shift. The extremely fast subsequent non-radiative relaxation via intramolecular vibrational modes populates lower energy, subradiant states, resulting in an effective inhibition of fluorescence. Our analytical treatment allows for the derivation of an approximate scaling law of this relaxation process, linear in the number of available low energy vibrational modes and directly proportional to the dipole-dipole interaction strength between neighbouring monomers.
How multiple observables mutually influence their dynamics has been a crucial issue in statistical mechanics. We introduce a new concept, "quantum velocity limits," to establish a quantitative and rigorous theory for non-equilibrium quantum dynamics for multiple observables. Quantum velocity limits are universal inequalities for a vector the describes velocities of multiple observables. They elucidate that the speed of an observable of our interest can be tighter bounded when we have knowledge of other observables, such as experimentally accessible ones or conserved quantities, compared with the conventional speed limits for a single observable. We first derive an information-theoretical velocity limit in terms of the generalized correlation matrix of the observables and the quantum Fisher information. The velocity limit has various novel consequences: (i) conservation law in the system, a fundamental ingredient of quantum dynamics, can improve the velocity limits through the correlation between the observables and conserved quantities; (ii) speed of an observable can be bounded by a nontrivial lower bound from the information on another observable; (iii) there exists a notable non-equilibrium tradeoff relation, stating that speeds of uncorrelated observables, e.g., anti-commuting observables, cannot be simultaneously large; (iv) velocity limits for any observables on a local subsystem in locally interacting many-body systems remain convergent even in the thermodynamic limit. Moreover, we discover another distinct velocity limit for multiple observables on the basis of the local conservation law of probability current, which becomes advantageous for macroscopic transitions of multiple quantities.
We study the effect that the SU(2) symmetry, and the rich Hilbert space structure that it generates in lattice spin systems, has on the average entanglement entropy of highly excited eigenstates of local Hamiltonians and of random pure states. Focusing on the zero total magnetization sector ($J_z=0$) for different fixed total spin $J$, we argue that the average entanglement entropy of highly excited eigenstates of quantum-chaotic Hamiltonians and of random pure states has a leading volume-law term whose coefficient $s_A$ depends on the spin density $j=J/(\mathfrak{j}L)$, with $s_A(j \rightarrow 0)=\ln (2\mathfrak{j}+1)$ and $s_A(j \rightarrow 1)=0$, where $\mathfrak{j}$ is the microscopic spin. We provide numerical evidence that $s_A$ is smaller in highly excited eigenstates of integrable interacting Hamiltonians, which lends support to the expectation that the average eigenstate entanglement entropy can be used as a diagnostic of quantum chaos and integrability for Hamiltonians with non-Abelian symmetries. In the context of Hamiltonian eigenstates we consider spins $\mathfrak{j}=\frac12$ and $1$, while for our calculations based on random pure states we focus on the spin $\mathfrak{j}=\frac12$ case.
We present a first-principles study of defect formation and electronic structure of erbium (Er)-doped yttria (Y$_2$O$_3$). This is an emerging material for spin-photon interfaces in quantum information science due to the narrow linewidth optical emission from Er dopants at standard telecommunication wavelengths and their potential for quantum memories and transducers. We calculate formation energies of neutral, negatively, and positively charged Er dopants and find the charge neutral configuration to be the most stable, consistent with experiment. Of the two substitutional sites of Er for Y, the $C_2$ (more relevant for quantum memories) and $C_{3i}$ (more relevant for quantum transduction), we identify the former as possessing the lowest formation energy. The electronic properties are calculated using the Perdew-Burke-Ernzerhof (PBE) functional along with the Hubbard $U$ parameter and spin-orbit coupling (SOC), which yields a $\sim$ 6 $\mu_B$ orbital and a $\sim$ 3 $\mu_B$ spin magnetic moment, and 11 electrons in the Er $4f$ shell, confirming the formation of charge-neutral Er$^{3+}$. This standard density functional theory (DFT) approach underestimates the band gap of the host and lacks a first-principles justification for $U$. To overcome these issues, we performed screened hybrid functional (HSE) calculations, including a negative $U$ for the $4f$ orbitals, with mixing ($\alpha$) and screening ($w$) parameters. These produced robust electronic features with slight modifications in the band gap and the $4f$ splittings depending on the choice of tuning parameters. We also computed the many-particle electronic excitation energies and compared them with experimental values from photoluminescence.
Power system dynamics are generally modeled by high dimensional nonlinear differential-algebraic equations due to a large number of generators, loads, and transmission lines. Thus, its computational complexity grows exponentially with the system size. In this paper, we aim to evaluate the alternative computing approach, particularly the use of quantum computing algorithms to solve the power system dynamics. Leveraging a symbolic programming framework, we convert the power system dynamics' DAEs into an equivalent set of ordinary differential equations (ODEs). Their data can be encoded into quantum computers via amplitude encoding. The system's nonlinearity is captured by Taylor polynomial expansion and the quantum state tensor whereas state variables can be updated by a quantum linear equation solver. Our results show that quantum computing can solve the dynamics of the power system with high accuracy whereas its complexity is polynomial in the logarithm of the system dimension.
We describe a numerical time-domain approach for high-accuracy calculations of Casimir-Polder forces near micro-structured materials. The use of a time-domain formulation enables the investigation of a broad range of materials described by advanced material models, including nonlocal response functions. We validate the method by a number of example calculations for which we thoroughly investigate the convergence properties of the method, and comparing to analytical reference calculations, we find average relative errors as low as a few parts in a million. As an application example, we investigate the anisotropy-induced repulsive behavior of the Casimir-Polder force near a sharp gold wedge described by a hydrodynamic Drude model.
Superconducting nanowire single-photon detectors (SNSPDs) have been widely used to study the discrete nature of quantum states of light in the form of photon-counting experiments. We show that SNSPDs can also be used to study continuous variables of optical quantum states by performing homodyne detection at a bandwidth of $400~\mathrm{kHz}$. By measuring the interference of a continuous-wave field of a local oscillator with the field of the vacuum state using two SNSPDs, we show that the variance of the difference in count rates is linearly proportional to the photon flux of the local oscillator over almost five orders of magnitude. The resulting shot-noise clearance of $(46.0\pm1.1)~\mathrm{dB}$ is the highest reported clearance for a balanced optical homodyne detector, demonstrating their potential for measuring highly squeezed states in the continuous-wave regime. In addition, we measured a $\mathrm{CMRR}=22.4~\mathrm{dB}$. From the joint click counting statistics, we also measure the phase-dependent quadrature of a weak coherent state to demonstrate our device's functionality as a homodyne detector.
Quantum contextuality describes scenarios in which it is impossible to explain the experimental evidence in terms of a measurement independent reality. Here, I introduce a three-path interferometer in which all five contexts needed for a demonstration of contextuality are realized in sequence. It is then possible to observe a paradoxical situation where the paths connecting input ports to their corresponding output ports appear to be blocked by destructive interference. It is shown that the conditional currents observed in weak measurements provide a consistent explanation of the paradox, indicating that weak values might help to bridge the gap between wavelike propagation effects and local particle detection.
The probability distribution of the occupancy numbers of a system at the equilibrium composed by an arbitrary number of non-interacting bosons is obtained by tracing out the environment from the ''universe'', that is the union of environment and system of interest. The first new result presented in the paper is that, when the Bayesian prior of the universe is the multinomial distribution, also the marginal of the system of interest is the multinomial distribution. This self-consistency reveals the coherence of the Bayesian-multinomial approach. The second new result presented in the paper is that, when the universe is in a bosonic eigenstate, the distribution of the system of interest is the multivariate hypergeometric distribution and it is well-known that, when the number of bosons of the universe tends to infinity, the multivariate hypergeometric distribution tends to the multinomial distribution. Furthermore, the paper proposes to identify the physical entropy of the bosonic system with the Shannon entropy of the occupancy numbers, fixing certain contradictions that arise in the classical analysis of thermodynamic entropy. Finally, by leveraging an information-theoretic inequality between the entropy of the multinomial distribution and the entropy of the multivariate hypergeometric distribution, both the Bayesianism of information theory and the empiricism of statistical mechanics are integrated into a common ''infomechanical'' framework.
Here we revisit one of the prototypical tasks for characterizing the structure of noise in quantum devices: estimating every eigenvalue of an $n$-qubit Pauli noise channel to error $\epsilon$. Prior work (Chen et al., 2022) proved no-go theorems for this task in the practical regime where one has a limited amount of quantum memory, e.g. any protocol with $\le 0.99n$ ancilla qubits of quantum memory must make exponentially many measurements, provided it is non-concatenating. Such protocols can only interact with the channel by repeatedly preparing a state, passing it through the channel, and measuring immediately afterward. This left open a natural question: does the lower bound hold even for general protocols, i.e. ones which chain together many queries to the channel, interleaved with arbitrary data-processing channels, before measuring? Surprisingly, in this work we show the opposite: there is a protocol that can estimate the eigenvalues of a Pauli channel to error $\epsilon$ using only $O(\log n/\epsilon^2)$ ancilla qubits and $\tilde{O}(n^2/\epsilon^2)$ measurements. In contrast, we show that any protocol with zero ancilla, even a concatenating one, must make $\Omega(2^n/\epsilon^2)$ measurements, which is tight. Our results imply, to our knowledge, the first quantum learning task where logarithmically many qubits of quantum memory suffice for an exponential statistical advantage.
The Quantum Approximate Optimization Algorithm (QAOA) is one of the most promising Noisy Intermediate Quantum Algorithms (NISQ) in solving combinatorial optimizations and displays potential over classical heuristic techniques. Unfortunately, QAOA performance depends on the choice of parameters and standard optimizers often fail to identify key parameters due to the complexity and mystery of these optimization functions. In this paper, we benchmark QAOA circuits modified with metaheuristic optimizers against classical and quantum heuristics to identify QAOA parameters. The experimental results reveal insights into the strengths and limitations of both Quantum Annealing and metaheuristic-integrated QAOA across different problem domains. The findings suggest that the hybrid approach can leverage classical optimization strategies to enhance the solution quality and convergence speed of QAOA, particularly for problems with rugged landscapes and limited quantum resources. Furthermore, the study provides guidelines for selecting the most appropriate approach based on the specific characteristics of the optimization problem at hand.
By replacing the spatial derivative with the Dunkl derivative, we generalize the Fokker-Planck equation in (1+1) dimensions. We obtain the Dunkl-Fokker-Planck eigenvalues equation and solve it for the harmonic oscillator plus a centrifugal-type potential. Furthermore, when the drift function is odd, we reduce our results to those of the recently developed Wigner-Dunkl supersymmetry.
Introductory courses on quantum mechanics usually include lectures on uncertainty relations, typically the inequality derived by Robertson and, perhaps, other statements. For the benefit of the lecturers, we present a unified approach -- well suited for undergraduate teaching -- for deriving all standard uncertainty relations: those for products of variances by Kennard, Robertson, and Schr\"odinger, as well as those for sums of variances by Maccone and Pati. We also give a brief review of the early history of this topic and try to answer why the use of variances for quantifying uncertainty is so widespread, while alternatives are available that can be more natural and more fitting. It is common to regard the states that saturate the Robertson inequality as "minimum uncertainty states" although they do not minimize the variance of one observable, given the variance of another, incompatible observable. The states that achieve this objective are different and can be found systematically.
We present a quantum algorithm that analyzes time series data simulated by a quantum differential equation solver. The proposed algorithm is a quantum version of the dynamic mode decomposition algorithm used in diverse fields such as fluid dynamics and epidemiology. Our quantum algorithm can also extract matrix eigenvalues by analyzing the corresponding linear dynamical system. Our algorithm handles a broader range of matrices with complex eigenvalues, unlike existing efficient quantum eigensolvers limited to specific matrix types. The complexity of our quantum algorithm is $O(\operatorname{poly}\log N)$ for an $N$-dimensional system. This is an exponential speedup over known classical algorithms with at least $O(N)$ complexity. Thus, our quantum algorithm is expected to enable high-dimensional dynamical systems analysis and large matrix eigenvalue decomposition, intractable for classical computers.
Color centers at the low-dispersion O-band wavelengths are an essential resource for long-lifetime quantum network nodes toward memory-assisted quantum communications using energy-time entanglement. In this work, we explore the process of developing T centers and other color center defects to improve qubit storage and radiative efficiency while examining the photoluminescence dynamics. We have extended the $TX_{0}$ lifetime of T centers by 65% to 1.56 $\mu$s. Furthermore, we discover the presence of a $^*Cu_n^m$ related doublet emission around 1312 nm close to the zero-dispersion wavelength, with a spin degeneracy resulting in a magnetic-field induced broadening by 25% under 0.5 T, which can be an alternative to T centers as a high-fidelity spin-photon interface.
We give a complete classification of the anyon sectors of Kitaev's quantum double model on the infinite triangular lattice and for finite gauge group $G$, including the non-abelian case. As conjectured, the anyon sectors of the model correspond precisely to the irreducible representations of the quantum double algebra of $G$. Our proof consists of two main parts. In the first part, we construct for each irreducible representation of the quantum double algebra a pure state and show that the GNS representations of these pure states are pairwise disjoint anyon sectors. In the second part we show that any anyon sector is unitarily equivalent to one of the anyon sectors constructed in the first part. Purity of the states constructed in the first part is shown by characterising these states as the unique states that satisfy appropriate local constraints. These constraints are of two types, namely flux constraints and gauge constraints. The flux constraints single out certain string-net states, while the gauge constraints fix the way in which these string-nets condense. At the core of the proof is the fact that certain groups of local gauge transformations act freely and transitively on collections of local string-nets. The proof that the GNS representations of these states are anyon sectors relies on showing that they are unitarily equivalent to amplimorphism representations which are much easier to compare to the ground state representation. For the second part, we show that any anyon sector contains a pure state that satisfies all but a finite number of the constraints characterising the pure states of the first part. Using known techniques we can then construct a pure state in the anyon sector that satisfies all but one of these constraints. Finally, we show that any such state must be a vector state in one of the anyon sectors constructed in the first part.
Open quantum many-body systems with controllable dissipation can exhibit novel features in their dynamics and steady states. A paradigmatic example is the dissipative transverse field Ising model. It has been shown recently that the steady state of this model with all-to-all interactions is genuinely non-equilibrium near criticality, exhibiting a modified time-reversal symmetry and violating the fluctuation-dissipation theorem. Experimental study of such non-equilibrium steady-state phase transitions is however lacking. Here we propose realistic experimental setups and measurement schemes for current trapped-ion quantum simulators to demonstrate this phase transition, where controllable dissipation is engineered via a continuous weak optical pumping laser. With extensive numerical calculations, we show that strong signatures of this dissipative phase transition and its non-equilibrium properties can be observed with a small system size across a wide range of system parameters. In addition, we show that the same signatures can also be seen if the dissipation is instead achieved via Floquet dynamics with periodic and probabilistic resetting of the spins. Dissipation engineered in this way may allow the simulation of more general types of driven-dissipative systems or facilitate the dissipative preparation of useful many-body entangled states.
Machine learning tasks are an exciting application for quantum computers, as it has been proven that they can learn certain problems more efficiently than classical ones. Applying quantum machine learning algorithms to classical data can have many important applications, as qubits allow for dealing with exponentially more data than classical bits. However, preparing the corresponding quantum states usually requires an exponential number of gates and therefore may ruin any potential quantum speedups. Here, we show that classical data with a sufficiently quickly decaying Fourier spectrum after being mapped to a quantum state can be well-approximated by states with a small Schmidt rank (i.e., matrix product states) and we derive explicit error bounds. These approximated states can, in turn, be prepared on a quantum computer with a linear number of nearest-neighbor two-qubit gates. We confirm our results numerically on a set of $1024\times1024$-pixel images taken from the 'Imagenette' dataset. Additionally, we consider different variational circuit ans\"atze and demonstrate numerically that one-dimensional sequential circuits achieve the same compression quality as more powerful ans\"atze.
In this work, we employ differential evolution algorithm to identify the optimal configurations of small atomic ensembles supporting quantum states with maximal radiative lifetime. We demonstrate that atoms mostly tend to assemble in quasi-regular structures with specific geometry strongly depending on the minimal interatomic distance $r_{min}$. We identified the clear underlying physics that governs the suppression of the radiative losses in particular geometries. However, we reveal that the specific configurations in small ensembles are not easily predictable based on the knowledge established for the arrays of large size. In particular, the states that inherit their properties from bound states in continuum in infinite lattices turn out to be the most subradiant in a wide range of $r_{min}$ values. We also show that for small interatomic distance the chains with modulated interatomic distances exhibit fast exponential decrease of the radiative losses with the size of the ensemble.
Coherently dressed spins have shown promising results as building blocks for future quantum computers owing to their resilience to environmental noise and their compatibility with global control fields. This mode of operation allows for more amenable qubit architecture requirements and simplifies signal routing on the chip. However, multi-qubit operations, such as qubit addressability and two-qubit gates, are yet to be demonstrated to establish global control in combination with dressed qubits as a viable path to universal quantum computing. Here we demonstrate simultaneous on-resonance driving of degenerate qubits using a global field while retaining addressability for qubits with equal Larmor frequencies. Furthermore, we implement SWAP oscillations during on-resonance driving, constituting the demonstration of driven two-qubit gates. Significantly, our findings highlight the fragility of entangling gates between superposition states and how dressing can increase the noise robustness. These results represent a crucial milestone towards global control operation with dressed qubits. It also opens a door to interesting spin physics on degenerate spins.
In this study, we present Q-Seg, a novel unsupervised image segmentation method based on quantum annealing, tailored for existing quantum hardware. We formulate the pixel-wise segmentation problem, which assimilates spectral and spatial information of the image, as a graph-cut optimization task. Our method efficiently leverages the interconnected qubit topology of the D-Wave Advantage device, offering superior scalability over existing quantum approaches and outperforming state-of-the-art classical methods. Our empirical evaluations on synthetic datasets reveal that Q-Seg offers better runtime performance against the classical optimizer Gurobi. Furthermore, we evaluate our method on segmentation of Earth Observation images, an area of application where the amount of labeled data is usually very limited. In this case, Q-Seg demonstrates near-optimal results in flood mapping detection with respect to classical supervised state-of-the-art machine learning methods. Also, Q-Seg provides enhanced segmentation for forest coverage compared to existing annotated masks. Thus, Q-Seg emerges as a viable alternative for real-world applications using available quantum hardware, particularly in scenarios where the lack of labeled data and computational runtime are critical.
Interacting lattice Hamiltonians at high temperature generically give rise to energy transport governed by the classical diffusion equation; however, predicting the rate of diffusion requires numerical simulation of the microscopic quantum dynamics. For the purpose of predicting such transport properties, computational time evolution methods must be paired with schemes to control the growth of entanglement to tractably simulate for sufficiently long times. One such truncation scheme -- dissipation-assisted operator evolution (DAOE) -- controls entanglement by damping out components of operators with large Pauli weight. In this paper, we generalize DAOE to treat fermionic systems. Our method instead damps out components of operators with large fermionic weight. We investigate the performance of DAOE, the new fermionic DAOE (FDAOE), and another simulation method, density matrix truncation (DMT), in simulating energy transport in an interacting one-dimensional Majorana chain. The chain is found to have a diffusion coefficient scaling like interaction strength to the fourth power, contrary to naive expectations based on Fermi's Golden rule -- but consistent with recent predictions based on the theory of \emph{weak integrability breaking}. In the weak interaction regime where the fermionic nature of the system is most relevant, FDAOE is found to simulate the system more efficiently than DAOE.
Building on decades of fundamental research, new applications of Quantum Science have started to emerge in the fields of computing, sensing and networks. In the current phase of deployment, in which quantum technology is not yet in routine use but is still transitioning out of the laboratory, Venture Capital (VC) is critical. In association with public funding programs, VC supports startups born in academic institutions and has a role to play in structuring the priorities of the ecosystem, guiding it towards applications with the greatest impact on society. This paper illustrates this thesis with a case-study: the experience of the first dedicated quantum fund, Quantonation I, chronicling its impacts on the production of scientific knowledge, job creation and funding of the industry. The paper introduces concepts to support the emergence of new startups and advocates for funding of scale-up quantum companies. The paper concludes with proposals to improve the impact of the industry by taking steps to better involve society-at-large and with a call for collaboration on projects focused on the applications with a large societal benefit.
We examine the production of the Hoyle and associated excited states from the viewpoint of pocket resonances in the reaction of an $\alpha$-particle on a ground state prolate $^8$Be nucleus within the optical model coupled-channel framework. The predicted reaction cross sections, as a function of the center-of-mass energy $E_{\rm cm}$, show prominent resonances, including the Hoyle resonance. The positions and widths of these resonances are sensitive to the target deformation ($\beta_2$ parameter) and the parity of the nuclear surface potential $-$ deeper for the even-parity $L$ partial waves relative to those for the odd-parity $L$ partial waves at the surface region because of the Bose-Einstein exchange of the $\alpha$-bosons. Decomposing the reaction cross sections to different partial waves, we find that the resonance energies and widths reasonably agree with the available experimental data and previous hyperspherical calculations for the $0_2^+$ (Hoyle state), $0_3^+$, $1_1^-$ and $3_1^-$ states of $^{12}$C, except for the narrow theoretical width of the $2_2^+$ state. Analyzing the wavefunctions and the resonance widths, we identify the narrow and sharp $0_2^+$, $3_1^-$ and $2_2^+$ resonances as pocket resonances -- resonances which occur below the potential barrier, while the broad $0_3^+$ and $1_1^-$ resonances as above-the-barrier resonances. For astrophysical applications, we also evaluate the astrophysical $S(E_{\rm cm})$-factor for $E_{\rm cm}$ $<$ 1.0 MeV, for the fusion of $\alpha$+$^8$Be into the $^{12}$C$(2^+)$ state based on our estimated $s$-wave $\alpha$+$^8$Be reaction cross section and the associated $\gamma$- and $\alpha$-decay widths for the decay of $^{12}$C excited states in the potential pocket.