-
From high-dimensional committors to reactive insights
Authors:
Nils E. Strand,
Schuyler B. Nicholson,
Hadrien Vroylandt,
Todd R. Gingrich
Abstract:
Transition path theory (TPT) offers a powerful formalism for extracting the rate and mechanism of rare dynamical transitions between metastable states. Most applications of TPT either focus on systems with modestly sized state spaces or use collective variables to try to tame the curse of dimensionality. Increasingly, expressive function approximators like neural networks and tensor networks have…
▽ More
Transition path theory (TPT) offers a powerful formalism for extracting the rate and mechanism of rare dynamical transitions between metastable states. Most applications of TPT either focus on systems with modestly sized state spaces or use collective variables to try to tame the curse of dimensionality. Increasingly, expressive function approximators like neural networks and tensor networks have shown promise in computing the central object of TPT, the committor function, even in very high dimensional systems. That progress prompts our consideration of how one could use such a high dimensional function to extract mechanistic insight. Here, we present and illustrate a straightforward but powerful way to track how individual dynamical coordinates evolve during a reactive event. The strategy, which involves marginalizing the reactive ensemble, naturally captures the evolution of the dynamical coordinate's distribution, not just its mean reactive behavior.
△ Less
Submitted 9 August, 2024; v1 submitted 3 June, 2024;
originally announced June 2024.
-
Computing time-periodic steady-state currents via the time evolution of tensor network states
Authors:
Nils E. Strand,
Hadrien Vroylandt,
Todd R. Gingrich
Abstract:
We present an approach based upon binary tree tensor network (BTTN) states for computing steady-state current statistics for a many-particle 1D ratchet subject to volume exclusion interactions. The ratcheted particles, which move on a lattice with periodic boundary conditions subject to a time-periodic drive, can be stochastically evolved in time to sample representative trajectories via a Gillesp…
▽ More
We present an approach based upon binary tree tensor network (BTTN) states for computing steady-state current statistics for a many-particle 1D ratchet subject to volume exclusion interactions. The ratcheted particles, which move on a lattice with periodic boundary conditions subject to a time-periodic drive, can be stochastically evolved in time to sample representative trajectories via a Gillespie method. In lieu of generating realizations of trajectories, a BTTN state can variationally approximate a distribution over the vast number of many-body configurations. We apply the density matrix renormalization group (DMRG) algorithm to initialize BTTN states, which are then propagated in time via the time-dependent variational principle (TDVP) algorithm to yield the steady-state behavior, including the effects of both typical and rare trajectories. The application of the methods to ratchet currents is highlighted in a companion letter, but the approach extends naturally to other interacting lattice models with time-dependent driving. Though trajectory sampling is conceptually and computationally simpler, we discuss situations for which the BTTN TDVP strategy could be more favorable.
△ Less
Submitted 11 January, 2022;
originally announced January 2022.
-
Using tensor network states for multi-particle Brownian ratchets
Authors:
Nils E. Strand,
Hadrien Vroylandt,
Todd R. Gingrich
Abstract:
The study of Brownian ratchets has taught how time-periodic driving supports a time-periodic steady state that generates nonequilibrium transport. When a single particle is transported in one dimension, it is possible to rationalize the current in terms of the potential, but experimental efforts have ventured beyond that single-body case to systems with many interacting carriers. Working with a la…
▽ More
The study of Brownian ratchets has taught how time-periodic driving supports a time-periodic steady state that generates nonequilibrium transport. When a single particle is transported in one dimension, it is possible to rationalize the current in terms of the potential, but experimental efforts have ventured beyond that single-body case to systems with many interacting carriers. Working with a lattice model of volume-excluding particles in one dimension, we analyze the impact of interactions on a flashing ratchet's current. To surmount the many-body problem, we employ the time-dependent variational principle with a binary tree tensor network, methods discussed at length in a companion paper. Rather than propagating individual trajectories, the tensor network approach propagates a distribution over many-body configurations via a controllable variational approximation. The calculations, which reproduce Gillespie trajectory sampling, identify and explain a shift in the frequency of maximum current to higher driving frequency as the lattice occupancy increases.
△ Less
Submitted 12 January, 2022; v1 submitted 10 January, 2022;
originally announced January 2022.
-
Current inversion in a periodically driven two-dimensional Brownian ratchet
Authors:
Nils E. Strand,
Rueih-Sheng Fu,
Todd R. Gingrich
Abstract:
It is well-known that Brownian ratchets can exhibit current reversals, wherein the sign of the current switches as a function of the driving frequency. We introduce a spatial discretization of such a two-dimensional Brownian ratchet to enable spectral methods that efficiently compute those currents. These discrete-space models provide a convenient way to study the Markovian dynamics conditioned up…
▽ More
It is well-known that Brownian ratchets can exhibit current reversals, wherein the sign of the current switches as a function of the driving frequency. We introduce a spatial discretization of such a two-dimensional Brownian ratchet to enable spectral methods that efficiently compute those currents. These discrete-space models provide a convenient way to study the Markovian dynamics conditioned upon generating particular values of the currents. By studying such conditioned processes, we demonstrate that low-frequency negative values of current arise from typical events and high-frequency positive values of current arises from rare events. We demonstrate how these observations can inform the sculpting of time-dependent potential landscapes with a specific frequency response.
△ Less
Submitted 6 August, 2020; v1 submitted 31 March, 2020;
originally announced March 2020.
-
Imaginary shift in CASPT2 nuclear gradient and derivative coupling theory
Authors:
Jae Woo Park,
Rachael Al-Saadon,
Nils E. Strand,
Toru Shiozaki
Abstract:
We report the analytical nuclear gradient theory for complete active space second-order perturbation theory (CASPT2) with imaginary shift, which is commonly used to avoid divergence of the perturbation expression. Our formulation is based on the Lagrangian approach and is an extension of the algorithm for CASPT2 nuclear gradients with real shift. The working equations are derived and implemented i…
▽ More
We report the analytical nuclear gradient theory for complete active space second-order perturbation theory (CASPT2) with imaginary shift, which is commonly used to avoid divergence of the perturbation expression. Our formulation is based on the Lagrangian approach and is an extension of the algorithm for CASPT2 nuclear gradients with real shift. The working equations are derived and implemented into an efficient parallel program. Numerical examples are presented for the ground- and excited-state geometries and conical intersections of a green fluorescent protein model chromophore, $p$-HBDI$^-$. We also report timing benchmarks with adenine, $p$-HBDI$^-$, and iron porphyrin. It is demonstrated that the energies and geometries obtained with the imaginary shift improve accuracy at a minor additional cost which is mainly associated with evaluating the effective density matrix elements for the imaginary shift term.
△ Less
Submitted 14 April, 2019;
originally announced April 2019.
-
The Relationship of Active Galactic Nuclei & Quasars With Their Local Galaxy Environment
Authors:
Natalie Erin Strand
Abstract:
We explore how the local environment is related to properties of active galactic nuclei (AGNs) of various luminosities. Recent simulations and observations are converging on the view that the extreme luminosity of quasars, the brightest of AGNs, is fueled in major mergers of gas-rich galaxies. In such a picture, quasars, the highest luminosity AGNs, are expected to be located in regions with a h…
▽ More
We explore how the local environment is related to properties of active galactic nuclei (AGNs) of various luminosities. Recent simulations and observations are converging on the view that the extreme luminosity of quasars, the brightest of AGNs, is fueled in major mergers of gas-rich galaxies. In such a picture, quasars, the highest luminosity AGNs, are expected to be located in regions with a higher density of galaxies on small scales where mergers are more likely to take place. However, in this picture, the activity observed in low-luminosity AGNs is due to secular processes that are less dependent on the local galaxy density. To test this hypothesis, we compare the local photometric galaxy density on kiloparsec scales around spectroscopic type I and type II quasars to the local density around lower-luminosity spectroscopic type I and type II AGNs. To minimize projection effects and evolution in the photometric galaxy sample we use to characterize AGN environments, we place our random control sample at the same redshift as our AGNs and impose a narrow redshift window around both the AGNs and control targets. Our results support these merger models for bright AGN origins. We find that the brightest sources have overdensities that increase on the smallest scales compared to dimmer sources. In addition, we investigate the nature of the quasar and AGN environments themselves and find that the increased overdensity of early-type galaxies in the environments of bright type I sources suggests that they are located in richer cluster environments than dim sources. We measure increased environment overdensity with increased quasar black hole mass, consistent with the well-known M_{DMH}-M_{BH} relationship, and find evidence for quenching in the environments of high accretion efficiency type I quasars.
△ Less
Submitted 14 October, 2009;
originally announced October 2009.
-
Robust Machine Learning Applied to Astronomical Datasets III: Probabilistic Photometric Redshifts for Galaxies and Quasars in the SDSS and GALEX
Authors:
Nicholas M. Ball,
Robert J. Brunner,
Adam D. Myers,
Natalie E. Strand,
Stacey L. Alberts,
David Tcheng
Abstract:
We apply machine learning in the form of a nearest neighbor instance-based algorithm (NN) to generate full photometric redshift probability density functions (PDFs) for objects in the Fifth Data Release of the Sloan Digital Sky Survey (SDSS DR5). We use a conceptually simple but novel application of NN to generate the PDFs - perturbing the object colors by their measurement error - and using the…
▽ More
We apply machine learning in the form of a nearest neighbor instance-based algorithm (NN) to generate full photometric redshift probability density functions (PDFs) for objects in the Fifth Data Release of the Sloan Digital Sky Survey (SDSS DR5). We use a conceptually simple but novel application of NN to generate the PDFs - perturbing the object colors by their measurement error - and using the resulting instances of nearest neighbor distributions to generate numerous individual redshifts. When the redshifts are compared to existing SDSS spectroscopic data, we find that the mean value of each PDF has a dispersion between the photometric and spectroscopic redshift consistent with other machine learning techniques, being sigma = 0.0207 +/- 0.0001 for main sample galaxies to r < 17.77 mag, sigma = 0.0243 +/- 0.0002 for luminous red galaxies to r < ~19.2 mag, and sigma = 0.343 +/- 0.005 for quasars to i < 20.3 mag. The PDFs allow the selection of subsets with improved statistics. For quasars, the improvement is dramatic: for those with a single peak in their probability distribution, the dispersion is reduced from 0.343 to sigma = 0.117 +/- 0.010, and the photometric redshift is within 0.3 of the spectroscopic redshift for 99.3 +/- 0.1% of the objects. Thus, for this optical quasar sample, we can virtually eliminate 'catastrophic' photometric redshift estimates. In addition to the SDSS sample, we incorporate ultraviolet photometry from the Third Data Release of the Galaxy Evolution Explorer All-Sky Imaging Survey (GALEX AIS GR3) to create PDFs for objects seen in both surveys. For quasars, the increased coverage of the observed frame UV of the SED results in significant improvement over the full SDSS sample, with sigma = 0.234 +/- 0.010. We demonstrate that this improvement is genuine. [Abridged]
△ Less
Submitted 21 April, 2008;
originally announced April 2008.
-
AGN Environments in the Sloan Digital Sky Survey I: Dependence on Type, Redshift, and Luminosity
Authors:
Natalie E. Strand,
Robert J. Brunner,
Adam D. Myers
Abstract:
We explore how the local environment is related to the redshift, type, and luminosity of active galactic nuclei (AGN). Recent simulations and observations are converging on the view that the extreme luminosity of quasars is fueled in major mergers of gas-rich galaxies. In such a picture, quasars are expected to be located in regions with a higher density of galaxies on small scales where mergers…
▽ More
We explore how the local environment is related to the redshift, type, and luminosity of active galactic nuclei (AGN). Recent simulations and observations are converging on the view that the extreme luminosity of quasars is fueled in major mergers of gas-rich galaxies. In such a picture, quasars are expected to be located in regions with a higher density of galaxies on small scales where mergers are more likely to take place. However, in this picture, the activity observed in low-luminosity AGN is due to secular processes that are less dependent on the local galaxy density. To test this hypothesis, we compare the local photometric galaxy density on kiloparsec scales around spectroscopic Type I and Type II quasars to the local density around lower luminosity spectroscopic Type I and Type II AGN. To minimize projection effects and evolution in the photometric galaxy sample we use to characterize AGN environments, we place our random control sample at the same redshift as our AGN and impose a narrow redshift window around both the AGN and control targets. We find that higher luminosity AGN have more overdense environments compared to lower luminosity AGN on all scales out to our $2\Mpchseventy$ limit. Additionally, in the range $0.3\leqslant z\leqslant 0.6$, Type II quasars have similarly overdense environments to those of bright Type I quasars on all scales out to our $2\Mpchseventy$ limit, while the environment of dimmer Type I quasars appears to be less overdense than the environment of Type II quasars. We see increased overdensity for Type II AGN compared to Type I AGN on scales out to our limit of $2\Mpchseventy$ in overlapping redshift ranges. We also detect marginal evidence for evolution in the number of galaxies within $2\Mpchseventy$ of a quasar with redshift.
△ Less
Submitted 24 July, 2008; v1 submitted 14 December, 2007;
originally announced December 2007.
-
Quasar Clustering at $25\kpch$ from a Complete Sample of Binaries
Authors:
Adam D. Myers,
Gordon T. Richards,
Robert J. Brunner,
Donald P. Schneider,
Natalie E. Strand,
Patrick B. Hall,
Jeffrey A. Blomquist,
Donald G. York
Abstract:
We present spectroscopy of binary quasar candidates selected from Data Release 4 of the Sloan Digital Sky Survey (SDSS DR4) using Kernel Density Estimation (KDE). We present 27 new sets of observations, 10 of which are binary quasars, roughly doubling the number of known $g < 21$ binaries with component separations of 3 to 6". Only 3 of 49 spectroscopically identified objects are non-quasars, co…
▽ More
We present spectroscopy of binary quasar candidates selected from Data Release 4 of the Sloan Digital Sky Survey (SDSS DR4) using Kernel Density Estimation (KDE). We present 27 new sets of observations, 10 of which are binary quasars, roughly doubling the number of known $g < 21$ binaries with component separations of 3 to 6". Only 3 of 49 spectroscopically identified objects are non-quasars, confirming that the quasar selection efficiency of the KDE technique is $\sim95$%. Several of our observed binaries are wide-separation lens candidates that merit additional higher-resolution observations. One interesting pair may be an M star binary, or an M star-binary quasar superposition. Our candidates are initially selected by UV-excess ($u-g < 1$), but are otherwise selected irrespective of the relative colors of the quasar pair, and we thus use them to suggest optimal color similarity and photometric redshift approaches for targeting binary quasars, or projected quasar pairs. From a sample that is complete on proper scales of $23.7 < R_{prop} < 29.7\kpch$, we determine the projected quasar correlation function to be $W_p=24.0 \pm^{16.9}_{10.8}$, which is $2σ$ lower than recent estimates. We argue that our low $W_p$ estimates may indicate redshift evolution in the quasar correlation function from $z\sim1.9$ to $z\sim1.4$ on scales of $R_{prop} \sim25\kpch$. The size of this evolution broadly tracks quasar clustering on larger scales, consistent with merger-driven models of quasar origin. Although our sample alone is insufficient to detect evolution in quasar clustering on small scales, an $i$-selected DR6 KDE quasar catalog, which will contain several hundred $z \leqsim 5$ binary quasars, could easily constrain any clustering evolution at $R_{prop} \sim25\kpch$.
△ Less
Submitted 21 September, 2007;
originally announced September 2007.
-
Robust Machine Learning Applied to Astronomical Datasets II: Quantifying Photometric Redshifts for Quasars Using Instance-Based Learning
Authors:
Nicholas M. Ball,
Robert J. Brunner,
Adam D. Myers,
Natalie E. Strand,
Stacey L. Alberts,
David Tcheng,
Xavier Llorà
Abstract:
We apply instance-based machine learning in the form of a k-nearest neighbor algorithm to the task of estimating photometric redshifts for 55,746 objects spectroscopically classified as quasars in the Fifth Data Release of the Sloan Digital Sky Survey. We compare the results obtained to those from an empirical color-redshift relation (CZR). In contrast to previously published results using CZRs,…
▽ More
We apply instance-based machine learning in the form of a k-nearest neighbor algorithm to the task of estimating photometric redshifts for 55,746 objects spectroscopically classified as quasars in the Fifth Data Release of the Sloan Digital Sky Survey. We compare the results obtained to those from an empirical color-redshift relation (CZR). In contrast to previously published results using CZRs, we find that the instance-based photometric redshifts are assigned with no regions of catastrophic failure. Remaining outliers are simply scattered about the ideal relation, in a similar manner to the pattern seen in the optical for normal galaxies at redshifts z < ~1. The instance-based algorithm is trained on a representative sample of the data and pseudo-blind-tested on the remaining unseen data. The variance between the photometric and spectroscopic redshifts is sigma^2 = 0.123 +/- 0.002 (compared to sigma^2 = 0.265 +/- 0.006 for the CZR), and 54.9 +/- 0.7%, 73.3 +/- 0.6%, and 80.7 +/- 0.3% of the objects are within delta z < 0.1, 0.2, and 0.3 respectively. We also match our sample to the Second Data Release of the Galaxy Evolution Explorer legacy data and the resulting 7,642 objects show a further improvement, giving a variance of sigma^2 = 0.054 +/- 0.005, and 70.8 +/- 1.2%, 85.8 +/- 1.0%, and 90.8 +/- 0.7% of objects within delta z < 0.1, 0.2, and 0.3. We show that the improvement is indeed due to the extra information provided by GALEX, by training on the same dataset using purely SDSS photometry, which has a variance of sigma^2 = 0.090 +/- 0.007. Each set of results represents a realistic standard for application to further datasets for which the spectra are representative.
△ Less
Submitted 22 March, 2007; v1 submitted 17 December, 2006;
originally announced December 2006.