Search NASA⌕ Search

SEARCH · Search NASA

Results for “correlation analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Dynamical properties of a room temperature ionic liquid: Using molecular dynamics simulations to implement a dynamic ion cage model [plus supplemental information]

The transport behavior of ionic liquids (ILs) is pivotal for a variety of applications, especially when ILs are used as electrolytes. Many aspects of the transport dynamics of ILs remain to be understood. Here, a common ionic liquid, 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide (BmimNTf 2 ), was studied with molecular dynamics simulations. The results show that BmimNTf 2 displays typical structural relaxation, subdiffusive behavior, and a breakdown of the Stokes-Einstein diffusion relation as in glass-forming liquids. In addition, the simulations show that the translational dynamics, reorientation dynamics, and structural relaxation dynamics are well described by the Vogel-Fulcher-Tammann equation like fragile glass forming liquids. Building on previous work that employed ion cage models, it was found that the diffusion dynamics of the cations and anions were well described by a hopping process random walk where the step time is the ion cage lifetime obtained from the cage correlation function. Detailed analysis of the ion cage structures indicated that the electrostatic potential energy of the ion cage dominates the diffusion dynamics of the caged ion. The ion orientational relaxation dynamics showed that ion reorientation is a necessary step for ion cage restructuring. The dynamic ion cage model description of ion diffusion presented here may have implications for designing ILs to control their transport behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Effective Hamiltonians derived from equation-of-motion coupled-cluster wave functions: Theory and application to the Hubbard and Heisenberg Hamiltonians

Effective Hamiltonians, which are commonly used for fitting experimental observables, provide a coarse-grained representation of exact many-electron states obtained in quantum chemistry calculations; however, the mapping between the two is not trivial. In this contribution, we apply Bloch’s formalism to equation-of-motion coupled-cluster (EOM-CC) wave functions to rigorously derive effective Hamiltonians in the Bloch’s and des Cloizeaux’ forms. We report the key equations and illustrate the theory by application to systems with two or three unpaired electrons, which give rise to electronic states of covalent and ionic characters. We show that the Hubbard and Heisenberg Hamiltonians can be extracted directly from the so-obtained effective Hamiltonians. By establishing a quantitative connection between many-body states and simple models, the approach facilitates the analysis of the correlated wave functions. We propose a simple diagnostic for assessing the validity of the model space choice based on the overlaps between the target- and model-space states. Finally, artifacts affecting the quality of electronic structure calculations such as spin contamination are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Computational evolution of high-performing unfused non-fullerene acceptors for organic solar cells

Materials optimization for organic solar cells (OSCs) is a highly active field, with many approaches using empirical experimental synthesis, computational brute force to screen a subset of chemical space, or generative machine learning methods that often require significant training sets. While these methods may find high-performing materials, they can be inefficient and time-consuming. Genetic algorithms (GAs) are an alternative approach, allowing for the “virtual synthesis” of molecules and a prediction of their “fitness” for some property, with new candidates suggested based on good characteristics of previously generated molecules. In this work, a GA is used to discover high-performing unfused non-fullerene acceptors (NFAs) based on an empirical prediction of power conversion efficiency (PCE) and provides design rules for future work. The electron-withdrawing/donating strength, as well as the sequence and symmetry, of those units are examined. The utilization of a GA over a brute-force approach resulted in speedups up to 1.8 × 10 12 . New types of units, not frequently seen in OSCs, are suggested, and in total 5426 NFAs are discovered with the GA. Of these, 1087 NFAs are predicted to have a PCE greater than 18%, which is roughly the current record efficiency. While the symmetry of the sequence showed no correlation with PCE, analysis of the sequence arrangement revealed that higher performance can be achieved with a donor core and acceptor end groups. Future NFA designs should consider this strategy as an alternative to the current A-D-A'-D-A architecture.

14 SOLAR ENERGY↗
DOE OSTI2022

Installation of a solid state neutral particle analyzer array on mega ampere spherical tokamak upgrade

A compact solid state neutral particle analyzer (SSNPA) diagnostic, previously installed at NSTX-U, has been moved to MAST-U and successfully operated in the first physics campaign (MU01). The SSNPA operates by detecting the flux of fast neutral particles produced by charge exchange (CX) reactions to diagnose the fast ion distribution. The diagnostic consists of three 16-channel sensors, which provide a radial view of the plasma and have a sightline intersection with the South–South neutral beam line. From this radial geometry, active CX signals from mostly trapped particles are observed. Each channel has a spatial resolution of 3–4 cm, a temporal resolution of 200 kHz, and an average pitch angle resolution of a few degrees. The three-sensor configuration allows for coarse energy resolution of the CX signals; each sensor sees similar sightlines but different filter thicknesses alter the energy cutoffs by known amounts. Experimental data show that all channels are collecting data as intended. The signal to noise ratio is typically around 15. Preliminary data analysis shows a correlation between the SSNPA signal and magnetohydrodynamic activity in the plasma as expected.

Instruments & Instrumentation↗
DOE OSTI2024

Mining experimental magnetized liner inertial fusion data: Trends in stagnation morphology

In magnetized liner inertial fusion (MagLIF), a cylindrical liner filled with fusion fuel is imploded with the goal of producing a one-dimensional plasma column at thermonuclear conditions. However, structures attributed to three-dimensional effects are observed in self-emission x-ray images. Despite this, the impact of many experimental inputs on the column morphology has not been characterized. We demonstrate the use of a linear regression analysis to explore correlations between morphology and a wide variety of experimental inputs across 57 MagLIF experiments. Results indicate the possibility of several unexplored effects. For example, we demonstrate that increasing the initial magnetic field correlates with improved stability. Although intuitively expected, this has never been quantitatively assessed in integrated MagLIF experiments. We also demonstrate that azimuthal drive asymmetries resulting from the geometry of the “current return can” appear to measurably impact the morphology. In conjunction with several counterintuitive null results, we expect the observed correlations will encourage further experimental, theoretical, and simulation-based studies. Finally, we note that the method used in this work is general and may be applied to explore not only correlations between input conditions and morphology but also with other experimentally measured quantities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Topology and spectral entanglement in cavity-mediated photon scattering

Here, we develop a microscopic diagrammatic theory for cavity-mediated photon scattering in a topological one-dimensional insulator described by the Su–Schrieffer–Heeger model. Within the velocity-gauge formulation, we derive the photon self-energy and vertex corrections arising from virtual electron–hole excitations coupled to a quantized cavity mode, and we evaluate the resulting polariton dispersion and two-photon correlation spectra. Our analysis shows that vacuum fluctuations of the cavity field induce a momentum-resolved self-energy that mixes conduction and valence bands through virtual photon exchange, producing interband hybridization and avoided crossings in the electronic dispersion. This “cavity dressing” is symmetry-dependent, vanishing at the Brillouin-zone edge where the dipole matrix element is zero, and its strength is controlled by the spatial coherence range ζ ≈ (l c /a) 2 of virtual excitations. We further examine how the cavity modifies nonlinear optical observables, including the Kerr nonlinearity and biphoton spectral entanglement, and identify the regimes where these effects become sensitive to the underlying topological phase. The theoretical framework established here provides a unified description of light–matter coupling in topological and polaritonic systems, bridging solid-state cavity QED with the emerging field of cavity-modified quantum materials. Our results suggest that engineered photonic environments can coherently reshape the electronic landscape of topological insulators, offering new routes to control collective electronic and optical phenomena through vacuum-field fluctuations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Mechanical design of the highly porous cuttlebone: A bioceramic hard buoyancy tank for cuttlefish

Cuttlefish, a unique group of marine mollusks, produces an internal biomineralized shell, known as cuttlebone, which is an ultra-lightweight cellular structure (porosity, ~93 vol%) used as the animal’s hard buoyancy tank. Although cuttlebone is primarily composed of a brittle mineral, aragonite, the structure is highly damage tolerant and can withstand water pressure of about 20 atmospheres (atm) for the species Sepia officinalis . Currently, our knowledge on the structural origins for cuttlebone’s remarkable mechanical performance is limited. Combining quantitative three-dimensional (3D) structural characterization, four-dimensional (4D) mechanical analysis, digital image correlation, and parametric simulations, here we reveal that the characteristic chambered “wall–septa” microstructure of cuttlebone, drastically distinct from other natural or engineering cellular solids, allows for simultaneous high specific stiffness (8.4 MN∙m/kg) and energy absorption (4.4 kJ/kg) upon loading. We demonstrate that the vertical walls in the chambered cuttlebone microstructure have evolved an optimal waviness gradient, which leads to compression-dominant deformation and asymmetric wall fracture, accomplishing both high stiffness and high energy absorption. Moreover, the distribution of walls is found to reduce stress concentrations within the horizontal septa, facilitating a larger chamber crushing stress and a more significant densification. Overall, the design strategies revealed here can provide important lessons for the development of low-density, stiff, and damage-tolerant cellular ceramics.

42 ENGINEERING↗
DOE OSTI2011

The radio/gamma-ray connection in active galactic nuclei in the era of the Fermi Large Area Telescope

We present a detailed statistical analysis of the correlation between radio and gamma-ray emission of the active galactic nuclei (AGNs) detected by Fermi during its first year of operation, with the largest data sets ever used for this purpose. We use both archival interferometric 8.4 GHz data (from the Very Large Array and ATCA, for the full sample of 599 sources) and concurrent single-dish 15 GHz measurements from the Owens Valley Radio Observatory (OVRO, for a sub sample of 199 objects). Our unprecedentedly large sample permits us to assess with high accuracy the statistical significance of the correlation, using a surrogate data method designed to simultaneously account for common-distance bias and the effect of a limited dynamical range in the observed quantities. We find that the statistical significance of a positive correlation between the centimeter radio and the broadband (E > 100 MeV) gamma-ray energy flux is very high for the whole AGN sample, with a probability of <10 –7 for the correlation appearing by chance. Using the OVRO data, we find that concurrent data improve the significance of the correlation from 1.6 × 10 –6 to 9.0 × 10 –8 . Our large sample size allows us to study the dependence of correlation strength and significance on specific source types and gamma-ray energy band. As a result, we find that the correlation is very significant (chance probability < 10 –7 ) for both flat spectrum radio quasars and BL Lac objects separately; a dependence of the correlation strength on the considered gamma-ray energy band is also present, but additional data will be necessary to constrain its significance.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

THE RADIO/GAMMA-RAY CONNECTION IN ACTIVE GALACTIC NUCLEI IN THE ERA OF THE FERMI LARGE AREA TELESCOPE

We present a detailed statistical analysis of the correlation between radio and gamma-ray emission of the active galactic nuclei (AGNs) detected by Fermi during its first year of operation, with the largest data sets ever used for this purpose. We use both archival interferometric 8.4 GHz data (from the Very Large Array and ATCA, for the full sample of 599 sources) and concurrent single-dish 15 GHz measurements from the Owens Valley Radio Observatory (OVRO, for a sub sample of 199 objects). Our unprecedentedly large sample permits us to assess with high accuracy the statistical significance of the correlation, using a surrogate data method designed to simultaneously account for common-distance bias and the effect of a limited dynamical range in the observed quantities. We find that the statistical significance of a positive correlation between the centimeter radio and the broadband (E > 100 MeV) gamma-ray energy flux is very high for the whole AGN sample, with a probability of <10{sup -7} for the correlation appearing by chance. Using the OVRO data, we find that concurrent data improve the significance of the correlation from 1.6 x 10{sup -6} to 9.0 x 10{sup -8}. Our large sample size allows us to study the dependence of correlation strength and significance on specific source types and gamma-ray energy band. We find that the correlation is very significant (chance probability < 10{sup -7}) for both flat spectrum radio quasars and BL Lac objects separately; a dependence of the correlation strength on the considered gamma-ray energy band is also present, but additional data will be necessary to constrain its significance.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Comparing simulated emission from molecular clouds using experimental design

We propose a new approach to comparing simulated observations that enables us to determine the significance of the underlying physical effects. We utilize the methodology of experimental design, a subfield of statistical analysis, to establish a framework for comparing simulated position-position-velocity data cubes to each other. We propose three similarity metrics based on methods described in the literature: principal component analysis, the spectral correlation function, and the Cramer multi-variate two-sample similarity statistic. Using these metrics, we intercompare a suite of mock observational data of molecular clouds generated from magnetohydrodynamic simulations with varying physical conditions. Using this framework, we show that all three metrics are sensitive to changing Mach number and temperature in the simulation sets, but cannot detect changes in magnetic field strength and initial velocity spectrum. We highlight the shortcomings of one-factor-at-a-time designs commonly used in astrophysics and propose fractional factorial designs as a means to rigorously examine the effects of changing physical properties while minimizing the investment of computational resources.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Temperature-driven growth of antiferromagnetic domains in thin-film FeRh

The evolution of the antiferromagnetic phase across the temperature-driven ferromagnetic (FM) to antiferromagnetic (AF) phase transition in epitaxial FeRh thin films was studied by x-ray magnetic linear and circular dichroism (XMLD and XMCD) and photoemission electron microscopy. By comparing XMLD and XMCD images recorded at the same temperature, the AF phase was identified, its structure directly imaged, and its evolution studied across the transition. A quantitative analysis of the correlation length of the images shows differences between the characteristic length scale of the two phases with the AF phase having a finer feature size. The asymmetry of the transition from FM to AF upon cooling and AF-FM upon heating is evidenced: upon cooling the formation of AF phase is dominated by nucleation at defects, with little subsequent growth, resulting in a small and non-random final AF domain structure, while upon heating, heterogeneous nucleation at different sites followed by significant domain size growth of the FM phase is observed, resulting in a non-reproducible final FM large domain structure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Characteristics of infrasound signals from North Korean underground nuclear explosions on 2016 January 6 and September 9

This study investigates infrasound signals from North Korean underground nuclear explosions on 2016 January 6 at 01:30:01 UTC (UNEJ16) and on 2016 September 9 at 00:30:01 UTC (UNES16). Infrasound observations from eight, seismo-acoustic arrays at distances of 300–600 km from the explosions and aligned approximately east–west across South Korea are used in the analysis. Progressive multichannel correlation was used followed by analyst review to identify the infrasound arrivals with estimates of azimuth, phase velocity, celerity, F-statistic, signal-to-noise ratio, correlation and signal duration. These observations are compared with model predictions of travel path, ray turning height and relative amplitudes, using Ground-to-Space atmospheric specifications at the time of each explosion. Stratospheric conditions at the time of the UNEJ16 are favourable to the Korean arrays when all arrays detected infrasound signals, while the few detections for UNES16 are consistent with unfavourable stratospheric winds and reduced amplitude thermospheric arrivals. Infrasound locations of the explosions were estimated using the Bayesian Infrasonic Source Location method (BISL), with wind-corrected backazimuth estimates based on the ray tracing calculations. Adding wind corrections to BISL resulted in infrasound locations up to 50 per cent closer to the seismic epicentre over uncorrected locations. The UNES16 infrasound location estimate using atmospheric corrections had larger differences (~67 km difference from seismic location) than the UNEJ16 location (~11 km difference from seismic location) as a result of the reduced number of observations and poor azimuth coverage for UNES16. Based on stratospheric phases recorded at six arrays and the empirical yield-scaling relations, wind corrected amplitudes from UNEJ16 lead to equivalent infrasound source energy estimates between 0.9 and 16.1 tons of TNT detonated in the atmosphere with an average of 6.4 tons and a standard deviation of 4.6 tons.

58 GEOSCIENCES↗
DOE OSTI2025

Unveiling shared genetic regulators of plant architectural and biomass yield traits in the Sorghum Association Panel

Abstract Sorghum is emerging as an ideal genetic model for designing high-biomass bioenergy crops. Biomass yield, a complex trait influenced by various plant architectural characteristics, is typically regulated by numerous genes. This study aimed to dissect the genetic regulators underlying 14 plant architectural traits and 10 biomass yield traits in the Sorghum Association Panel across two growing seasons. We identified 321 associated loci through genome-wide association studies (GWAS), involving 234 264 single nucleotide polymorphisms (SNPs). These loci include genes with known associations to biomass traits, such as maturity, dwarfing (Dw), and leafbladeless1, as well as several uncharacterized loci not previously linked to these traits. We also identified 22 pleiotropic loci associated with variation in multiple phenotypes. Three of these loci, located on chromosomes 3 (S03_15463061), 6 (S06_42790178; Dw2), and 9 (S09_57005346; Dw1), exerted significant and consistent effects on multiple traits across both growing seasons. Additionally, we identified three genomic hotspots on chromosomes 6, 7, and 9, each containing multiple SNPs associated with variation in plant architecture and biomass yield traits. Chromosome-wise correlation analyses revealed multiple blocks of positively associated SNPs located near or within the same genomic regions. Finally, genome-wide correlation-based network analysis showed that loci associated with flowering, plant height, leaf traits, plant density, and tiller number per plant were highly interconnected with other genetic loci influencing plant architectural and biomass yield traits. The pyramiding of favorable alleles related to these traits holds promise for enhancing the future development of bioenergy sorghum crops.

Singh, Anuradha (ORCID:0000000197149095)↗
DOE OSTI2022

Nuclear properties with semilocal momentum-space regularized chiral interactions beyond N 2 LO

We present a comprehensive investigation of few-nucleon systems as well as light and medium-mass nuclei up to A = 48 using the current Low Energy Nuclear Physics International Collaboration two-nucleon interactions in combination with the third-order (N 2 LO) three-nucleon forces. To address the systematic overbinding of nuclei starting from A ≈ 10 found in our earlier study utilizing the N 2 LO two- and three-nucleon forces, we take into account higher-order corrections to the two-nucleon potentials up through fifth order in chiral effective field theory. Here, the resulting Hamiltonian can be completely determined using the A = 3 binding energies and selected nucleon-deuteron cross sections as input. It is then shown to predict other nucleon-deuteron scattering observables and spectra of light p-shell nuclei, for which a detailed correlated truncation error analysis is performed, in agreement with experimental data. Moreover, the predicted ground state energies of nuclei in the oxygen isotopic chain from 14 O to 26 O as well as 40 Ca and 48 Ca show a remarkably good agreement with experimental values, given that the Hamiltonian is fixed completely from the A ≤ 3 data, once the fourth-order (N 3 LO) corrections to the two-nucleon interactions are taken into account. On the other hand, the charge radii are found to be underpredicted by ≈10% for the oxygen isotopes and by almost 20% for 40 Ca and 48 Ca.

20 ≤ A ≤ 38↗
DOE OSTI2013

Energy Density Functional for Nuclei and Neutron Stars

Background: Recent observational data on neutron star masses and radii provide stringent constraints on the equation of state of neutron rich matter [ Annu. Rev. Nucl. Part. Sci. 62 485 (2012)]. Purpose: We aim to develop a nuclear energy density functional that can be simultaneously applied to finite nuclei and neutron stars. Methods: We use the self-consistent nuclear density functional theory (DFT) with Skyrme energy density functionals and covariance analysis to assess correlations between observables for finite nuclei and neutron stars. In a first step two energy functionals a high density energy functional giving reasonable neutron properties, and a low density functional fitted to nuclear properties are matched. In a second step, we optimize a new functional using exactly the same protocol as in earlier studies pertaining to nuclei but now including neutron star data. This allows direct comparisons of performance of the new functional relative to the standard one. Results: The new functional TOV-min yields results for nuclear bulk properties (energy, rms radius, diffraction radius, and surface thickness) that are of the same quality as those obtained with the established Skyrme functionals, including SV-min. When comparing SV-min and TOV-min, isoscalar nuclear matter indicators vary slightly while isovector properties are changed considerably. We discuss neutron skins, dipole polarizability, separation energies of the heaviest elements, and proton and neutron drip lines. We confirm a correlation between the neutron skin of 208Pb and the neutron star radius. Conclusions: We demonstrate that standard energy density functionals optimized to nuclear data do not carry information on the expected maximum neutron star mass, and that predictions can only be made within an extremely broad uncertainty band. For atomic nuclei, the new functional TOV-min performs at least as well as the standard nuclear functionals, but it also reproduces expected neutron star data within assumed error bands. This functional is expected to yield more reliable predictions in the region of very neutron rich heavy nuclei.

Erler, J.↗
DOE OSTI2024

Testing hadronic-model predictions of depth of maximum of air-shower profiles and ground-particle signals using hybrid data of the Pierre Auger Observatory

We test the predictions of hadronic interaction models regarding the depth of maximum of air-shower profiles, X max , and ground-particle signals in water-Cherenkov detectors at 1000 m from the shower core, S ( 1000 ) , using the data from the fluorescence and surface detectors of the Pierre Auger Observatory. The test consists of fitting the measured two-dimensional ( S ( 1000 ) , X max ) distributions using templates for simulated air showers produced with hadronic interaction models pos-, et--04, 2.3d and leaving the scales of predicted X max and the signals from hadronic component at ground as free-fit parameters. The method relies on the assumption that the mass composition remains the same at all zenith angles, while the longitudinal shower development and attenuation of ground signal depend on the mass composition in a correlated way. The analysis was applied to 2239 events detected by both the fluorescence and surface detectors of the Pierre Auger Observatory with energies between 10 18.5 eV to 10 19.0 eV and zenith angles below 60°. We found, that within the assumptions of the method, the best description of the data is achieved if the predictions of the hadronic interaction models are shifted to deeper X max values and larger hadronic signals at all zenith angles. Given the magnitude of the shifts and the data sample size, the statistical significance of the improvement of data description using the modifications considered in the paper is larger than 5 σ even for any linear combination of experimental systematic uncertainties. Published by the American Physical Society 2024

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

Revisiting no-scale supergravity inspired scenarios: Updated theoretical and phenomenological constraints

We consider no-scale supergravity inspired scenarios, where the gravitino mass and related soft supersymmetry-breaking parameters are determined dynamically by radiative corrections to an essentially flat tree-level potential in the supersymmetry-breaking hidden sector. We examine the theoretical and phenomenological viability of such a mechanism, when including up-to-date calculations of the low-energy sparticle spectrum and taking into account the latest LHC results and other experimental constraints. We (re)emphasize the role of the scale-dependent vacuum energy contribution to the effective potential, in obtaining realistic no-scale electroweak minima, examining carefully the impact of boundary conditions and of variants of the minimization procedure. We also discuss and implement the B{sub 0} (soft breaking Higgs mixing parameter) input boundary condition at high scale, therefore fixing tan{beta}(B{sub 0}) at low scales. For general high-scale boundary conditions with B{sub 0},m{sub 0},{center_dot}{center_dot}{center_dot}{ne}0, our analysis provides theoretical correlations among the supersymmetric, soft, and vacuum energy parameters and related phenomenological consequences at the LHC. For instance, a zero vacuum energy at the grand unified theory scale would lead to a decoupled supersymmetric spectrum, together with a light standard-model-like Higgs boson at the electroweak scale. Given the experimental exclusion limits, a substantial class of the boundary conditions, and, in particular, the strict no-scale with m{sub 0}=A{sub 0}=B{sub 0}=0, are only compatible with a stau being the lightest Minimal Supersymmetric Standard Model particle. Then, an enlarged allowed parameter space emerges when assuming a gravitino lightest supersymmetric particle to account for the observed dark matter relic density.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Biexciton Condensation in Electron-Hole-Doped Hubbard Bilayers: A Sign-Problem-Free Quantum Monte Carlo Study

The bilayer Hubbard model with electron-hole doping is an ideal platform to study excitonic orders due to suppressed recombination via spatial separation of electrons and holes. However, suffering from the sign problem, previous quantum Monte Carlo studies could not arrive at an unequivocal conclusion regarding the presence of phases with clear signatures of excitonic condensation in bilayer Hubbard models. Here, we develop a determinant quantum Monte Carlo algorithm for the bilayer Hubbard model that is sign-problem-free for equal and opposite doping in the two layers and study excitonic order and charge and spin density modulations as a function of chemical potential difference between the two layers, on-site Coulomb repulsion, and interlayer interaction. In the intermediate coupling regime and in proximity to the SU(4)-symmetric point, we find a biexcitonic condensate phase at finite electron-hole doping, as well as a competing (π, π) charge density wave state. In this study, we extract the Berezinskii-Kosterlitz-Thouless transition temperature from superfluid density and a finite-size scaling analysis of the correlation functions and explain our results in terms of an effective biexcitonic hard-core boson model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗