Showing 1–50 of 170 results for author: Alavi, A

Guidelines for Augmentation Selection in Contrastive Learning for Time Series Classification

Authors: Ziyu Liu, Azadeh Alavi, Minyi Li, Xiang Zhang

Abstract: Self-supervised contrastive learning has become a key technique in deep learning, particularly in time series analysis, due to its ability to learn meaningful representations without explicit supervision. Augmentation is a critical component in contrastive learning, where different augmentations can dramatically impact performance, sometimes influencing accuracy by over 30%. However, the selection… ▽ More Self-supervised contrastive learning has become a key technique in deep learning, particularly in time series analysis, due to its ability to learn meaningful representations without explicit supervision. Augmentation is a critical component in contrastive learning, where different augmentations can dramatically impact performance, sometimes influencing accuracy by over 30%. However, the selection of augmentations is predominantly empirical which can be suboptimal, or grid searching that is time-consuming. In this paper, we establish a principled framework for selecting augmentations based on dataset characteristics such as trend and seasonality. Specifically, we construct 12 synthetic datasets incorporating trend, seasonality, and integration weights. We then evaluate the effectiveness of 8 different augmentations across these synthetic datasets, thereby inducing generalizable associations between time series characteristics and augmentation efficiency. Additionally, we evaluated the induced associations across 6 real-world datasets encompassing domains such as activity recognition, disease diagnosis, traffic monitoring, electricity usage, mechanical fault prognosis, and finance. These real-world datasets are diverse, covering a range from 1 to 12 channels, 2 to 10 classes, sequence lengths of 14 to 1280, and data frequencies from 250 Hz to daily intervals. The experimental results show that our proposed trend-seasonality-based augmentation recommendation algorithm can accurately identify the effective augmentations for a given time series dataset, achieving an average Recall@3 of 0.667, outperforming baselines. Our work provides guidance for studies employing contrastive learning in time series analysis, with wide-ranging applications. All the code, datasets, and analysis results will be released at https://github.com/DL4mHealth/TS-Contrastive-Augmentation-Recommendation. △ Less

Submitted 12 July, 2024; originally announced July 2024.

Comments: 20 pages, 11 figures

Leveraging SPD Matrices on Riemannian Manifolds in Quantum Classical Hybrid Models for Structural Health Monitoring

Authors: Azadeh Alavi, Sanduni Jayasinghe

Abstract: Realtime finite element modeling of bridges assists modern structural health monitoring systems by providing comprehensive insights into structural integrity. This capability is essential for ensuring the safe operation of bridges and preventing sudden catastrophic failures. However, FEM computational cost and the need for realtime analysis pose significant challenges. Additionally, the input data… ▽ More Realtime finite element modeling of bridges assists modern structural health monitoring systems by providing comprehensive insights into structural integrity. This capability is essential for ensuring the safe operation of bridges and preventing sudden catastrophic failures. However, FEM computational cost and the need for realtime analysis pose significant challenges. Additionally, the input data is a 7 dimensional vector, while the output is a 1017 dimensional vector, making accurate and efficient analysis particularly difficult. In this study, we propose a novel hybrid quantum classical Multilayer Perceptron pipeline leveraging Symmetric Positive Definite matrices and Riemannian manifolds for effective data representation. To maintain the integrity of the qubit structure, we utilize SPD matrices, ensuring data representation is well aligned with the quantum computational framework. Additionally, the method leverages polynomial feature expansion to capture nonlinear relationships within the data. The proposed pipeline combines classical fully connected neural network layers with quantum circuit layers to enhance model performance and efficiency. Our experiments focused on various configurations of such hybrid models to identify the optimal structure for accurate and efficient realtime analysis. The best performing model achieved a Mean Squared Error of 0.00031, significantly outperforming traditional methods. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: 3 pages, 1 figure

Aperiodic fragments in periodic solids: Eliminating the need for supercells and background charges in electronic structure calculations of defects

Authors: Robert H. Lavroff, Daniel Kats, Lorenzo Maschio, Nikolay Bogdanov, Ali Alavi, Anastassia N. Alexandrova, Denis Usvyat

Abstract: To date, computational methods for modeling defects (vacancies, adsorbates, etc.) rely on periodic supercells in which the defect is far enough from its repeated image such that they can be assumed non-interacting. Defects in real solids, however, can be spaced microns or more apart, whereas affordable supercells for density functional theory calculations are no wider than a few nanometers. The re… ▽ More To date, computational methods for modeling defects (vacancies, adsorbates, etc.) rely on periodic supercells in which the defect is far enough from its repeated image such that they can be assumed non-interacting. Defects in real solids, however, can be spaced microns or more apart, whereas affordable supercells for density functional theory calculations are no wider than a few nanometers. The relative proximity and periodic repetition of the defect's images may lead to unphysical artifacts, especially if the defect is charged and/or open-shell. Furthermore, to avoid divergence of the periodic electrostatics, a compensating background charge must be introduced if the defect is charged. Even if post-hoc corrections are used, this is a source of unquantifiable error and, in some cases, renders energies useless. In this communication, we introduce an embedding formalism such that a pristine, primitive unit cell may be used for the periodic mean field, after which atoms may be moved or charged within an embedded fragment. This fragment can then be treated with a post-Hartree Fock method to capture important electron correlations pertaining to the defect. By eliminating the need for compensating background charges and periodicity of the defect, we circumvent all associated unphysicalities and numerical issues. Furthermore, the primitive cell calculations drastically reduce computational expense compared to supercell approaches. This method is size-intensive with respect to energy differences and can be routinely applied even to multireference defects, localized excited states, etc. using a variety of fragment solvers. In examining with this approach bond-breaking in a fluorine-substituted graphane monolayer, a difficult testing ground for condensed-phase electronic structure methods, we observe key aspects of the dissociation pathway, specifically a covalent-to-ionic avoided crossing. △ Less

Submitted 5 June, 2024; originally announced June 2024.

Comments: This manuscript has not yet been submitted to a journal

Leveraging Knowlegde Graphs for Interpretable Feature Generation

Authors: Mohamed Bouadi, Arta Alavi, Salima Benbernou, Mourad Ouziri

Abstract: The quality of Machine Learning (ML) models strongly depends on the input data, as such Feature Engineering (FE) is often required in ML. In addition, with the proliferation of ML-powered systems, especially in critical contexts, the need for interpretability and explainability becomes increasingly important. Since manual FE is time-consuming and requires case specific knowledge, we propose KRAFT,… ▽ More The quality of Machine Learning (ML) models strongly depends on the input data, as such Feature Engineering (FE) is often required in ML. In addition, with the proliferation of ML-powered systems, especially in critical contexts, the need for interpretability and explainability becomes increasingly important. Since manual FE is time-consuming and requires case specific knowledge, we propose KRAFT, an AutoFE framework that leverages a knowledge graph to guide the generation of interpretable features. Our hybrid AI approach combines a neural generator to transform raw features through a series of transformations and a knowledge-based reasoner to evaluate features interpretability using Description Logics (DL). The generator is trained through Deep Reinforcement Learning (DRL) to maximize the prediction accuracy and the interpretability of the generated features. Extensive experiments on real datasets demonstrate that KRAFT significantly improves accuracy while ensuring a high level of interpretability. △ Less

Submitted 1 June, 2024; originally announced June 2024.

Galaxy Rest-Frame UV Colors at z ~ 2-4 with HST UVCANDELS

Authors: Alexa Morales, Steven Finkelstein, Micaela Bagley, Anahita Alavi, Norman Grogin, Nimish Hathi, Anton Koekemoer, Kalina Nedkova, Laura Prichard, Marc Rafelski, Ben Sunnquist, Sina Taamoli, Harry Teplitz, Xin Wang, Rogier Windhorst, L. Y. Aaron Yung

Abstract: We present an analysis of rest-frame UV colors of 17,243 galaxies at $z\sim2-4$ in the HST UVCANDELS fields: GOODS-N, GOODS-S, COSMOS, and EGS. Here, we study the rest-frame UV spectral slope, $β$, measured via model spectra obtained via spectral energy distribution (SED) fitting, $β_{SED}$, and explore its correlation with various galaxy parameters (photometric redshift, UV magnitude, stellar mas… ▽ More We present an analysis of rest-frame UV colors of 17,243 galaxies at $z\sim2-4$ in the HST UVCANDELS fields: GOODS-N, GOODS-S, COSMOS, and EGS. Here, we study the rest-frame UV spectral slope, $β$, measured via model spectra obtained via spectral energy distribution (SED) fitting, $β_{SED}$, and explore its correlation with various galaxy parameters (photometric redshift, UV magnitude, stellar mass, dust attenuation, star formation rate [SFR], and specific SFR) obtained via SED fitting with Dense Basis. We also obtain measurements for $β$ via photometric power-law fitting and compare them to our SED-fit-based results, finding good agreement on average. While we find little evolution in $β$ with redshift from $z=2-4$ for the full population, there are clear correlations between $β$ (and related parameters) when binned by stellar mass. For this sample, lower stellar mass galaxies (log[$M_*$] = 7.5-8.5 $M_\odot$) are typically bluer ($β_{SED}=-2.0\pm 0.2$ / $β_{PL} = -2.1\pm0.4$), fainter ($MUV = -17.8^{+0.7}_{-0.6}$) less dusty ($A{v}=0.4\pm0.1$ mag), exhibit lower rates of star formation (log[SFR]=$0.1\pm0.2 M_\odot/$ yr) and higher specific star formation rates (log[sSFR]=$-8.2\pm0.2 \ \mathrm{yr}^{-1}$) than their high-mass counterparts. Higher-mass galaxies (log[$M_*$] $=10.0-12.0 \ M_\odot$) are on average redder ($β_{SED}=-0.9^{+0.8}_{-0.5}$ / $β_{PL}=-1.0^{+0.8}_{-0.5}$), brighter ($MUV=-19.6^{+1.0}_{-1.2}$), dustier ($Av = 0.9^{+0.5}_{-0.4}$ mag), have higher SFRs (log[SFR]=$1.2^{+0.6}_{-1.1} M_\odot$ yr), and lower sSFRs (log[sSFR]=$-9.1^{+0.5}_{-1.1} {yr}^{-1}$). This study's substantial sample size provides a benchmark for demonstrating that the rest-frame UV spectral slope correlates with stellar mass-dependent galaxy characteristics at $z\sim2-4$, a relationship less discernible with smaller datasets typically available at higher redshifts. △ Less

Submitted 31 May, 2024; originally announced May 2024.

Comments: 16 pages, 5 tables, 6 figures, submitted to ApJ

Euclid. I. Overview of the Euclid mission

Authors: Euclid Collaboration, Y. Mellier, Abdurro'uf, J. A. Acevedo Barroso, A. Achúcarro, J. Adamek, R. Adam, G. E. Addison, N. Aghanim, M. Aguena, V. Ajani, Y. Akrami, A. Al-Bahlawan, A. Alavi, I. S. Albuquerque, G. Alestas, G. Alguero, A. Allaoui, S. W. Allen, V. Allevato, A. V. Alonso-Tetilla, B. Altieri, A. Alvarez-Candal, A. Amara, L. Amendola , et al. (1086 additional authors not shown)

Abstract: The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14… ▽ More The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue`Euclid on Sky'

UVCANDELS: The role of dust on the stellar mass-size relation of disk galaxies at 0.5 $\leq z \leq$ 3.0

Authors: Kalina V. Nedkova, Marc Rafelski, Harry I. Teplitz, Vihang Mehta, Laura DeGroot, Swara Ravindranath, Anahita Alavi, Alexander Beckett, Norman A. Grogin, Boris Häußler, Anton M. Koekemoer, Grecco A. Oyarzún, Laura Prichard, Mitchell Revalski, Gregory F. Snyder, Ben Sunnquist, Xin Wang, Rogier A. Windhorst, Nima Chartab, Christopher J. Conselice, Yicheng Guo, Nimish Hathi, Matthew J. Hayes, Zhiyuan Ji, Keunho J. Kim , et al. (8 additional authors not shown)

Abstract: We use the Ultraviolet Imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey fields (UVCANDELS) to measure half-light radii in the rest-frame far-UV for $\sim$16,000 disk-like galaxies over $0.5\leq z \leq 3$. We compare these results to rest-frame optical sizes that we measure in a self-consistent way and find that the stellar mass-size relation of disk galaxies is steeper… ▽ More We use the Ultraviolet Imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey fields (UVCANDELS) to measure half-light radii in the rest-frame far-UV for $\sim$16,000 disk-like galaxies over $0.5\leq z \leq 3$. We compare these results to rest-frame optical sizes that we measure in a self-consistent way and find that the stellar mass-size relation of disk galaxies is steeper in the rest-frame UV than in the optical across our entire redshift range. We show that this is mainly driven by massive galaxies ($\gtrsim10^{10}$M$_\odot$), which we find to also be among the most dusty. Our results are consistent with the literature and have commonly been interpreted as evidence of inside-out growth wherein galaxies form their central structures first. However, they could also suggest that the centers of massive galaxies are more heavily attenuated than their outskirts. We distinguish between these scenarios by modeling and selecting galaxies at $z=2$ from the VELA simulation suite in a way that is consistent with UVCANDELS. We show that the effects of dust alone can account for the size differences we measure at $z=2$. This indicates that, at different wavelengths, size differences and the different slopes of the stellar mass-size relation do not constitute evidence for inside-out growth. △ Less

Submitted 28 June, 2024; v1 submitted 17 May, 2024; originally announced May 2024.

Comments: Accepted for publication in ApJ. 22 pages, 12 figures, and 4 tables

Orbital optimisation in xTC transcorrelated methods

Authors: Daniel Kats, Evelin M. C. Christlmaier, Thomas Schraivogel, Ali Alavi

Abstract: We present a combination of the bi-orthogonal orbital optimisation framework with the recently introduced xTC version of transcorrelation. This allows us to implement non-iterative perturbation based methods on top of the transcorrelated Hamiltonian. Besides, the orbital optimisation influences results of other truncated methods, such as the distinguishable cluster with singles and doubles. The ac… ▽ More We present a combination of the bi-orthogonal orbital optimisation framework with the recently introduced xTC version of transcorrelation. This allows us to implement non-iterative perturbation based methods on top of the transcorrelated Hamiltonian. Besides, the orbital optimisation influences results of other truncated methods, such as the distinguishable cluster with singles and doubles. The accuracy of these methods in comparison to standard xTC methods is demonstrated, and the advantages and disadvantages of the orbital optimisation are discussed. △ Less

Submitted 11 April, 2024; originally announced April 2024.

Journal ref: Faraday Discussions, 2024

WFC3 Infrared Spectroscopic Parallel (WISP) Survey: Photometric and Emission Line Data Release

Authors: A. J. Battisti, M. B. Bagley, M. Rafelski, I. Baronchelli, Y. S. Dai, A. L. Henry, H. Atek, J. Colbert, M. A. Malkan, P. J. McCarthy, C. Scarlata, B. Siana, H. I. Teplitz, A. Alavi, K. Boyett, A. J. Bunker, J. P. Gardner, N. P. Hathi, D. Masters, V. Mehta, M. Rutkowski, K. Shahinyan, B. Sunnquist, X. Wang

Abstract: We present reduced images and catalogues of photometric and emission line data ($\sim$230,000 and $\sim$8,000 sources, respectively) for the WFC3 Infrared Spectroscopic Parallel (WISP) Survey. These data are made publicly available on the Mikulski Archive for Space Telescopes (MAST) and include reduced images from various facilities: ground-based $ugri$, HST WFC3, and Spitzer IRAC (Infrared Array… ▽ More We present reduced images and catalogues of photometric and emission line data ($\sim$230,000 and $\sim$8,000 sources, respectively) for the WFC3 Infrared Spectroscopic Parallel (WISP) Survey. These data are made publicly available on the Mikulski Archive for Space Telescopes (MAST) and include reduced images from various facilities: ground-based $ugri$, HST WFC3, and Spitzer IRAC (Infrared Array Camera). Coverage in at least one additional filter beyond the WFC3/IR data are available for roughly half of the fields (227 out of 483), with $\sim$20% (86) having coverage in six or more filters from $u$-band to IRAC 3.6$μ$m (0.35-3.6$μ$m). For the lower spatial resolution (and shallower) ground-based and IRAC data, we perform PSF-matched, prior-based, deconfusion photometry (i.e., forced-photometry) using the TPHOT software to optimally extract measurements or upper limits. We present the methodology and software used for the WISP emission line detection and visual inspection. The former adopts a continuous wavelet transformation that significantly reduces the number of spurious sources as candidates before the visual inspection stage. We combine both WISP catalogues and perform SED fitting on galaxies with reliable spectroscopic redshifts and multi-band photometry to measure their stellar masses. We stack WISP spectra as functions of stellar mass and redshift and measure average emission line fluxes and ratios. We find that WISP emission line sources are typically `normal' star-forming galaxies based on the Mass-Excitation diagram ([OIII]/H$β$ vs. $M_\star$; $0.74<z_\mathrm{grism}<2.31$), the galaxy main sequence (SFR vs. $M_\star$; $0.30<z_\mathrm{grism}<1.45$), $S_{32}$ ratio vs. $M_\star$ ($0.30<z_\mathrm{grism}<0.73$), and $O_{32}$ and $R_{23}$ ratios vs. $M_\star$ ($1.27<z_\mathrm{grism}<1.45$). △ Less

Submitted 6 April, 2024; originally announced April 2024.

Comments: 36 pages, 21 figures, 17 tables. Accepted for publication in MNRAS. The WISP Photometric and Emission Line catalogues and reduced images are in the process of being added as HLSPs to the WISP MAST website (https://archive.stsci.edu/prepds/wisp/). Please email the first-author (provided in paper) to request access to files prior to the MAST release

Self-Supervised Learning for Time Series: Contrastive or Generative?

Authors: Ziyu Liu, Azadeh Alavi, Minyi Li, Xiang Zhang

Abstract: Self-supervised learning (SSL) has recently emerged as a powerful approach to learning representations from large-scale unlabeled data, showing promising results in time series analysis. The self-supervised representation learning can be categorized into two mainstream: contrastive and generative. In this paper, we will present a comprehensive comparative study between contrastive and generative m… ▽ More Self-supervised learning (SSL) has recently emerged as a powerful approach to learning representations from large-scale unlabeled data, showing promising results in time series analysis. The self-supervised representation learning can be categorized into two mainstream: contrastive and generative. In this paper, we will present a comprehensive comparative study between contrastive and generative methods in time series. We first introduce the basic frameworks for contrastive and generative SSL, respectively, and discuss how to obtain the supervision signal that guides the model optimization. We then implement classical algorithms (SimCLR vs. MAE) for each type and conduct a comparative analysis in fair settings. Our results provide insights into the strengths and weaknesses of each approach and offer practical recommendations for choosing suitable SSL methods. We also discuss the implications of our findings for the broader field of representation learning and propose future research directions. All the code and data are released at \url{https://github.com/DL4mHealth/SSL_Comparison}. △ Less

Submitted 14 March, 2024; originally announced March 2024.

Comments: Published at the AI4TS Workshop, IJCAI 2023

X marks the spot: accurate energies from intersecting extrapolations of continuum quantum Monte Carlo data

Authors: Seyed Mohammadreza Hosseini, Ali Alavi, Pablo Lopez Rios

Abstract: We explore the application of an extrapolative method that yields very accurate total and relative energies from variational and diffusion quantum Monte Carlo (VMC and DMC) results. For a trial wave function consisting of a small configuration interaction (CI) wave function obtained from full CI quantum Monte Carlo and reoptimized in the presence of a Jastrow factor and an optional backflow transf… ▽ More We explore the application of an extrapolative method that yields very accurate total and relative energies from variational and diffusion quantum Monte Carlo (VMC and DMC) results. For a trial wave function consisting of a small configuration interaction (CI) wave function obtained from full CI quantum Monte Carlo and reoptimized in the presence of a Jastrow factor and an optional backflow transformation, we find that the VMC and DMC energies are smooth functions of the sum of the squared coefficients of the initial CI wave function, and that quadratic extrapolations of the non-backflow VMC and backflow DMC energies intersect within uncertainty of the exact total energy. With adequate statistical treatment of quasi-random fluctuations, the extrapolate and intersect with polynomials of order two (XSPOT) method is shown to yield results in agreement with benchmark-quality total and relative energies for the C2, N2, CO2, and H2O molecules, as well as for the C2 molecule in its first electronic singlet excited state, using only small CI expansion sizes. △ Less

Submitted 24 April, 2024; v1 submitted 1 March, 2024; originally announced March 2024.

Comments: 10 pages, 8 figures (+ supplementary: 7 pages, 3 figures)

The Rest-UV Spectral Properties of Dwarf Galaxies at $z\sim2$

Authors: Christopher Snapp-Kolas, Brian Siana, Timothy Gburek, Anahita Alavi, Najmeh Emami, Johan Richard, Daniel P. Stark

Abstract: Rest-UV spectroscopy can constrain properties of the stellar populations, outflows, covering fractions, and can indirectly constrain the Lyman continuum escape fraction of galaxies. Many works have studied the rest-UV spectra of more massive star forming galaxies and low-mass galaxies selected via strong nebular line emission or via Ly$α$ emission. However, studies of rest-UV spectroscopy have yet… ▽ More Rest-UV spectroscopy can constrain properties of the stellar populations, outflows, covering fractions, and can indirectly constrain the Lyman continuum escape fraction of galaxies. Many works have studied the rest-UV spectra of more massive star forming galaxies and low-mass galaxies selected via strong nebular line emission or via Ly$α$ emission. However, studies of rest-UV spectroscopy have yet to be done on an unbiased sample at low mass during the epoch of peak star formation ($z\sim2$). We present a stacked rest-UV spectrum of a complete sample of 16 dwarf galaxies ($\rm \langle log(M^*/M_\odot)\rangle_{median} = 8.2$) at $z\sim2$. The rest-UV Keck/LRIS spectroscopy is complemented by rest-optical Keck/MOSFIRE spectroscopy and Hubble photometry. We find generally larger Ly$α$ equivalent widths ($\rm EW_{Lyα} = 11.2\;$Å) when compared with higher mass ($\rm \langle log(M^*/M_\odot)\rangle_{median} = 10.3$) composites from KBSS ($\rm EW_{Lyα} = -5\;$Å). The average low- and high-ionization absorption line EWs ($\rm EW_{LIS}$ and $\rm EW_{HIS}$, respectively) are weaker ($\rm EW_{LIS}$=-1.18 Å, $\rm EW_{HIS}=$-0.99 Å) in dwarf galaxies than in higher mass galaxies ($\rm EW_{LIS}$=-2.04 Å, $\rm EW_{HIS}=$-1.42 Å). The LIS absorption lines are optically thick and is thus a good tracer of the neutral hydrogen covering fraction. Both higher $\rm EW_{Lyα}$ and lower $\rm EW_{LIS}$ measurements imply that the escape fraction of ionizing radiation is larger in lower-mass galaxies at $z\sim2$. △ Less

Submitted 10 January, 2024; originally announced January 2024.

Comments: 9 pages, 8 figures, 1 table. To be published in MNRAS

Lyman Continuum Emission from AGN at 2.3$\lesssim$z$\lesssim$3.7 in the UVCANDELS Fields

Authors: Brent M. Smith, Rogier A. Windhorst, Harry Teplitz, Matthew Hayes, Marc Rafelski, Mark Dickinson, Vihang Mehta, Nimish P. Hathi, John MacKenty, L. Y. Aaron Yung, Anton M. Koekemoer, Emmaris Soto, Christopher J. Conselice, Ray A. Lucas, Xin Wang, Keunho J. Kim, Anahita Alavi, Norman A. Grogin, Ben Sunnquist, Laura Prichard, Rolf A. Jansen, the UVCANDELS team

Abstract: We present the results of our search for Lyman continuum (LyC) emitting AGN at redshifts 2.3$\lesssim$z$\lesssim$4.9 from HST WFC3 F275W observations in the UVCANDELS fields. We also include LyC emission from AGN using HST WFC3 F225W, F275W, and F336W found in the ERS and HDUV data. We performed exhaustive queries of the Vizier database to locate AGN with high quality spectroscopic redshifts. In t… ▽ More We present the results of our search for Lyman continuum (LyC) emitting AGN at redshifts 2.3$\lesssim$z$\lesssim$4.9 from HST WFC3 F275W observations in the UVCANDELS fields. We also include LyC emission from AGN using HST WFC3 F225W, F275W, and F336W found in the ERS and HDUV data. We performed exhaustive queries of the Vizier database to locate AGN with high quality spectroscopic redshifts. In total, we found 51 AGN that met our criteria within the UVCANDELS and ERS footprints. Of these 51, we find 12 AGN had $\geq$4$σ$ detected LyC flux in the WFC3/UVIS images. Using space- and ground-based data from X-ray to radio, we fit the multi-wavelength photometric data of each AGN to a CIGALE SED and correlate various SED parameters to the LyC flux. KS-tests of the SED parameter distributions for the LyC-detected and non-detected AGN showed they are likely not distinct samples. However, we find that X-ray luminosity, star-formation onset age, and disk luminosity show strong correlations relative to their emitted LyC flux. We also find strong correlation of the LyC flux to several dust parameters, i.e., polar and toroidal dust emission, 6 $μm$ luminosity, and anti-correlation with metallicity and $A_{FUV}$. We simulate the LyC escape fraction ($f_{esc}$) using the CIGALE and IGM transmission models for the LyC-detected AGN and find an average $f_{esc}$$\simeq$18%, weighted by uncertainties. We stack the LyC flux of subsamples of AGN according to the wavelength continuum region in which they are detected and find no significant distinctions in their LyC emission, although our $sub-mm\ detected$ F336W sample shows the brightest stacked LyC flux. These findings indicate that LyC-production and -escape in AGN is more complicated than the simple assumption of thermal emission and a 100% escape fraction. Further testing of AGN models with larger samples than presented here is needed. △ Less

Submitted 5 January, 2024; originally announced January 2024.

Comments: 21 pages, 6 figures, 3 tables. Accepted for publication in The Astrophysical Journal

Neutrino spin oscillations in gravitational fields in noncommutative higher dimensions

Authors: S. A. Alavi, T. Fallahi Serish

Abstract: Investigation of neutrino spin oscillation in the gravitational fields of black holes(BH) is one of the interesting topics in neutrino physics. On the other hand, in recent years, many studies have been devoted to the exploration of different physical phenomena in higher dimensions. Noncommutative geometry has also been in the focus of researchers in the past years to explore deeper and more accur… ▽ More Investigation of neutrino spin oscillation in the gravitational fields of black holes(BH) is one of the interesting topics in neutrino physics. On the other hand, in recent years, many studies have been devoted to the exploration of different physical phenomena in higher dimensions. Noncommutative geometry has also been in the focus of researchers in the past years to explore deeper and more accurate the structure of space time. In this work, the neutrino spin oscillation in the noncommutative higher dimensions gravitational fields of Schwarzschild and Reissner-Nordstrom metrics are studied. The effects of noncommutativity of space are calculated and its role in different dimensions are discussed. Finally upper bounds on noncommutativity parameter are obtained. △ Less

Submitted 20 December, 2023; originally announced December 2023.

Comments: 29 pages, 19 Figures, 1 table

The UV luminosity function at 0.6 < z < 1 from UVCANDELS

Authors: Lei Sun, Xin Wang, Harry I. Teplitz, Vihang Mehta, Anahita Alavi, Marc Rafelski, Rogier A. Windhorst, Claudia Scarlata, Jonathan P. Gardner, Brent M. Smith, Ben Sunnquist, Laura Prichard, Yingjie Cheng, Norman Grogin, Nimish P. Hathi, Matthew Hayes, Anton M. Koekemoer, Bahram Mobasher, Kalina V. Nedkova, Robert O'Connell, Brant Robertson, Sina Taamoli, L. Y. Aaron Yung, Gabriel Brammer, James Colbert , et al. (53 additional authors not shown)

Abstract: UVCANDELS is a HST Cycle-26 Treasury Program awarded 164 orbits of primary ultraviolet (UV) F275W imaging and coordinated parallel optical F435W imaging in four CANDELS fields: GOODS-N, GOODS-S, EGS, and COSMOS, covering a total area of $\sim426$ arcmin$^2$. This is $\sim2.7$ times larger than the area covered by previous deep-field space UV data combined, reaching a depth of about 27 and 28 ABmag… ▽ More UVCANDELS is a HST Cycle-26 Treasury Program awarded 164 orbits of primary ultraviolet (UV) F275W imaging and coordinated parallel optical F435W imaging in four CANDELS fields: GOODS-N, GOODS-S, EGS, and COSMOS, covering a total area of $\sim426$ arcmin$^2$. This is $\sim2.7$ times larger than the area covered by previous deep-field space UV data combined, reaching a depth of about 27 and 28 ABmag ($5σ$ in $0.2"$ apertures) for F275W and F435W, respectively. Along with the new photometric catalogs, we present an analysis of the rest-frame UV luminosity function (LF), relying on our UV-optimized aperture photometry method yielding a factor of $1.5\times$ increase than the H-isophot aperture photometry in the signal-to-noise ratios of galaxies in our F275W imaging. Using well tested photometric redshift measurements we identify 5810 galaxies at redshifts $0.6<z<1$, down to an absolute magnitude of $M_\text{UV} = -14.2$. In order to minimize the effect of uncertainties in estimating the completeness function, especially at the faint-end, we restrict our analysis to sources above $30\%$ completeness, which provides a final sample of 4726 galaxies at $-21.5<M_\text{UV}<-15.5$. We performed a maximum likelihood estimate to derive the best-fit parameters of the UV LF. We report a best-fit faint-end slope of $α= -1.359^{+0.041}_{-0.041}$ at $z \sim 0.8$. Creating sub-samples at $z\sim0.7$ and $z\sim0.9$, we observe a possible evolution of $α$ with redshift. The unobscured UV luminosity density at $M_\text{UV}<-10$ is derived as $ρ_\text{UV}=1.339^{+0.027}_{-0.030}\ (\times10^{26} \text{ergs/s/Hz/Mpc}^3)$ using our best-fit LF parameters. The new F275W and F435 photometric catalogs from UVCANDELS have been made publicly available on the Barbara A. Mikulski Archive for Space Telescopes (MAST). △ Less

Submitted 2 May, 2024; v1 submitted 27 November, 2023; originally announced November 2023.

Comments: 17 pages, 8 figures, Accepted for publication in ApJ

Some basic features of canonical ensemble in noncommutative spaces

Authors: S. A. Alavi

Abstract: We calculate the corrections due to noncommutativity of space on the Hamiltonian and then partition function of the canonical ensemble. We study some basic features of statistical mechanics and thermodynamics including equipartition and virial theorem and energy fluctuations: correspondence with microcanonical ensemble, in the framework of non-commutative canonical ensemble. The corrections impose… ▽ More We calculate the corrections due to noncommutativity of space on the Hamiltonian and then partition function of the canonical ensemble. We study some basic features of statistical mechanics and thermodynamics including equipartition and virial theorem and energy fluctuations: correspondence with microcanonical ensemble, in the framework of non-commutative canonical ensemble. The corrections imposed by noncommutativity of space are derived and the results are discussed. △ Less

Submitted 9 November, 2023; originally announced November 2023.

Comments: 6 pages, no figures

CEERS: 7.7 $μ$m PAH Star Formation Rate Calibration with JWST MIRI

Authors: Kaila Ronayne, Casey Papovich, Guang Yang, Lu Shen, Mark Dickinson, Robert Kennicutt, Anahita Alavi, Pablo Arrabal Haro, Micaela Bagley, Denis Burgarella, Aurélien Le Bail, Eric Bell, Nikko Cleri, Justin Cole, Luca Costantin, Alexander de la Vega, Emanuele Daddi, David Elbaz, Steven Finkelstein, Norman Grogin, Benne Holwerda, Jeyhan Kartaltepe, Allison Kirkpatrick, Anton Koekemoer, Ray Lucas , et al. (11 additional authors not shown)

Abstract: We test the relationship between UV-derived star formation rates (SFRs) and the 7.7 $μ$m polycyclic aromatic hydrocarbon (PAH) luminosities from the integrated emission of galaxies at z ~ 0 - 2. We utilize multi-band photometry covering 0.2 - 160 $μ$m from HST, CFHT, JWST, Spitzer, and Herschel for galaxies in the Cosmic Evolution Early Release Science (CEERS) Survey. We perform spectral energy di… ▽ More We test the relationship between UV-derived star formation rates (SFRs) and the 7.7 $μ$m polycyclic aromatic hydrocarbon (PAH) luminosities from the integrated emission of galaxies at z ~ 0 - 2. We utilize multi-band photometry covering 0.2 - 160 $μ$m from HST, CFHT, JWST, Spitzer, and Herschel for galaxies in the Cosmic Evolution Early Release Science (CEERS) Survey. We perform spectral energy distribution (SED) modeling of these data to measure dust-corrected far-UV (FUV) luminosities, $L_{FUV}$, and UV-derived SFRs. We then fit SED models to the JWST/MIRI 7.7 - 21 $μ$m CEERS data to derive rest-frame 7.7 $μ$m luminosities, $L_{770}$, using the average flux density in the rest-frame MIRI F770W bandpass. We observe a correlation between $L_{770}$ and $L_{FUV}$, where log $L_{770}$ is proportional to (1.27+/-0.04) log $L_{FUV}$. $L_{770}$ diverges from this relation for galaxies at lower metallicities, lower dust obscuration, and for galaxies dominated by evolved stellar populations. We derive a "single-wavelength" SFR calibration for $L_{770}$ which has a scatter from model estimated SFRs (${σ_{ΔSFR}}$) of 0.24 dex. We derive a "multi-wavelength" calibration for the linear combination of the observed FUV luminosity (uncorrected for dust) and the rest-frame 7.7 $μ$m luminosity, which has a scatter of ${σ_{ΔSFR}}$ = 0.21 dex. The relatively small decrease in $σ$ suggests this is near the systematic accuracy of the total SFRs using either calibration. These results demonstrate that the rest-frame 7.7 $μ$m emission constrained by JWST/MIRI is a tracer of the SFR for distant galaxies to this accuracy, provided the galaxies are dominated by star-formation with moderate-to-high levels of attenuation and metallicity. △ Less

Submitted 13 October, 2023; v1 submitted 11 October, 2023; originally announced October 2023.

Comments: 20 pages, 11 figures, 2 tables, submitted to ApJ

Design Verification of the Quantum Control Stack

Authors: Seyed Amir Alavi, Samin Ishtiaq, Nick Johnson, Rojalin Mishra, Dwaraka Oruganti Nagalakshmi, Asher Pearl, Jan Snoeijs

Abstract: This paper describes the verification of the classical software and hardware stack that is used to control cold atom- and superconducting-based quantum computing hardware. The paper serves both as an introduction to quantum computing and to how classical device verification techniques can be employed there. Two main challenges in building a quantum control stack are generating precise deterministi… ▽ More This paper describes the verification of the classical software and hardware stack that is used to control cold atom- and superconducting-based quantum computing hardware. The paper serves both as an introduction to quantum computing and to how classical device verification techniques can be employed there. Two main challenges in building a quantum control stack are generating precise deterministic-timing operations at the edge and scaled-out processing in the middle layer. Both challenges are to do with a certain kind of functional performance correctness. And, as usual, the design lives under tight power, memory and latency constraints. The quantum control stack is a complex interaction of algorithms, software runtimes and digital hardware. We take inspiration from modern software approaches to engineering, such as continuous integration and hardware automation, to quickly ship experimental features to customers in the field. △ Less

Submitted 8 October, 2023; originally announced October 2023.

Comments: In DVCon Europe 2023

ACM Class: D.1; C.1

The Lyman Continuum Escape Fraction of Star-forming Galaxies at $2.4\lesssim z\lesssim3.7$ from UVCANDELS

Authors: Xin Wang, Harry I. Teplitz, Brent M. Smith, Rogier A. Windhorst, Marc Rafelski, Vihang Mehta, Anahita Alavi, Gabriel Brammer, James Colbert, Norman Grogin, Nimish P. Hathi, Anton M. Koekemoer, Laura Prichard, Claudia Scarlata, Ben Sunnquist, Pablo Arrabal Haro, Christopher Conselice, Eric Gawiser, Yicheng Guo, Matthew Hayes, Rolf A. Jansen, Zhiyuan Ji, Ray A. Lucas, Robert O'Connell, Brant Robertson , et al. (52 additional authors not shown)

Abstract: The UltraViolet Imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey Fields (UVCANDELS) survey is a Hubble Space Telescope (HST) Cycle-26 Treasury Program, allocated in total 164 orbits of primary Wide-Field Camera 3 Ultraviolet and Visible light F275W imaging with coordinated parallel Advanced Camera for Surveys F435W imaging, on four of the five premier extragalactic sur… ▽ More The UltraViolet Imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey Fields (UVCANDELS) survey is a Hubble Space Telescope (HST) Cycle-26 Treasury Program, allocated in total 164 orbits of primary Wide-Field Camera 3 Ultraviolet and Visible light F275W imaging with coordinated parallel Advanced Camera for Surveys F435W imaging, on four of the five premier extragalactic survey fields: GOODS-N, GOODS-S, EGS, and COSMOS. We introduce this survey by presenting a thorough search for galaxies at $z\gtrsim2.4$ that leak significant Lyman continuum (LyC) radiation, as well as a stringent constraint on the LyC escape fraction ($f_{\rm esc}$) from stacking the UV images of a population of star-forming galaxies with secure redshifts. Our extensive search for LyC emission and stacking analysis benefit from the catalogs of high-quality spectroscopic redshifts compiled from archival ground-based data and HST slitless spectroscopy, carefully vetted by dedicated visual inspection efforts. We report a sample of five galaxies as individual LyC leaker candidates, showing $f_{\rm esc}^{\rm rel}\gtrsim60\%$ estimated using detailed Monte Carlo analysis of intergalactic medium attenuation. We develop a robust stacking method to apply to five samples of in total 85 non-detection galaxies in the redshift range of $z\in[2.4,3.7]$. Most stacks give tight 2-$σ$ upper limits below $f_{\rm esc}^{\rm rel}<6\%$. A stack for a subset of 32 emission-line galaxies shows tentative LyC leakage detected at 2.9-$σ$, indicating $f_{\rm esc}^{\rm rel}=5.7\%$ at $z\sim2.65$, supporting the key role of such galaxies in contributing to the cosmic reionization and maintaining the UV ionization background. These new F275W and F435W imaging mosaics from UVCANDELS have been made publicly available on the Barbara A. Mikulski Archive for Space Telescopes. △ Less

Submitted 17 August, 2023; originally announced August 2023.

Comments: 33 pages, 21 figures, and 5 tables. Resubmitted after addressing the referee report

UV-Bright Star-Forming Clumps and Their Host Galaxies in UVCANDELS at 0.5 $\leq$ z $\leq$ 1

Authors: Alec Martin, Yicheng Guo, Xin Wang, Anton M. Koekemoer, Marc Rafelski, Harry I. Teplitz, Rogier A. Windhorst, Anahita Alavi, Norman A. Grogin, Laura Prichard, Ben Sunnquist, Daniel Ceverino, Nima Chartab, Christopher J. Conselice, Y. Sophia Dai, Avishai Dekel, Johnathan P. Gardner, Eric Gawiser, Nimish P. Hathi, Matthew J. Hayes, Rolf A. Jansen, Zhiyuan Ji, David C. Koo, Ray A. Lucas, Nir Mandelker , et al. (10 additional authors not shown)

Abstract: Giant star-forming clumps are a prominent feature of star-forming galaxies (SFGs) and contain important clues on galaxy formation and evolution. However, basic demographics of clumps and their host galaxies remain uncertain. Using the HST/WFC3 F275W images from the Ultraviolet Imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (UVCANDELS), we detect and analyze giant sta… ▽ More Giant star-forming clumps are a prominent feature of star-forming galaxies (SFGs) and contain important clues on galaxy formation and evolution. However, basic demographics of clumps and their host galaxies remain uncertain. Using the HST/WFC3 F275W images from the Ultraviolet Imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey (UVCANDELS), we detect and analyze giant star-forming clumps in galaxies at 0.5 $\leq$ z $\leq$ 1, connecting two epochs when clumps are common (at cosmic high-noon, z $\sim$ 2) and rare (in the local universe). We construct a clump sample whose rest-frame 1600 Å luminosity is 3 times higher than the most luminous local HII regions (M$_{UV} \leq -$16 AB). In our sample, 35 $\pm$ 3$\%$ of low-mass galaxies (log[M$_{*}$/M$_{\odot}$] $<$ 10) are clumpy (i.e., containing at least one off-center clump). This fraction changes to 22 $\pm$ 3$\%$ and 22 $\pm$ 4$\%$ for intermediate (10 $\leq$ log[M$_{*}$/M$_{\odot}$] $\leq$ 10.5) and high-mass (log[M$_{*}$/M$_{\odot}$] $>$ 10.5) galaxies in agreement with previous studies. When compared to similar-mass non-clumpy SFGs, low- and intermediate-mass clumpy SFGs tend to have higher SFRs and bluer rest-frame U-V colors, while high-mass clumpy SFGs tend to be larger than non-clumpy SFGs. However, clumpy and non-clumpy SFGs have similar Sérsic index, indicating a similar underlying density profile. Furthermore, we investigate how UV luminosity of star-forming regions correlates with the physical properties of host galaxies. On average, more luminous star-forming regions reside in more luminous, smaller, and/or higher-specific SFR galaxies and are found closer to their hosts' galactic center. △ Less

Submitted 2 October, 2023; v1 submitted 31 July, 2023; originally announced August 2023.

Comments: 21 pages, 13 figures, accepted for publication in ApJ

Journal ref: ApJ 955 106 (2023)

Fraction of Clumpy Star-Forming Galaxies at $0.5\leq z\leq 3$ in UVCANDELS: Dependence on Stellar Mass and Environment

Authors: Zahra Sattari, Bahram Mobasher, Nima Chartab, Daniel D. Kelson, Harry I. Teplitz, Marc Rafelski, Norman A. Grogin, Anton M. Koekemoer, Xin Wang, Rogier A. Windhorst, Anahita Alavi, Laura Prichard, Ben Sunnquist, Jonathan P. Gardner, Eric Gawiser, Nimish P. Hathi, Matthew J. Hayes, Zhiyuan Ji, Vihang Mehta, Brant E. Robertson, Claudia Scarlata, L. Y. Aaron Yung, Christopher J. Conselice, Y. Sophia Dai, Yicheng Guo , et al. (3 additional authors not shown)

Abstract: High-resolution imaging of galaxies in rest-frame UV has revealed the existence of giant star-forming clumps prevalent in high redshift galaxies. Studying these sub-structures provides important information about their formation and evolution and informs theoretical galaxy evolution models. We present a new method to identify clumps in galaxies' high-resolution rest-frame UV images. Using imaging… ▽ More High-resolution imaging of galaxies in rest-frame UV has revealed the existence of giant star-forming clumps prevalent in high redshift galaxies. Studying these sub-structures provides important information about their formation and evolution and informs theoretical galaxy evolution models. We present a new method to identify clumps in galaxies' high-resolution rest-frame UV images. Using imaging data from CANDELS and UVCANDELS, we identify star-forming clumps in an HST/F160W$\leq 25$ AB mag sample of 6767 galaxies at $0.5\leq z\leq 3$ in four fields, GOODS-N, GOODS-S, EGS, and COSMOS. We use a low-pass band filter in Fourier space to reconstruct the background image of a galaxy and detect small-scale features (clumps) on the background-subtracted image. Clumpy galaxies are defined as those having at least one off-center clump that contributes a minimum of 10$\%$ of the galaxy's total rest-frame UV flux. We measure the fraction of clumpy galaxies ($\rm f_{clumpy}$) as a function of stellar mass, redshift, and galaxy environment. Our results indicate that $\rm f_{clumpy}$ increases with redshift, reaching $\sim 65\%$ at $z\sim 1.5$. We also find that $\rm f_{clumpy}$ in low-mass galaxies ($\rm 9.5\leq log(M_*/M_\odot)\leq 10$) is 10$\%$ higher compared to that of their high-mass counterparts ($\rm log(M_*/M_\odot)>10.5$). Moreover, we find no evidence of significant environmental dependence of $\rm f_{clumpy}$ for galaxies at the redshift range of this study. Our results suggest that the fragmentation of gas clouds under violent disk instability remains the primary driving mechanism for clump formation, and incidents common in dense environments, such as mergers, are not the dominant processes. △ Less

Submitted 15 May, 2023; originally announced May 2023.

Comments: 16 pages, 11 figures, 2 tables, accepted for publication in ApJ

xTC: An efficient treatment of three-body interactions in transcorrelated methods

Authors: Evelin Martine Christlmaier, Thomas Schraivogel, Pablo López Ríos, Ali Alavi, Daniel Kats

Abstract: An efficient implementation for approximate inclusion of the three-body operator arising in transcorrelated methods via exclusion of explicit three-body components (xTC) is presented and tested against results in the "HEAT" benchmark set [A. Tajti et al., J. Chem. Phys. 121, 11599 (2004)]. Using relatively modest basis sets and computationally simple methods, total, atomization, and formation ener… ▽ More An efficient implementation for approximate inclusion of the three-body operator arising in transcorrelated methods via exclusion of explicit three-body components (xTC) is presented and tested against results in the "HEAT" benchmark set [A. Tajti et al., J. Chem. Phys. 121, 11599 (2004)]. Using relatively modest basis sets and computationally simple methods, total, atomization, and formation energies within near-chemical accuracy from HEAT results were obtained. The xTC ansatz reduces the nominal scaling of the three-body part of transcorrelation by two orders of magnitude to O(N^5) and can readily be used with almost any quantum chemical correlation method. △ Less

Submitted 3 August, 2023; v1 submitted 14 April, 2023; originally announced April 2023.

Transcorrelated coupled cluster methods. II. Molecular systems

Authors: Thomas Schraivogel, Evelin Martine Christlmaier, Pablo López Ríos, Ali Alavi, Daniel Kats

Abstract: We demonstrate the accuracy of ground-state energies of the transcorrelated Hamiltonian, employing sophisticated Jastrow factors obtained from variational Monte Carlo, together with the coupled cluster and distinguishable cluster methods at the level of singles and doubles excitations. Our results show that already with the cc-pVTZ basis the transcorrelated distinguishable cluster method gets clos… ▽ More We demonstrate the accuracy of ground-state energies of the transcorrelated Hamiltonian, employing sophisticated Jastrow factors obtained from variational Monte Carlo, together with the coupled cluster and distinguishable cluster methods at the level of singles and doubles excitations. Our results show that already with the cc-pVTZ basis the transcorrelated distinguishable cluster method gets close to complete basis limit and near full configuration interaction quality values for relative energies of over thirty atoms and molecules. To gauge the performance in different correlation regimes we also investigate the breaking of the nitrogen molecule with transcorrelated coupled cluster methods. Numerical evidence is presented to further justify an efficient way to incorporate the major effects coming from the three-body integrals without explicitly introducing them into the amplitude equations. △ Less

Submitted 17 May, 2023; v1 submitted 22 March, 2023; originally announced March 2023.

Comments: peer-reviewed and accepted manuscript (J. Chem. Phys.)

Ab Initio Transcorrelated Method enabling accurate Quantum Chemistry on near-term Quantum Hardware

Authors: Werner Dobrautz, Igor O. Sokolov, Ke Liao, Pablo López Ríos, Martin Rahm, Ali Alavi, Ivano Tavernelli

Abstract: Quantum computing is emerging as a new computational paradigm with the potential to transform several research fields, including quantum chemistry. However, current hardware limitations (including limited coherence times, gate infidelities, and limited connectivity) hamper the straightforward implementation of most quantum algorithms and call for more noise-resilient solutions. In quantum chemistr… ▽ More Quantum computing is emerging as a new computational paradigm with the potential to transform several research fields, including quantum chemistry. However, current hardware limitations (including limited coherence times, gate infidelities, and limited connectivity) hamper the straightforward implementation of most quantum algorithms and call for more noise-resilient solutions. In quantum chemistry, the limited number of available qubits and gate operations is particularly restrictive since, for each molecular orbital, one needs, in general, two qubits. In this study, we propose an explicitly correlated Ansatz based on the transcorrelated (TC) approach, which transfers -- without any approximation -- correlation from the wavefunction directly into the Hamiltonian, thus reducing the number of resources needed to achieve accurate results with noisy, near-term quantum devices. In particular, we show that the exact transcorrelated approach not only allows for more shallow circuits but also improves the convergence towards the so-called basis set limit, providing energies within chemical accuracy to experiment with smaller basis sets and, therefore, fewer qubits. We demonstrate our method by computing bond lengths, dissociation energies, and vibrational frequencies close to experimental results for the hydrogen dimer and lithium hydride using just 4 and 6 qubits, respectively. Conventional methods require at least ten times more qubits for the same accuracy. △ Less

Submitted 17 April, 2024; v1 submitted 3 March, 2023; originally announced March 2023.

TREXIO: A File Format and Library for Quantum Chemistry

Authors: Evgeny Posenitskiy, Vijay Gopal Chilkuri, Abdallah Ammar, Michał Hapka, Katarzyna Pernal, Ravindra Shinde, Edgar Josué Landinez Borda, Claudia Filippi, Kosuke Nakano, Otto Kohulák, Sandro Sorella, Pablo de Oliveira Castro, William Jalby, Pablo López Rıós, Ali Alavi, Anthony Scemama

Abstract: TREXIO is an open-source file format and library developed for the storage and manipulation of data produced by quantum chemistry calculations. It is designed with the goal of providing a reliable and efficient method of storing and exchanging wave function parameters and matrix elements, making it an important tool for researchers in the field of quantum chemistry. In this work, we present an ove… ▽ More TREXIO is an open-source file format and library developed for the storage and manipulation of data produced by quantum chemistry calculations. It is designed with the goal of providing a reliable and efficient method of storing and exchanging wave function parameters and matrix elements, making it an important tool for researchers in the field of quantum chemistry. In this work, we present an overview of the TREXIO file format and library. The library consists of a front-end implemented in the C programming language and two different back-ends: a text back-end and a binary back-end utilizing the HDF5 library which enables fast read and write operations. It is compatible with a variety of platforms and has interfaces for the Fortran, Python, and OCaml programming languages. In addition, a suite of tools has been developed to facilitate the use of the TREXIO format and library, including converters for popular quantum chemistry codes and utilities for validating and manipulating data stored in TREXIO files. The simplicity, versatility, and ease of use of TREXIO make it a valuable resource for researchers working with quantum chemistry data. △ Less

Submitted 30 March, 2023; v1 submitted 28 February, 2023; originally announced February 2023.

Comments: 13 pages, 2 figures

Multi-configurational nature of electron correlation within nitrogen vacancy centers in diamond

Authors: Yilin Chen, Tonghuan Jiang, Haoxiang Chen, Erxun Han, Ali Alavi, Kuang Yu, En-Ge Wang, Ji Chen

Abstract: Diamond is a solid-state platform to develop quantum technologies, but it has been a long-standing problem that the current understanding of quantum states in diamond is mostly limited to single-electron pictures. Here, we combine the full configuration interaction quantum Monte Carlo method and the density-matrix functional embedding theory, to achieve unprecedented accuracy in describing the man… ▽ More Diamond is a solid-state platform to develop quantum technologies, but it has been a long-standing problem that the current understanding of quantum states in diamond is mostly limited to single-electron pictures. Here, we combine the full configuration interaction quantum Monte Carlo method and the density-matrix functional embedding theory, to achieve unprecedented accuracy in describing the many-body quantum states of nitrogen vacancy (NV) centers in diamond. More than 30 electrons and 130 molecular orbitals are correlated, which reveals the multi-configurational wavefunction of the many-body quantum states in diamond. The multi-configurational description explains puzzling experimental measurements in intersystem crossing and charge state transition in NV centers in diamond. The calculations not only reproduce the available experimental measurements of the energy gaps between quantum states but also provide new benchmarks for states that are still subject to considerable uncertainty. This study highlights the importance of multi-configurational wavefunction in the many-body quantum states in solids. △ Less

Submitted 27 February, 2023; originally announced February 2023.

Optimizing Jastrow factors for the transcorrelated method

Authors: J. Philip Haupt, Seyed Mohammadreza Hosseini, Pablo Lopez Rios, Werner Dobrautz, Aron Cohen, Ali Alavi

Abstract: We investigate the optimization of flexible tailored real-space Jastrow factors for use in the transcorrelated (TC) method in combination with highly accurate quantum chemistry methods such as initiator full configuration interaction quantum Monte Carlo (FCIQMC). Jastrow factors obtained by minimizing the variance of the TC reference energy are found to yield better, more consistent results than t… ▽ More We investigate the optimization of flexible tailored real-space Jastrow factors for use in the transcorrelated (TC) method in combination with highly accurate quantum chemistry methods such as initiator full configuration interaction quantum Monte Carlo (FCIQMC). Jastrow factors obtained by minimizing the variance of the TC reference energy are found to yield better, more consistent results than those obtained by minimizing the variational energy. We compute all-electron atomization energies for the challenging first-row molecules C2 , CN, N2 , and O2 and find that the TC method yields chemically accurate results using only the cc-pVTZ basis set, roughly matching the accuracy of non-TC calculations with the much larger cc-pV5Z basis set. We also investigate an approximation in which pure three-body excitations are neglected from the TC-FCIQMC dynamics, saving storage and computational cost, and show that it affects relative energies negligibly. Our results demonstrate that the combination of tailored real-space Jastrow factors with the multi-configurational TC-FCIQMC method provides a route to obtaining chemical accuracy using modest basis sets, obviating the need for basis-set extrapolation and composite techniques. △ Less

Submitted 12 May, 2023; v1 submitted 27 February, 2023; originally announced February 2023.

Comments: Submitted to J Chem Phys

CEERS: Spatially Resolved UV and mid-IR Star Formation in Galaxies at 0.2 < z < 2.5: The Picture from the Hubble and James Webb Space Telescopes

Authors: Lu Shen, Casey Papovich, Guang Yang, Jasleen Matharu, Xin Wang, Benjamin Magnelli, David Elbaz, Shardha Jogee, Anahita Alavi, Pablo Arrabal Haro, Bren E. Backhaus, Micaela B. Bagley, Eric F. Bell, Laura Bisigello, Antonello Calabrò, M. C. Cooper, Luca Costantin, Emanuele Daddi, Mark Dickinson, Steven L. Finkelstein, Seiji Fujimoto, Mauro Giavalisco, Norman A. Grogin, Yuchen Guo, Benne W. Holwerda , et al. (16 additional authors not shown)

Abstract: We present the mid-IR (MIR) morphologies for 64 star-forming galaxies at $0.2<z<2.5$ with stellar mass $\rm{M_*>10^{9}~M_\odot}$ using JWST MIRI observations from the Cosmic Evolution Early Release Science survey (CEERS). The MIRI bands span the MIR (7.7--21~$μ$m), enabling us to measure the effective radii ($R_{\rm{eff}}$) and Sérsic indexes of these SFGs at rest-frame 6.2 and 7.7 $μ$m, which con… ▽ More We present the mid-IR (MIR) morphologies for 64 star-forming galaxies at $0.2<z<2.5$ with stellar mass $\rm{M_*>10^{9}~M_\odot}$ using JWST MIRI observations from the Cosmic Evolution Early Release Science survey (CEERS). The MIRI bands span the MIR (7.7--21~$μ$m), enabling us to measure the effective radii ($R_{\rm{eff}}$) and Sérsic indexes of these SFGs at rest-frame 6.2 and 7.7 $μ$m, which contains strong emission from Polycyclic aromatic hydrocarbon (PAH) features, a well-established tracer of star formation in galaxies. We define a ``PAH-band'' as the MIRI bandpass that contains these features at the redshift of the galaxy. We then compare the galaxy morphologies in the PAH-bands to those in rest-frame Near-UV (NUV) using HST ACS/F435W or ACS/F606W and optical/near-IR using HST WFC3/F160W imaging from UVCANDELS and CANDELS, where the NUV-band and F160W trace the profile of (unobscured) massive stars and the stellar continuum, respectively. The $R_{\rm{eff}}$ of galaxies in the PAH-band are slightly smaller ($\sim$10\%) than those in F160W for galaxies with $\rm{M_*\gtrsim10^{9.5}~M_\odot}$ at $z\leq1.2$, but the PAH-band and F160W have a similar fractions of light within 1 kpc. In contrast, the $R_{\rm{eff}}$ of galaxies in the NUV-band are larger, with lower fractions of light within 1 kpc compared to F160W for galaxies at $z\leq1.2$. Using the MIRI data to estimate the $\rm{SFR_{\rm{IR}}}$ surface density, we find the correlation between the $\rm{SFR_{\rm{IR}}}$ surface density and stellar mass has a steeper slope than that of the $\rm{SFR_{\rm{UV}}}$ surface density and stellar mass, suggesting more massive galaxies having increasing amounts of obscured fraction of star formation in their inner regions. This paper demonstrates how the high-angular resolution data from JWST/MIRI can reveal new information about the morphology of obscured-star formation. △ Less

Submitted 2 April, 2023; v1 submitted 13 January, 2023; originally announced January 2023.

Comments: 28 pages, 11 figures, Accepted by ApJ

A spatially resolved analysis of star-formation burstiness by comparing UV and H$α$ in galaxies at z$\sim$1 with UVCANDELS

Authors: Vihang Mehta, Harry I. Teplitz, Claudia Scarlata, Xin Wang, Anahita Alavi, James Colbert, Marc Rafelski, Norman Grogin, Anton Koekemoer, Laura Prichard, Rogier Windhorst, Justin M. Barber, Christopher J. Conselice, Y. Sophia Dai, Jonathan P. Gardner, Eric Gawiser, Yicheng Guo, Nimish Hathi, Pablo Arrabal Haro, Matthew Hayes, Kartheik G. Iyer, Rolf A. Jansen, Zhiyuan Ji, Peter Kurczynski, Maxwell Kuschel , et al. (8 additional authors not shown)

Abstract: The UltraViolet imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey Fields (UVCANDELS) program provides HST/UVIS F275W imaging for four CANDELS fields. We combine this UV imaging with existing HST/near-IR grism spectroscopy from 3D-HST$+$AGHAST to directly compare the resolved rest-frame UV and H$α$ emission for a sample of 979 galaxies at $0.7<z<1.5$ spanning a range in… ▽ More The UltraViolet imaging of the Cosmic Assembly Near-infrared Deep Extragalactic Legacy Survey Fields (UVCANDELS) program provides HST/UVIS F275W imaging for four CANDELS fields. We combine this UV imaging with existing HST/near-IR grism spectroscopy from 3D-HST$+$AGHAST to directly compare the resolved rest-frame UV and H$α$ emission for a sample of 979 galaxies at $0.7<z<1.5$ spanning a range in stellar mass of $10^{8-11.5}$ M$_\odot$. Using a stacking analysis, we perform a resolved comparison between homogenized maps of rest-UV and H$α$ to compute the average UV-to-H$α$ luminosity ratio (an indicator of burstiness in star-formation) as a function of galactocentric radius. We find that galaxies below stellar mass of $\sim$10$^{9.5}$ M$_\odot$, at all radii, have a UV-to-H$α$ ratio higher than the equilibrium value expected from constant star-formation, indicating a significant contribution from bursty star-formation. Even for galaxies with stellar mass $\gtrsim$10$^{9.5}$ M$_\odot$, the UV-to-H$α$ ratio is elevated towards in their outskirts ($R/R_{eff}>1.5$), suggesting that bursty star-formation is likely prevalent in the outskirts of even the most massive galaxies but is likely over-shadowed by their brighter cores. Furthermore, we present the UV-to-H$α$ ratio as a function of galaxy surface brightness, a proxy for stellar mass surface density, and find that regions below $\sim$10$^{7.5}$ M$_\odot$ kpc$^{-2}$ are consistent with bursty star-formation, regardless of their galaxy stellar mass, potentially suggesting that local star-formation is independent of global galaxy properties at the smallest scales. Lastly, we find galaxies at $z>1.1$ to have bursty star-formation regardless of radius or surface brightness. △ Less

Submitted 15 June, 2023; v1 submitted 3 November, 2022; originally announced November 2022.

Comments: 22 pages, 10 figures; accepted to ApJ

Density Matrix Renormalization Group for Transcorrelated Hamiltonians: Ground and Excited States in \emph{ab initio} Systems

Authors: Ke Liao, Huanchen Zhai, Evelin Martine Christlmaier, Thomas Schraivogel, Pablo López Ríos, Daniel Kats, Ali Alavi

Abstract: We present the theory of a density matrix renormalization group (DMRG) algorithm which can solve for both the ground and excited states of non-Hermitian transcorrelated Hamiltonians, and show applications in \emph{ab initio} molecular systems. Transcorrelation (TC) accelerates the basis set convergence rate by including known physics (such as, but not limited to, the electron-electron cusp) in the… ▽ More We present the theory of a density matrix renormalization group (DMRG) algorithm which can solve for both the ground and excited states of non-Hermitian transcorrelated Hamiltonians, and show applications in \emph{ab initio} molecular systems. Transcorrelation (TC) accelerates the basis set convergence rate by including known physics (such as, but not limited to, the electron-electron cusp) in the Jastrow factor used for the similarity transformation. It also improves the accuracy of approximate methods such as coupled cluster singles and doubles (CCSD) as shown by some recent studies. However, the non-Hermiticity of the TC Hamiltonians poses challenges for variational methods like DMRG. Imaginary-time evolution on the matrix product state (MPS) in the DMRG framework has been proposed to circumvent this problem; but this is currently limited to treating the ground state, and has lower efficiency than the time-independent DMRG (TI-DMRG), due to the need to eliminate Trotter errors. In this work, we show that with minimal changes to the existing TI-DMRG algorithm, namely replacing the original Davidson solver with the general Davidson solver to solve the non-Hermitian effective Hamiltonians at each site for a few low-lying right eigenstates, and following the rest of the original DMRG recipe, one can find the ground and excited states with improved efficiency compared to the original DMRG when extrapolating to the infinite bond dimension limit in the same basis set. Accelerated basis set convergence rate is also observed, as expected, within the TC framework. △ Less

Submitted 3 November, 2022; v1 submitted 31 October, 2022; originally announced November 2022.

The Origin of the Observed Lyman alpha EW Distribution of Dwarf Galaxies at z~2

Authors: C. Snapp-Kolas, B. Siana, T. Gburek, A. Alavi, N. Emami, J. Richard, D. P. Stark, C. Scarlata

Abstract: We present a rest-UV selected sample of 32 lensed galaxies at $z\sim 2$ observed with joint Keck/LRIS rest-UV and Keck/MOSFIRE rest-optical spectra behind the clusters Abell 1689, MACS J0717, and MACS J1149. The sample pushes towards the faintest UV luminosities observed ($-19 \le {\rm M_{\rm UV}} \le -17$) at this redshift. The fraction of dwarf galaxies identified as Ly$α$ emitters (… ▽ More We present a rest-UV selected sample of 32 lensed galaxies at $z\sim 2$ observed with joint Keck/LRIS rest-UV and Keck/MOSFIRE rest-optical spectra behind the clusters Abell 1689, MACS J0717, and MACS J1149. The sample pushes towards the faintest UV luminosities observed ($-19 \le {\rm M_{\rm UV}} \le -17$) at this redshift. The fraction of dwarf galaxies identified as Ly$α$ emitters ($\rm EW \ge 20\ \overset{\lower.5em\circ}{\mathrm{A}}$) is ${\rm X_{\rm LAE}}=25^{+15}_{-10}\%$. We use the Balmer lines and UV continuum to estimate the intrinsic EW allowing us to distinguish the effects of the ionizing spectra and Ly$α$ escape fraction on the observed EW distribution. Fainter galaxies ($\rm M_{\rm UV} > -19$) show larger intrinsic EWs and escape fractions than brighter galaxies. Only galaxies with intrinsic EWs greater than 40$\ \overset{\lower.5em\circ}{\mathrm{A}}$ have escape fractions larger than 0.05. We find an anti-correlation between the escape fraction and $\rm A_V$ as well as UV spectral slope. The volumetric escape fraction of our sample is $f_{\rm esc}^{\rm Lyα} = 4.59^{+2.0}_{-1.4}\%$ in agreement with measurements found elsewhere in the literature. About half of the total integrated Ly$α$ luminosity density comes from galaxies with ${\rm EW}_{\rm obs}>20\ \overset{\lower.5em\circ}{\mathrm{A}}$. △ Less

Submitted 31 October, 2022; originally announced November 2022.

Comments: 11 pages, 13 figures

Ferromagnetic domains in the large-$U$ Hubbard model with a few holes: an FCIQMC study

Authors: Sujun Yun, Werner Dobrautz, Hongjun Luo, Vamshi Katukuri, Niklas Libermann, Ali Alavi

Abstract: Two-dimensional Hubbard lattices with two or three holes are investigated as a function of $U$ in the large-$U$ limit. In the so-called Nagaoka limit (one-hole system at infinite $U$), it is known that the Hubbard model exhibits a ferromagnetic ground-state. Here, by means of exact FCIQMC simulations applied to periodic lattices up to 24 sites, we compute spin-spin correlation functions as a funct… ▽ More Two-dimensional Hubbard lattices with two or three holes are investigated as a function of $U$ in the large-$U$ limit. In the so-called Nagaoka limit (one-hole system at infinite $U$), it is known that the Hubbard model exhibits a ferromagnetic ground-state. Here, by means of exact FCIQMC simulations applied to periodic lattices up to 24 sites, we compute spin-spin correlation functions as a function of increasing $U$. The correlation functions clearly demonstrate the onset of ferromagnetic domains, centred on individual holes. The overall total spin of the wavefunctions remain the lowest possible (0, or $\frac{1}{2}$, depending on the number of holes). The ferromagnetic domains appear at interaction strengths comparable to the critical interaction strengths of the Nagaoka transition in finite systems with strictly one hole. The existence of such ferromagnetic domains is the signature of Nagaoka physics in Hubbard systems with a small (but greater than one) number of holes. △ Less

Submitted 21 October, 2022; originally announced October 2022.

The Direct-Method Oxygen Abundance of Typical Dwarf Galaxies at Cosmic High-Noon

Authors: Timothy Gburek, Brian Siana, Anahita Alavi, Najmeh Emami, Johan Richard, William R. Freeman, Daniel P. Stark, Christopher Snapp-Kolas

Abstract: We present a Keck/MOSFIRE, rest-optical, composite spectrum of 16 typical, gravitationally-lensed, star-forming, dwarf galaxies at $1.7 \lesssim z \lesssim 2.6$ ($z_{\rm{mean}}=2.30$), all chosen independent of emission-line strength. These galaxies have a median stellar mass of $\log(M_\ast/\rm{M_\odot})_{\rm{med}} = 8.29^{+0.51}_{-0.43}$ and a median star formation rate of… ▽ More We present a Keck/MOSFIRE, rest-optical, composite spectrum of 16 typical, gravitationally-lensed, star-forming, dwarf galaxies at $1.7 \lesssim z \lesssim 2.6$ ($z_{\rm{mean}}=2.30$), all chosen independent of emission-line strength. These galaxies have a median stellar mass of $\log(M_\ast/\rm{M_\odot})_{\rm{med}} = 8.29^{+0.51}_{-0.43}$ and a median star formation rate of $\rm{SFR_{Hα}^{med} = 2.25^{+2.15}_{-1.26}\ M_\odot\ yr^{-1}}$. We measure the faint, electron-temperature-sensitive, [O III] $λ$4363 emission line at $2.5σ$ ($4.1σ$) significance when considering a bootstrapped (statistical-only) uncertainty spectrum. This yields a direct-method oxygen abundance of $12+\log(\rm{O/H})_{\rm{direct}}=7.88^{+0.25}_{-0.22}$ ($0.15^{+0.12}_{-0.06}\ \rm{Z_\odot}$). We investigate the applicability at high-$z$ of locally-calibrated, oxygen-based, strong-line metallicity relations, finding that the local reference calibrations of arXiv:1805.08224 best reproduce ($\lesssim 0.12$ dex) our composite metallicity at fixed strong-line ratio. At fixed $M_\ast$, our composite is well-represented by the $z \sim 2.3$ direct-method stellar mass$\,-\,$gas-phase metallicity relation (MZR) of arXiv:1907.00013. When comparing to predicted MZRs from the IllustrisTNG and FIRE simulations, we find excellent agreement with the FIRE MZR. Our composite is consistent with no metallicity evolution, at fixed $M_\ast$ and SFR, of the locally-defined fundamental metallicity relation. We measure the doublet ratio [O II] $λ$3729/[O II] $\lambda3726 = 1.56 \pm 0.32$ ($1.51 \pm 0.12$) and a corresponding electron density of $n_e = 1^{+215}_{-0}\ \rm{cm^{-3}}$ ($n_e = 1^{+74}_{-0}\ \rm{cm^{-3}}$) when considering the bootstrapped (statistical-only) error spectrum. This result suggests that lower-mass galaxies have lower densities than higher-mass galaxies at $z \sim 2$. △ Less

Submitted 12 January, 2023; v1 submitted 11 August, 2022; originally announced August 2022.

Comments: 30 pages, 6 figures, 4 tables, abstract slightly abridged, use of PyNeb adopted (replacing use of IRAF) to calculate electron temperatures and densities; change in results is negligible, accepted for publication in ApJ

An excursion into the core of the cluster lens Abell 1689

Authors: Agniva Ghosh, Dominic Adams, Liliya L. R. Williams, Jori Liesenborgs, Anahita Alavi, Claudia Scarlata

Abstract: Abell 1689 is a well studied cluster of galaxies and one of the largest gravitational lens systems ever observed. We have obtained a reconstruction of the cluster Abell 1689 using Grale, a free-form lens inversion method that relies exclusively on the multiple image data. Non-inclusion of any data related to cluster member galaxies ensures an unbiased measure of the mass distribution, which is the… ▽ More Abell 1689 is a well studied cluster of galaxies and one of the largest gravitational lens systems ever observed. We have obtained a reconstruction of the cluster Abell 1689 using Grale, a free-form lens inversion method that relies exclusively on the multiple image data. Non-inclusion of any data related to cluster member galaxies ensures an unbiased measure of the mass distribution, which is the most notable feature of free-form methods like Grale. We used two different sets of multiple image systems from the available strong lensing data - one containing only the secure systems (107 images), and the other containing all available systems, only excluding some very non-secure systems (151 images). For the very well-constrained central $\sim$100 kpc region of the cluster we made detailed comparison of the Grale reconstructed lensing mass and stellar mass retrieved by the Spectral Energy Distribution (SED) fitting software FAST++. We found a light-unaccompanied mass peak in this region, whose existence, while tentative, is favored by the distribution of nearby images that are local maxima in the Fermat potential. However, further tests, using different methodologies are needed to confirm the reality of this feature. If it shown to be real, this light-unaccompanied mass peak is consistent with dark matter self-interaction cross-section $σ\lesssim 1$cm$^2$/g, while being in tension with larger cross-sections. △ Less

Submitted 7 August, 2023; v1 submitted 17 June, 2022; originally announced June 2022.

Comments: 14 pages, 11 figures

Perturbation Calculation of the Uniform Electron Gas with a Transcorrelated Hamiltonian

Authors: Hongjun Luo, Ali Alavi

Abstract: With a transcorrelated Hamiltonian, we perform a many body perturbation (MBPT) calculation on the uniform electron gas in the high density regime. By using a correlation factor optimised for a single determinant Jastrow ansatz, the second order correlation energy is calculated as $\frac{1-\ln2}{π^{2}}\ln(r_{s})-0.05075$. This already reproduces the exact logarithmic term of the random phase approx… ▽ More With a transcorrelated Hamiltonian, we perform a many body perturbation (MBPT) calculation on the uniform electron gas in the high density regime. By using a correlation factor optimised for a single determinant Jastrow ansatz, the second order correlation energy is calculated as $\frac{1-\ln2}{π^{2}}\ln(r_{s})-0.05075$. This already reproduces the exact logarithmic term of the random phase approximation (RPA) result, while the constant term is roughly $7\%$ larger than the RPA one. The close agreement with the RPA method demonstrates that the transcorrelated method offers a viable and potentially efficient method for treating metallic systems. △ Less

Submitted 7 June, 2022; originally announced June 2022.

Neutrino spin oscillations in gravitational fields in higher dimensions

Authors: S. A. Alavi, T. Fallahi. Serish

Abstract: Neutrino physics in one of the most active fields of research with important implications for particle physics, cosmology and astrophysics. On the other hand, motivated by some theories including string theory, formulation of physical theories in more than four space-time dimensions has been the subject of increasing attention in recent years. Interaction of neutrinos with gravitational fields is… ▽ More Neutrino physics in one of the most active fields of research with important implications for particle physics, cosmology and astrophysics. On the other hand, motivated by some theories including string theory, formulation of physical theories in more than four space-time dimensions has been the subject of increasing attention in recent years. Interaction of neutrinos with gravitational fields is one of the interesting phenomena which can lead to transition between different helicity states (spin oscillations). We study neutrino spin oscillations in Schwarzschild and RN backgrounds in higher dimensional gravitational fields. We calculate the transition probability as a function of time and also study the dependence of the oscillation frequency on the orbital radius. The results help us to better understand the behavior of gravity and neutrinos in higher dimensions. △ Less

Submitted 16 January, 2023; v1 submitted 4 June, 2022; originally announced June 2022.

Comments: 12 pages. Replaced with the version published in Modern Physics Letters A

Journal ref: Modern Physics Letters A, Vol. 37, No. 30, 2250200 (2022)

Investigating the Dominant Environmental Quenching Process in UVCANDELS/COSMOS Groups

Authors: Maxwell Kuschel, Claudia Scarlata, Vihang Mehta, Harry I. Teplitz, Marc Rafelski, Xin Wang, Ben Sunnquist, Laura Prichard, Norman Grogin, Rogier Windhorst, Michael Rutkowski, Anahita Alavi, Nima Chartab, Christopher J. Conselice, Y. Sophia Dai, Eric Gawiser, Mauro Giavalisco, Pablo Arrabal Haro, Nimish Hathi, Rolf Jansen, Zhiyuan Ji, Anton Koekemoer, Ray A. Lucas, Kameswara Mantha, Bahram Mobasher , et al. (14 additional authors not shown)

Abstract: We explore how the fraction of quenched galaxies changes in groups of galaxies with respect to the distance to the center of the group, redshift, and stellar mass to determine the dominant process of environmental quenching in $0.2 < z < 0.8$ groups. We use new UV data from the UVCANDELS project in addition to existing multiband photometry to derive new galaxy physical properties of the group gala… ▽ More We explore how the fraction of quenched galaxies changes in groups of galaxies with respect to the distance to the center of the group, redshift, and stellar mass to determine the dominant process of environmental quenching in $0.2 < z < 0.8$ groups. We use new UV data from the UVCANDELS project in addition to existing multiband photometry to derive new galaxy physical properties of the group galaxies from the zCOSMOS 20k Group Catalog. Limiting our analysis to a complete sample of log$(M_*/M_{\odot})>10.56$ group galaxies we find that the probability of being quenched increases slowly with decreasing redshift, diverging from the stagnant field galaxy population. A corresponding analysis on how the probability of being quenched increases with time within groups suggests that the dominant environmental quenching process is characterized by slow ($\sim$Gyr) timescales. We find a quenching time of approximately $4.91^{+0.91}_{-1.47} $Gyrs, consistent with the slow processes of strangulation (Larson et al. 1980) and delayed-then-rapid quenching (Wetzel et al. 2013 arXiv:1206.3571v2 [astro-ph.CO]), although more data are needed to confirm this result. △ Less

Submitted 20 June, 2023; v1 submitted 24 May, 2022; originally announced May 2022.

Journal ref: ApJ 947 17 (2023)

General analytical nuclear force and molecular potential energy surface from full configuration interaction quantum Monte Carlo

Authors: Tonghuan Jiang, Wei Fang, Ali Alavi, Ji Chen

Abstract: Full configuration interaction quantum Monte Carlo (FCIQMC) is a state-of-the-art stochastic electronic structure method, providing a methodology to compute FCI-level state energies of molecular systems within a quantum chemical basis. However, especially to probe {\em dynamics} at the FCIQMC level, it is necessary to devise more efficient schemes to produce nuclear forces and potential energy sur… ▽ More Full configuration interaction quantum Monte Carlo (FCIQMC) is a state-of-the-art stochastic electronic structure method, providing a methodology to compute FCI-level state energies of molecular systems within a quantum chemical basis. However, especially to probe {\em dynamics} at the FCIQMC level, it is necessary to devise more efficient schemes to produce nuclear forces and potential energy surfaces (PES) from FCIQMC. In this work, we derive the general formula for nuclear force from FCIQMC, and clarify different contributions of the total force. This method to obtain FCIQMC forces eliminates previous restrictions, and can be used with frozen core approximation and free selection of orbitals, making it promising for more efficient nuclear force calculations. After numerical check of this procedure on the binding curve of N$_2$ molecule, we use the FCIQMC energy and force to obtain the full-dimensional ground state PES of water molecule via Gaussian processes regression. The new water FCIQMC PES can be used as the basis for H$_2$O ground state nuclear dynamics, structure optimization, and rotation-vibrational spectrum calculation. △ Less

Submitted 28 April, 2022; originally announced April 2022.

The average dust attenuation curve at z~1.3 based on HST grism surveys

Authors: A. J. Battisti, M. B. Bagley, I. Baronchelli, Y. -S. Dai, A. L. Henry, M. A. Malkan, A. Alavi, D. Calzetti, J. Colbert, P. J. McCarthy, V. Mehta, M. Rafelski, C. Scarlata, I. Shivaei, E. Wisnioski

Abstract: We present the first characterisation of the average dust attenuation curve at $z\sim1.3$ by combining rest-frame ultraviolet through near-IR photometry with Balmer decrement ($\mathrm{H}α$/$\mathrm{H}β$) constraints for $\sim$900 galaxies with $8\lesssim\log (M_\star /M_\odot)<10.2$ at $0.75<z<1.5$ in the HST WFC3 IR Spectroscopic Parallel (WISP) and 3D-HST grism surveys. Using galaxies in SDSS,… ▽ More We present the first characterisation of the average dust attenuation curve at $z\sim1.3$ by combining rest-frame ultraviolet through near-IR photometry with Balmer decrement ($\mathrm{H}α$/$\mathrm{H}β$) constraints for $\sim$900 galaxies with $8\lesssim\log (M_\star /M_\odot)<10.2$ at $0.75<z<1.5$ in the HST WFC3 IR Spectroscopic Parallel (WISP) and 3D-HST grism surveys. Using galaxies in SDSS, we establish that the ($\mathrm{H}α$+[NII])/[OIII] line ratio and stellar mass are good proxies for the Balmer decrement in low-spectral resolution grism data when only upper-limits on $\mathrm{H}β$ are available and/or $\mathrm{H}α$ is blended with [NII]. The slope of the $z\sim1.3$ attenuation curve ($A(0.15μm)/A(V)=3.15$) and its normalization ($R_V=3.26$) lie in-between the values found for $z=0$ and $z\sim2$ dust attenuation curves derived with similar methods. These provide supporting evidence that the average dust attenuation curve of star forming galaxies evolves continuously with redshift. The $z\sim1.3$ curve has a mild 2175Å feature (bump amplitude, $E_b=0.83$; $\sim$25% that of the MW extinction curve), which is comparable to several other studies at $0<z\lesssim3$, and suggests that the average strength of this feature may not evolve significantly with redshift. The methods we develop to constrain dust attenuation from HST grism data can be applied to future grism surveys with JWST, Euclid, and RST. These new facilities will detect millions of emission line galaxies and offer the opportunity to significantly improve our understanding of how and why dust attenuation curves evolve. △ Less

Submitted 12 April, 2022; originally announced April 2022.

Comments: 21 pages, 13 figures, 4 tables. Accepted for publication in MNRAS

Analysis of quantum decay law: Is quantum tunneling really exponential?

Authors: M. S. Hosseini-Ghalehni, B. Azadegan, S. A. Alavi

Abstract: The exponential decay law is well established since its first derivation in 1928, however it is not exact but only an approximate description. In recent years some experimental and theoretical indications for non-exponential decay have been documented. First we solve analytically the time-dependent Schrödinger equation in one dimension for a potential consisting of an infinite wall plus a rectangu… ▽ More The exponential decay law is well established since its first derivation in 1928, however it is not exact but only an approximate description. In recent years some experimental and theoretical indications for non-exponential decay have been documented. First we solve analytically the time-dependent Schrödinger equation in one dimension for a potential consisting of an infinite wall plus a rectangular barrier with finite width and also a cut harmonic oscillator potential by considering it as a sequence of square potentials. Then using the staggered Leap-Frog method, we solve the time-dependent Schrödinger equation for the cut harmonic oscillator potential. In both methods, time dependence of the survival probability of the particle and the decay parameter λ are analyzed. The results exhibit non-exponential behavior for survival probability at short and intermediate times. △ Less

Submitted 28 January, 2023; v1 submitted 18 March, 2022; originally announced March 2022.

Comments: 16 pages, 16 Figures. Replaced with the version pulished in The European Physical Journal Plus

Journal ref: Eur. Phys. J. Plus 137, 1326 (2022)

Casimir effect and Lorentz invariance violation

Authors: S. A. Alavi

Abstract: The Casimir effect is one of the most direct manifestations of the existence of the vacuum quantum fluctuations, discovered by H. B Casimir in 1948. On the other hand, Lorentz invariance is one of the main and basic concepts in special relativity, which states that, the laws of physics are invariant under Lorentz transformation. In this work, we calculate the corrections imposed by LIV on Casimir… ▽ More The Casimir effect is one of the most direct manifestations of the existence of the vacuum quantum fluctuations, discovered by H. B Casimir in 1948. On the other hand, Lorentz invariance is one of the main and basic concepts in special relativity, which states that, the laws of physics are invariant under Lorentz transformation. In this work, we calculate the corrections imposed by LIV on Casimir effect (force). This may provide a direct probe to test LIV in nature. △ Less

Submitted 25 February, 2022; originally announced March 2022.

Comments: 4 pages

Journal ref: PoS(EPS-HEP2021)621. Contribution to The European Physical Society Conference on High Energy Physics (EPS_HEP2021) 26-30 July 2021

Lorentz invariance violation (LIV) in some basic phenomena in quantum physics

Authors: Z. Shafeei, S. A. Alavi

Abstract: Lorentz symmetry is one of the cornerstone of both general relativity and the standard model of particle physics. We study the violation of Lorentz symmetry in some basic phenomena in atomic physics. Using the Green's function, and the source 4-current, the differential equation of 4-vector of electromagnetic potential is solved and the modified coulomb potential is obtained by some researchers. U… ▽ More Lorentz symmetry is one of the cornerstone of both general relativity and the standard model of particle physics. We study the violation of Lorentz symmetry in some basic phenomena in atomic physics. Using the Green's function, and the source 4-current, the differential equation of 4-vector of electromagnetic potential is solved and the modified coulomb potential is obtained by some researchers. Using modified Coulomb potential, we find the corrections due to LIV on the spectrum of Hydrogen and Helium atoms. We also investigate the consequences of LIV on Stark, Zeeman and Spin orbit effects and obtain some upper bounds for the LIV coefficients. △ Less

Submitted 25 February, 2022; originally announced February 2022.

Comments: 4 pages. Contribution to The European Physical Society Conference on High Energy Physics (EPS_HEP2021) 26-30 July 2021

Journal ref: PoS(EPS-HEP2021)620

Performance of a one-parameter correlation factor for transcorrelation: the Li-Ne total energies and ionization potentials

Authors: Werner Dobrautz, Aron J. Cohen, Ali Alavi, Emmanuel Giner

Abstract: In this work we investigate the performance of a recently proposed transcorrelated (TC) approach based on a single-parameter correlation factor [JCP, 154, 8, 2021] for systems involving more than two electrons. The benefit of such an approach relies on its simplicity as efficient numerical-analytical schemes can be set up to compute the two- and three-body integrals occuring in the effective TC Ha… ▽ More In this work we investigate the performance of a recently proposed transcorrelated (TC) approach based on a single-parameter correlation factor [JCP, 154, 8, 2021] for systems involving more than two electrons. The benefit of such an approach relies on its simplicity as efficient numerical-analytical schemes can be set up to compute the two- and three-body integrals occuring in the effective TC Hamiltonian. To obtain accurate ground state energies within a given basis set, the present TC scheme is coupled to the recently proposed TC-full configuration interaction quantum Monte Carlo method [JCP, 151, 6, 2019]. We report ground state total energies on the Li-Ne series, together with their first cations, computed in increasing large basis sets and compare to more elaborate correlation factors involving electron-electron-nucleus coordinates. Numerical results on the Li-Ne ionization potentials show that the use of the single-parameter correlation factor brings on average only a slightly lower accuracy (1.2 mH) in a triple-zeta quality basis set with respect to a more sophisticated correlation factor. However, already using a quadruple-zeta quality basis set yields results within chemical accuracy to complete basis set limit results when using this novel single-parameter correlation factor. △ Less

Submitted 21 February, 2022; originally announced February 2022.

General embedded cluster protocol for accurate modeling of oxygen vacancies in metal-oxides

Authors: Benjamin Xu Shi, Venkat Kapil, Andrea Zen, Ji Chen, Ali Alavi, Angelos Michaelides

Abstract: The O vacancy (Ov) formation energy, $E_\textrm{Ov}$, is an important property of a metal-oxide, governing its performance in applications such as fuel cells or heterogeneous catalysis. These defects are routinely studied with density functional theory (DFT). However, it is well-recognized that standard DFT formulations (e.g. the generalized gradient approximation) are insufficient for modeling th… ▽ More The O vacancy (Ov) formation energy, $E_\textrm{Ov}$, is an important property of a metal-oxide, governing its performance in applications such as fuel cells or heterogeneous catalysis. These defects are routinely studied with density functional theory (DFT). However, it is well-recognized that standard DFT formulations (e.g. the generalized gradient approximation) are insufficient for modeling the Ov, requiring higher levels of theory. The embedded cluster method offers a promising approach to compute $E_\textrm{Ov}$ accurately, giving access to all electronic structure methods. Central to this approach is the construction of quantum(-mechanically treated) clusters placed within suitable embedding environments. Unfortunately, current approaches to constructing the quantum clusters either require large system sizes, preventing application of high-level methods, or require significant manual input, preventing investigations of multiple systems simultaneously. In this work, we present a systematic and general quantum cluster design protocol that can determine small converged quantum clusters for studying the Ov in metal-oxides with accurate methods such as local coupled cluster with single, double and perturbative triple excitations [CCSD(T)]. We apply this protocol to study the Ov in the bulk and surface planes of rutile TiO2 and rocksalt MgO, producing the first accurate and well-converged determinations of $E_\textrm{Ov}$ with this method. These reference values are used to benchmark exchange-correlation functionals in DFT and we find that all studied functionals underestimate $E_\textrm{Ov}$, with the average error decreasing along the rungs of Jacob's ladder. This protocol is automatable for high-throughput calculations and can be generalized to study other point defects or adsorbates. △ Less

Submitted 31 March, 2022; v1 submitted 9 February, 2022; originally announced February 2022.

Comments: Accepted for publication in Journal of Chemical Physics on 2022-03-10

Ab initio wavefunction analysis of electron removal quasi-particle state of NdNiO$_2$ with fully correlated quantum chemical methods

Authors: Vamshi M. Katukuri, Nikolay A. Bogdanov, Ali Alavi

Abstract: The discovery of superconductivity in hole-doped infinite-layer NdNiO$_2$ -- a transition metal (TM) oxide that is both isostructural and isoelectronic to cuprate superconductors -- has lead to renewed enthusiasm in the hope of understanding the origin of unconventional superconductivity. Here, we investigate the electron-removal states in infinite-layered Ni$^{1+}$ oxide, NdNiO$_2$, which mimics… ▽ More The discovery of superconductivity in hole-doped infinite-layer NdNiO$_2$ -- a transition metal (TM) oxide that is both isostructural and isoelectronic to cuprate superconductors -- has lead to renewed enthusiasm in the hope of understanding the origin of unconventional superconductivity. Here, we investigate the electron-removal states in infinite-layered Ni$^{1+}$ oxide, NdNiO$_2$, which mimics hole-doping, with the state-of-the-art many-body multireference quantum chemistry methods. From the analysis of the many-body wavefunction, we find that the hole-doped $d^8$ ground state of NdNiO$_2$ is very different from the $d^8$ ground state in isostructural cuprate analog CaCuO$_2$, although the parent $d^9$ ground states are for the most part identical. We show that the doped hole in NdNiO$_2$ mainly localizes on the Ni $3d_{x^2-y^2}$ orbital to form a closed-shell singlet, and this singlet configuration contributes to $\sim$40% of the wavefunction. In contrast, in CaCuO$_2$ the Zhang-Rice singlet configurations contribute to $\sim$65% of the wavefunction. With the help of the quantum information concept of entanglement entropy, we quantify the different types of electronic correlations in the nickelate and cuprate compounds and find that the dynamic radial-type correlations within the Ni $d$ manifold are persistent in hole-doped NdNiO$_2$. As a result, the $d^8$ multiplet effects are stronger and the additional hole foot-print is more three-dimensional in NdNiO$_2$. Our analysis shows that the most commonly used three-band Hubbard model employed to express the doped scenario in cuprates represents $\sim$90% of the $d^8$ wavefunction for CaCuO$_2$, but such a model grossly approximates the $d^8$ wavefunction for NdNiO$_2$ as it only stands for $\sim$60% of the wavefunction. △ Less

Submitted 14 January, 2022; originally announced January 2022.

Comments: 10 pages with 4 figures

Orders of magnitude reduction in the computational overhead for quantum many-body problems on quantum computers via an exact transcorrelated method

Authors: Igor O. Sokolov, Werner Dobrautz, Hongjun Luo, Ali Alavi, Ivano Tavernelli

Abstract: Transcorrelated methods provide an efficient way of partially transferring the description of electronic correlations from the ground state wavefunction directly into the underlying Hamiltonian. In particular, Dobrautz et al. [Phys. Rev. B, 99(7), 075119, (2019)] have demonstrated that the use of momentum-space representation, combined with a non-unitary similarity transformation, results in a Hub… ▽ More Transcorrelated methods provide an efficient way of partially transferring the description of electronic correlations from the ground state wavefunction directly into the underlying Hamiltonian. In particular, Dobrautz et al. [Phys. Rev. B, 99(7), 075119, (2019)] have demonstrated that the use of momentum-space representation, combined with a non-unitary similarity transformation, results in a Hubbard Hamiltonian that possesses a significantly more compact ground state wavefunction, dominated by a single Slater determinant. This compactness/single-reference character greatly facilitates electronic structure calculations. As a consequence, however, the Hamiltonian becomes non-Hermitian, posing problems for quantum algorithms based on the variational principle. We overcome these limitations with the ansatz-based quantum imaginary time evolution algorithm and apply the transcorrelated method in the context of digital quantum computing. We demonstrate that this approach enables up to 4 orders of magnitude more accurate and compact solutions in various instances of the Hubbard model at intermediate interaction strength ($U/t=4$), enabling the use of shallower quantum circuits for wavefunction ansatzes. In addition, we propose a more efficient implementation of the quantum imaginary time evolution algorithm in quantum circuits that is tailored to non-Hermitian problems. To validate our approach, we perform hardware experiments on the ibmq_lima quantum computer. Our work paves the way for the use of exact transcorrelated methods for the simulations of ab initio systems on quantum computers. △ Less

Submitted 2 March, 2023; v1 submitted 9 January, 2022; originally announced January 2022.

Comments: 21 pages, 6 figures

Combined unitary and symmetric group approach applied to low-dimensional spin systems

Authors: Werner Dobrautz, Vamshi M. Katukuri, Nikolay A. Bogdanov, Daniel Kats, Giovanni Li Manni, Ali Alavi

Abstract: A novel combined unitary and symmetric group approach is used to study the spin-$\frac{1}{2}$ Heisenberg model and related Fermionic systems in a spin-adapted representation, using a linearly-parameterised Ansatz for the many-body wave function. We show that a more compact ground state wave function representation is obtained when combining the symmetric group, $\mathcal{S}_n$, in the form of perm… ▽ More A novel combined unitary and symmetric group approach is used to study the spin-$\frac{1}{2}$ Heisenberg model and related Fermionic systems in a spin-adapted representation, using a linearly-parameterised Ansatz for the many-body wave function. We show that a more compact ground state wave function representation is obtained when combining the symmetric group, $\mathcal{S}_n$, in the form of permutations of the underlying lattice site ordering, with the cumulative spin-coupling based on the unitary group, $\mathrm{U}(n)$. In one-dimensional systems the observed compression of the wave function is reminiscent of block-spin renormalization group approaches, and allows us to study larger lattices (here taken up to 80 sites) with the spin-adapted full configuration interaction quantum Monte Carlo method, which benefits from the sparsity of the Hamiltonian matrix and the corresponding sampled eigenstates that emerge from the reordering. We find that in an optimal lattice ordering the configuration state function with highest weight already captures with high accuracy the spin-spin correlation function of the exact ground state wave function. This feature is found for more general lattice models, such as the Hubbard model, and ab initio quantum chemical models, in this work exemplified by a one-dimensional hydrogen chain. We also provide numerical evidence that the optimal lattice ordering for the unitary group approach is not generally equivalent to the optimal ordering obtained for methods based on matrix-product states, such as the density-matrix renormalization group approach. △ Less

Submitted 17 December, 2021; originally announced December 2021.

Dynamical Noncommutative Graphene

Authors: Ilyas Haouam, S Ali Alavi

Abstract: We study graphene in a two-dimensional dynamical noncommutative space in the presence of a constant magnetic field. The model is solved using perturbation theory and to the second order of perturbation. The energy levels of the system are calculated and the corresponding eigenstates are obtained. For all cases, the energy shift depends on the dynamical noncommutative parameter τ. Using the accurac… ▽ More We study graphene in a two-dimensional dynamical noncommutative space in the presence of a constant magnetic field. The model is solved using perturbation theory and to the second order of perturbation. The energy levels of the system are calculated and the corresponding eigenstates are obtained. For all cases, the energy shift depends on the dynamical noncommutative parameter τ. Using the accuracy of energy measurement we put an upper bound on the noncommutativity parameter τ. In addition, we investigate some of the thermodynamic quantities of the system at zero temperature limit and extreme relativistic case, which reveals interesting differences between commutative and dynamical noncommutative spaces. △ Less

Submitted 21 October, 2021; originally announced October 2021.

Comments: 13 pages

Journal ref: International Journal of Modern Physics AVol. 37, No. 10, 2250054 (2022)

Spectroscopically Identified Emission Line Galaxy Pairs in the WISP survey

Authors: Y. Sophia Dai, Matthew M. Malkan, Harry I. Teplitz, Claudia Scarlata, Anahita Alavi, Hakim Atek, Micaela Bagley, Ivano Baronchelli, Andrew Battisti, Andrew J Bunker, Nimish P. Hathi, Alaina Henry, Jiasheng Huang, Gaoxiang Jin, Zijian Li, Crystal Martin, Vihang Mehta, John Phillips, Marc Rafelski, Michael Rutkowski, Hai Xu, Cong K Xu, Anita Zanella

Abstract: We identify a sample of spectroscopically measured emission line galaxy (ELG) pairs up to z=1.6 from the WFC3 Infrared Spectroscopic Parallels (WISP) survey. WISP obtained slitless, near-infrared grism spectroscopy along with direct imaging in the J and H bands by observing in the pure-parallel mode with the Wide Field Camera Three (WFC3) on the Hubble Space Telescope (HST). From our search of 419… ▽ More We identify a sample of spectroscopically measured emission line galaxy (ELG) pairs up to z=1.6 from the WFC3 Infrared Spectroscopic Parallels (WISP) survey. WISP obtained slitless, near-infrared grism spectroscopy along with direct imaging in the J and H bands by observing in the pure-parallel mode with the Wide Field Camera Three (WFC3) on the Hubble Space Telescope (HST). From our search of 419 WISP fields covering an area of ~0.5 deg$^{2}$, we find 413 ELG pair systems, mostly Halpha emitters. We then derive reliable star formation rates (SFRs) based on the attenuation-corrected Halpha fluxes. Compared to isolated galaxies, we find an average SFR enhancement of 40%-65%, which is stronger for major pairs and pairs with smaller velocity separations (Delta_v < 300 km/s). Based on the stacked spectra from various subsamples, we study the trends of emission line ratios in pairs, and find a general consistency with enhanced lower-ionization lines. We study the pair fraction among ELGs, and find a marginally significant increase with redshift $f \propto (1+z)^α$, where the power-law index α=0.58$\pm$0.17 from $z\sim$0.2 to $z\sim$1.6. The fraction of Active galactic Nuclei (AGNs), is found to be the same in the ELG pairs as compared to isolated ELGs. △ Less

Submitted 14 October, 2021; originally announced October 2021.

Comments: 22 pages, 9 figures, 6 tables, accepted for publication in ApJ

Lyman Continuum Galaxy Candidates in COSMOS

Authors: Laura J. Prichard, Marc Rafelski, Jeff Cooke, Uros Mestric, Robert Bassett, Emma V. Ryan-Weber, Ben Sunnquist, Anahita Alavi, Nimish Hathi, Xin Wang, Mitchell Revalski, Varun Bajaj, John M. O'Meara, Lee Spitler

Abstract: Star-forming galaxies are the sources likely to have reionized the universe. As we cannot observe them directly due to the opacity of the intergalactic medium at $z\gtrsim5$, we study $z\sim3\text{--}5$ galaxies as proxies to place observational constraints on cosmic reionization. Using new deep \textit{Hubble Space Telescope} rest-frame UV F336W and F435W imaging (30-orbit, $\sim40$~arcmin$^2$,… ▽ More Star-forming galaxies are the sources likely to have reionized the universe. As we cannot observe them directly due to the opacity of the intergalactic medium at $z\gtrsim5$, we study $z\sim3\text{--}5$ galaxies as proxies to place observational constraints on cosmic reionization. Using new deep \textit{Hubble Space Telescope} rest-frame UV F336W and F435W imaging (30-orbit, $\sim40$~arcmin$^2$, $\sim29\text{--}30$~mag depth at 5$σ$), we attempt to identify a sample of Lyman continuum galaxies (LCGs). These are individual sources that emit ionizing flux below the Lyman break ($<912~\textÅ$). This population would allow us to constrain cosmic reionization parameters such as the number density and escape fraction ($f_{\rm esc}$) of ionizing sources. We compile a comprehensive parent sample that does not rely on the Lyman-break technique for redshifts. We present three new spectroscopic candidates at $z\sim3.7\text{--}4.4$, and 32 new photometric candidates. The high-resolution multi-band HST imaging and new Keck/Low Resolution Imaging Spectrometer (LRIS) redshifts make these promising spectroscopic LCG candidates. Using both a traditional and probabilistic approach, we find the most likely $f_{\rm esc}$ values for the three spectroscopic LCG candidates are $>100\%$, and therefore not physical. We are unable to confirm the true nature of these sources with the best available imaging and direct blue Keck/LRIS spectroscopy. More spectra, especially from the new class of 30 m telescopes, will be required to build a statistical sample of LCGs to place firm observational constraints on cosmic reionization. △ Less

Submitted 14 February, 2022; v1 submitted 13 October, 2021; originally announced October 2021.

Comments: Accepted for publication in ApJ. 38 pages (28 in main body), 14 figures