Showing 1–46 of 46 results for author: Bowles, R

Search for CP-violating Neutrino Non-Standard Interactions with the NOvA Experiment

Authors: NOvA Collaboration, M. A. Acero, B. Acharya, P. Adamson, L. Aliaga, N. Anfimov, A. Antoshkin, E. Arrieta-Diaz, L. Asquith, A. Aurisano, A. Back, N. Balashov, P. Baldi, B. A. Bambah, A. Bat, K. Bays, R. Bernstein, T. J. C. Bezerra, V. Bhatnagar, D. Bhattarai, B. Bhuyan, J. Bian, A. C. Booth, R. Bowles, B. Brahma , et al. (182 additional authors not shown)

Abstract: This Letter reports a search for charge-parity (CP) symmetry violating non-standard interactions (NSI) of neutrinos with matter using the NOvA Experiment, and examines their effects on the determination of the standard oscillation parameters. Data from $ν_μ(\barν_μ)\rightarrowν_μ(\barν_μ)$ and $ν_μ(\barν_μ)\rightarrowν_{e}(\barν_{e})$ oscillation channels are used to measure the effect of the NSI… ▽ More This Letter reports a search for charge-parity (CP) symmetry violating non-standard interactions (NSI) of neutrinos with matter using the NOvA Experiment, and examines their effects on the determination of the standard oscillation parameters. Data from $ν_μ(\barν_μ)\rightarrowν_μ(\barν_μ)$ and $ν_μ(\barν_μ)\rightarrowν_{e}(\barν_{e})$ oscillation channels are used to measure the effect of the NSI parameters $\varepsilon_{eμ}$ and $\varepsilon_{eτ}$. With 90% C.L. the magnitudes of the NSI couplings are constrained to be $|\varepsilon_{eμ}| \, \lesssim 0.3$ and $|\varepsilon_{eτ}| \, \lesssim 0.4$. A degeneracy at $|\varepsilon_{eτ}| \, \approx 1.8$ is reported, and we observe that the presence of NSI limits sensitivity to the standard CP phase $δ_{\tiny\text{CP}}$. △ Less

Submitted 11 March, 2024; originally announced March 2024.

Report number: FERMILAB-PUB-24-0108-PPD

Expanding neutrino oscillation parameter measurements in NOvA using a Bayesian approach

Authors: NOvA Collaboration, M. A. Acero, B. Acharya, P. Adamson, N. Anfimov, A. Antoshkin, E. Arrieta-Diaz, L. Asquith, A. Aurisano, A. Back, N. Balashov, P. Baldi, B. A. Bambah, A. Bat, K. Bays, R. Bernstein, T. J. C. Bezerra, V. Bhatnagar, D. Bhattarai, B. Bhuyan, J. Bian, A. C. Booth, R. Bowles, B. Brahma, C. Bromberg , et al. (174 additional authors not shown)

Abstract: NOvA is a long-baseline neutrino oscillation experiment that measures oscillations in charged-current $ν_μ \rightarrow ν_μ$ (disappearance) and $ν_μ \rightarrow ν_{e}$ (appearance) channels, and their antineutrino counterparts, using neutrinos of energies around 2 GeV over a distance of 810 km. In this work we reanalyze the dataset first examined in our previous paper [Phys. Rev. D 106, 032004 (20… ▽ More NOvA is a long-baseline neutrino oscillation experiment that measures oscillations in charged-current $ν_μ \rightarrow ν_μ$ (disappearance) and $ν_μ \rightarrow ν_{e}$ (appearance) channels, and their antineutrino counterparts, using neutrinos of energies around 2 GeV over a distance of 810 km. In this work we reanalyze the dataset first examined in our previous paper [Phys. Rev. D 106, 032004 (2022)] using an alternative statistical approach based on Bayesian Markov Chain Monte Carlo. We measure oscillation parameters consistent with the previous results. We also extend our inferences to include the first NOvA measurements of the reactor mixing angle $θ_{13}$ and the Jarlskog invariant. We use these results to quantify the strength of our inferences about CP violation, as well as to examine the effects of constraints from short-baseline measurements of $θ_{13}$ using antineutrinos from nuclear reactors when making NOvA measurements of $θ_{23}$. Our long-baseline measurement of $θ_{13}$ is also shown to be consistent with the reactor measurements, supporting the general applicability and robustness of the PMNS framework for neutrino oscillations. △ Less

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

Comments: 20 pages, 17 figures; version accepted by Phys. Rev. D. Data associated with this paper is available at https://doi.org/10.15484/2349444

Report number: FERMILAB-PUB-23-667-AD-CSAID-ND

Journal ref: Phys.Rev.D 110 (2024) 1, 012005

A topologically stabilized metastable fluid in a system of cylindrically confined hard spheres

Authors: Mahdi Zarif, Richard K. Bowles

Abstract: Metastability in soft condensed matter systems usually results from the presence of a nucleation free energy barrier and/or slow dynamics caused by high density jamming phenomena. Here, we use molecular dynamics and Monte Carlo simulation to show that the interactions between topological defects stabilizes a chiral helical fluid in a confined quasi-one-dimensional hard sphere fluid dramatically sl… ▽ More Metastability in soft condensed matter systems usually results from the presence of a nucleation free energy barrier and/or slow dynamics caused by high density jamming phenomena. Here, we use molecular dynamics and Monte Carlo simulation to show that the interactions between topological defects stabilizes a chiral helical fluid in a confined quasi-one-dimensional hard sphere fluid dramatically slowing its decay toward the equilibrium achiral fluid state. Analysis of thermodynamic, structural and dynamic properties of the system show the equation of state bifurcates continuously at intermediate pressures into two distinct branches that are accessed from different initial conditions, but terminate at the same close packed single helix in the high pressure limit. The equilibrium fluid, which forms the high pressure branch as the system is compressed from low density, is characterized by helical sections separated by randomly distributed topological defects that change the helical twist direction giving rise to an achiral fluid. The low pressure metastable branch, formed by decompressing the system from the perfect helix, is characterized by the appearance of loosely paired defects that help retain the chiral excess of the original state and stabliize the fluids until it merges continuously with the equilibrium branch at intermediate pressures. △ Less

Submitted 7 June, 2023; originally announced June 2023.

Comments: 8 Figures, 2 figures in supplemental information

Crystal Growth Rates from Molecular Liquids: The Kinetics of Entropy Loss

Authors: Richard K Bowles, Peter Harrowell

Abstract: It has been established empirically that the rate of addition of molecules to the crystal during crystal growth from the melt is proportional to exp(-|ΔS_fus|/R) where ΔS_fus is the entropy of fusion. Here we show that this entropic slowdown arises directly from the separation of the entropy loss and energy loss processes associated with the freezing of the liquid. We present a theoretical treatme… ▽ More It has been established empirically that the rate of addition of molecules to the crystal during crystal growth from the melt is proportional to exp(-|ΔS_fus|/R) where ΔS_fus is the entropy of fusion. Here we show that this entropic slowdown arises directly from the separation of the entropy loss and energy loss processes associated with the freezing of the liquid. We present a theoretical treatment of the kinetics based on a model flat energy landscape and derive an explicit expression for the coupling magnitude in terms of the crystal-melt interfacial free energy. The implications of our work for nucleation kinetics are also discussed. △ Less

Submitted 15 April, 2023; originally announced April 2023.

The Profiled Feldman-Cousins technique for confidence interval construction in the presence of nuisance parameters

Authors: M. A. Acero, B. Acharya, P. Adamson, L. Aliaga, N. Anfimov, A. Antoshkin, E. Arrieta-Diaz, L. Asquith, A. Aurisano, A. Back, C. Backhouse, M. Baird, N. Balashov, P. Baldi, B. A. Bambah, S. Bashar, A. Bat, K. Bays, R. Bernstein, V. Bhatnagar, D. Bhattarai, B. Bhuyan, J. Bian, A. C. Booth, R. Bowles , et al. (196 additional authors not shown)

Abstract: Measuring observables to constrain models using maximum-likelihood estimation is fundamental to many physics experiments. The Profiled Feldman-Cousins method described here is a potential solution to common challenges faced in constructing accurate confidence intervals: small datasets, bounded parameters, and the need to properly handle nuisance parameters. This method achieves more accurate frequ… ▽ More Measuring observables to constrain models using maximum-likelihood estimation is fundamental to many physics experiments. The Profiled Feldman-Cousins method described here is a potential solution to common challenges faced in constructing accurate confidence intervals: small datasets, bounded parameters, and the need to properly handle nuisance parameters. This method achieves more accurate frequentist coverage than other methods in use, and is generally applicable to the problem of parameter estimation in neutrino oscillations and similar measurements. We describe an implementation of this method in the context of the NOvA experiment. △ Less

Submitted 1 August, 2022; v1 submitted 28 July, 2022; originally announced July 2022.

Comments: 19 pages, 12 figures

Report number: FERMILAB-PUB-22-476-ND

Measurement of the $ν_e-$Nucleus Charged-Current Double-Differential Cross Section at $\left< E_ν \right> = $ 2.4 GeV using NOvA

Authors: M. A. Acero, P. Adamson, L. Aliaga, N. Anfimov, A. Antoshkin, E. Arrieta-Diaz, L. Asquith, A. Aurisano, A. Back, C. Backhouse, M. Baird, N. Balashov, P. Baldi, B. A. Bambah, S. Bashar, K. Bays, R. Bernstein, V. Bhatnagar, D. Bhattarai, B. Bhuyan, J. Bian, A. C. Booth, R. Bowles, B. Brahma, C. Bromberg , et al. (190 additional authors not shown)

Abstract: The inclusive electron neutrino charged-current cross section is measured in the NOvA near detector using $8.02\times10^{20}$ protons-on-target (POT) in the NuMI beam. The sample of GeV electron neutrino interactions is the largest analyzed to date and is limited by $\simeq$ 17\% systematic rather than the $\simeq$ 7.4\% statistical uncertainties. The double-differential cross section in final-sta… ▽ More The inclusive electron neutrino charged-current cross section is measured in the NOvA near detector using $8.02\times10^{20}$ protons-on-target (POT) in the NuMI beam. The sample of GeV electron neutrino interactions is the largest analyzed to date and is limited by $\simeq$ 17\% systematic rather than the $\simeq$ 7.4\% statistical uncertainties. The double-differential cross section in final-state electron energy and angle is presented for the first time, together with the single-differential dependence on $Q^{2}$ (squared four-momentum transfer) and energy, in the range 1 GeV $ \leq E_ν < $6 GeV. Detailed comparisons are made to the predictions of the GENIE, GiBUU, NEUT, and NuWro neutrino event generators. The data do not strongly favor a model over the others consistently across all three cross sections measured, though some models have especially good or poor agreement in the single differential cross section vs. $Q^{2}$. △ Less

Submitted 21 June, 2022; originally announced June 2022.

Report number: FERMILAB-PUB-22-446-ND

Is there a Granular Potential?

Authors: Josh M. Gramlich, Mahdi Zarif, Richard K. Bowles

Abstract: Granular materials, such as sand or grain, exhibit many structural and dynamic characteristics similar to those observed in molecular systems, despite temperature playing no role in their properties. This has led to an effort to develop a statistical mechanics for granular materials that has focused on establishing an equivalent to the microcanonical ensemble and a temperature-like thermodynamic v… ▽ More Granular materials, such as sand or grain, exhibit many structural and dynamic characteristics similar to those observed in molecular systems, despite temperature playing no role in their properties. This has led to an effort to develop a statistical mechanics for granular materials that has focused on establishing an equivalent to the microcanonical ensemble and a temperature-like thermodynamic variable. Here, we expand on these ideas by introducing a granular potential into the Edwards ensemble, as an analogue to the chemical potential, and explore its properties using a simple model of a granular system. A simple kinetic Monte Carlo simulation of the model shows the effect of mass transport leading to equilibrium and how this is connected to the redistribution of volume in the system. An exact analytical treatment of the model shows that the compactivity and the ratio of the granular potential to the compactivity determine the equilibrium between two open systems that are able to exchange volume and particles, and that mass moves from high to low values of this ratio. Analysis of the granular potential shows that adding a particle to the system increases the entropy at high compactivity, but decreases the entropy at low compactivity. Finally, we demonstrate the use of a small system thermodynamics method for the measurement of granular potential differences. △ Less

Submitted 17 January, 2023; v1 submitted 9 June, 2022; originally announced June 2022.

Comments: 11 Pages, 7 Figures. Extended Introduction and Discussion. New section on Measuring granular potential

The Influence on Crystal Nucleation of an Order-Disorder Transition among the Subcritical Clusters

Authors: Richard K. Bowles, Peter Harrowell

Abstract: Studies of nucleation generally focus on the properties of the critical cluster, but the presence of defects within the crystal lattice means that the population of nuclei necessarily evolve through a distribution of pre-critical clusters with varying degrees of structural disorder on their way to forming a growing stable crystal. To investigate the role pre-critical clusters play in nucleation, w… ▽ More Studies of nucleation generally focus on the properties of the critical cluster, but the presence of defects within the crystal lattice means that the population of nuclei necessarily evolve through a distribution of pre-critical clusters with varying degrees of structural disorder on their way to forming a growing stable crystal. To investigate the role pre-critical clusters play in nucleation, we develop a simple thermodynamic model for crystal nucleation in terms of cluster size and the degree of cluster order that allows us to alter the work of forming the pre-critical clusters without effecting the properties of the critical cluster. The steady state and transient nucleation behaviour of the system are then studied numerically, for different microscopic ordering kinetics. We find that the models exhibits a generic order-disorder transition in the pre-critical clusters. Independent of the type of ordering kinetics, increasing the accessibility of disordered pre-critical clusters decreases both the steady state nucleation rate and the nucleation lag time. Furthermore, the interplay between the free energy surface and the microscopic ordering kinetics leads to three distinct nucleation pathways. △ Less

Submitted 24 September, 2021; originally announced September 2021.

Comments: 9 pages. 5 Figures in main text. 3 Figures in Supplemental Information

Measurement of the Double-Differential Muon-neutrino Charged-Current Inclusive Cross Section in the NOvA Near Detector

Authors: M. A. Acero, P. Adamson, L. Aliaga, N. Anfimov, A. Antoshkin, E. Arrieta-Diaz, L. Asquith, A. Aurisano, A. Back, C. Backhouse, M. Baird, N. Balashov, P. Baldi, B. A. Bambah, S. Bashar, K. Bays, B. Behera, R. Bernstein, V. Bhatnagar, D. Bhattarai, B. Bhuyan, J. Bian, J. Blair, A. C. Booth, R. Bowles , et al. (181 additional authors not shown)

Abstract: We report cross-section measurements of the final-state muon kinematics for \numu charged-current interactions in the NOvA near detector using an accumulated 8.09$\times10^{20}$ protons-on-target (POT) in the NuMI beam. We present the results as a double-differential cross section in the observed outgoing muon energy and angle, as well as single-differential cross sections in the derived neutrino… ▽ More We report cross-section measurements of the final-state muon kinematics for \numu charged-current interactions in the NOvA near detector using an accumulated 8.09$\times10^{20}$ protons-on-target (POT) in the NuMI beam. We present the results as a double-differential cross section in the observed outgoing muon energy and angle, as well as single-differential cross sections in the derived neutrino energy, $E_ν$, and square of the four-momentum transfer, $Q^2$. We compare the results to inclusive cross-section predictions from various neutrino event generators via $χ^2$ calculations using a covariance matrix that accounts for bin-to-bin correlations of systematic uncertainties. These comparisons show a clear discrepancy between the data and each of the tested predictions at forward muon angle and low $Q^2$, indicating a missing suppression of the cross section in current neutrino-nucleus scattering models. △ Less

Submitted 18 July, 2023; v1 submitted 24 September, 2021; originally announced September 2021.

Report number: Fermilab PUB-21-455-ND-PPD-SCD

Inherent Structure Landscape of Hard Spheres Confined to Narrow Cylindrical Channels

Authors: Mahdi Zarif, Raymond J. Spiteri, Richard K. Bowles

Abstract: The inherent structure landscape for a system of hard spheres confined to a hard cylindrical channel, such that spheres can only contact their first and second neighbours, is studied using an analytical model that extends previous results [Phys. Rev. Lett. 115, 025702 (2015)] to provide a comprehensive picture of jammed packings over a range of packing densities. In the model, a packing is describ… ▽ More The inherent structure landscape for a system of hard spheres confined to a hard cylindrical channel, such that spheres can only contact their first and second neighbours, is studied using an analytical model that extends previous results [Phys. Rev. Lett. 115, 025702 (2015)] to provide a comprehensive picture of jammed packings over a range of packing densities. In the model, a packing is described as an arrangement of $k$ helical sections, separated by defects, that have alternating helical twist directions and where all spheres satisfy local jamming constraints. The structure of each helical section is determined by a single helical twist angle, and a jammed packing is obtained by minimizing the length of the channel per particle with respect to the $k$ helical section angles. An analysis of a small system of $N=20$ spheres shows that the basins on the inherent structure landscape associated with these helical arrangements split into a number of distinct jammed states separated by low barriers giving rise to a degree of hierarchical organization. The model accurately predicts the geometric properties of packings generated using the Lubachevsky and Stillinger compression scheme ($N=10^4$) and provides insight into the nature of the probability distribution of helical section lengths. △ Less

Submitted 20 August, 2021; originally announced August 2021.

Comments: 11 Figures

An Improved Measurement of Neutrino Oscillation Parameters by the NOvA Experiment

Authors: M. A. Acero, P. Adamson, L. Aliaga, N. Anfimov, A. Antoshkin, E. Arrieta-Diaz, L. Asquith, A. Aurisano, A. Back, C. Backhouse, M. Baird, N. Balashov, P. Baldi, B. A. Bambah, S. Bashar, K. Bays, R. Bernstein, V. Bhatnagar, D. Bhattarai, B. Bhuyan, J. Bian, J. Blair, A. C. Booth, R. Bowles, C. Bromberg , et al. (180 additional authors not shown)

Abstract: We present new $ν_μ\rightarrowν_e$, $ν_μ\rightarrowν_μ$, $\overlineν_μ\rightarrow\overlineν_e$, and $\overlineν_μ\rightarrow\overlineν_μ$ oscillation measurements by the NOvA experiment, with a 50% increase in neutrino-mode beam exposure over the previously reported results. The additional data, combined with previously published neutrino and antineutrino data, are all analyzed using improved tech… ▽ More We present new $ν_μ\rightarrowν_e$, $ν_μ\rightarrowν_μ$, $\overlineν_μ\rightarrow\overlineν_e$, and $\overlineν_μ\rightarrow\overlineν_μ$ oscillation measurements by the NOvA experiment, with a 50% increase in neutrino-mode beam exposure over the previously reported results. The additional data, combined with previously published neutrino and antineutrino data, are all analyzed using improved techniques and simulations. A joint fit to the $ν_e$, $ν_μ$, $\overlineν_e$, and $\overlineν_μ$ candidate samples within the 3-flavor neutrino oscillation framework continues to yield a best-fit point in the normal mass ordering and the upper octant of the $θ_{23}$ mixing angle, with $Δm^{2}_{32} = (2.41\pm0.07)\times 10^{-3}$ eV$^2$ and $\sin^2θ_{23} = 0.57^{+0.03}_{-0.04}$. The data disfavor combinations of oscillation parameters that give rise to a large asymmetry in the rates of $ν_e$ and $\overlineν_e$ appearance. This includes values of the CP-violating phase in the vicinity of $δ_\text{CP} = π/2$ which are excluded by $>3σ$ for the inverted mass ordering, and values around $δ_\text{CP} = 3π/2$ in the normal ordering which are disfavored at 2$σ$ confidence. △ Less

Submitted 8 August, 2022; v1 submitted 18 August, 2021; originally announced August 2021.

Comments: 11 pages, 6 figures. Supplementary material attached (7 figures)

Report number: FERMILAB-PUB-21-373-ND

Journal ref: Phys. Rev. D 106, 032004 (2022)

Extended search for supernova-like neutrinos in NOvA coincident with LIGO/Virgo detections

Authors: M. A. Acero, P. Adamson, L. Aliaga, N. Anfimov, A. Antoshkin, E. Arrieta-Diaz, L. Asquith, A. Aurisano, A. Back, C. Backhouse, M. Baird, N. Balashov, P. Baldi, B. A. Bambah, S. Bashar, K. Bays, R. Bernstein, V. Bhatnagar, B. Bhuyan, J. Bian, J. Blair, A. C. Booth, R. Bowles, C. Bromberg, N. Buchanan , et al. (178 additional authors not shown)

Abstract: A search is performed for supernova-like neutrino interactions coincident with 76 gravitational wave events detected by the LIGO/Virgo Collaboration. For 40 of these events, full readout of the time around the gravitational wave is available from the NOvA Far Detector. For these events, we set limits on the fluence of the sum of all neutrino flavors of $F < 7(4)\times 10^{10}\mathrm{cm}^{-2}$ at 9… ▽ More A search is performed for supernova-like neutrino interactions coincident with 76 gravitational wave events detected by the LIGO/Virgo Collaboration. For 40 of these events, full readout of the time around the gravitational wave is available from the NOvA Far Detector. For these events, we set limits on the fluence of the sum of all neutrino flavors of $F < 7(4)\times 10^{10}\mathrm{cm}^{-2}$ at 90% C.L. assuming energy and time distributions corresponding to the Garching supernova models with masses 9.6(27)$\mathrm{M}_\odot$. Under the hypothesis that any given gravitational wave event was caused by a supernova, this corresponds to a distance of $r > 29(50)$kpc at 90% C.L. Weaker limits are set for other gravitational wave events with partial Far Detector data and/or Near Detector data. △ Less

Submitted 23 August, 2021; v1 submitted 10 June, 2021; originally announced June 2021.

Comments: 10 pages, 2 figures

Report number: FERMILAB-PUB-21-276-ND

Journal ref: Phys. Rev. D 104, 063024 (2021)

Search for active-sterile antineutrino mixing using neutral-current interactions with the NOvA experiment

Authors: M. A. Acero, P. Adamson, L. Aliaga, N. Anfimov, A. Antoshkin, E. Arrieta-Diaz, L. Asquith, A. Aurisano, A. Back, C. Backhouse, M. Baird, N. Balashov, P. Baldi, B. A. Bambah, S. Bashar, K. Bays, R. Bernstein, V. Bhatnagar, B. Bhuyan, J. Bian, J. Blair, A. C. Booth, R. Bowles, C. Bromberg, N. Buchanan , et al. (174 additional authors not shown)

Abstract: This Letter reports results from the first long-baseline search for sterile antineutrinos mixing in an accelerator-based antineutrino-dominated beam. The rate of neutral-current interactions in the two NOvA detectors, at distances of 1 km and 810 km from the beam source, is analyzed using an exposure of $12.51\times10^{20}$ protons-on-target from the NuMI beam at Fermilab running in antineutrino m… ▽ More This Letter reports results from the first long-baseline search for sterile antineutrinos mixing in an accelerator-based antineutrino-dominated beam. The rate of neutral-current interactions in the two NOvA detectors, at distances of 1 km and 810 km from the beam source, is analyzed using an exposure of $12.51\times10^{20}$ protons-on-target from the NuMI beam at Fermilab running in antineutrino mode. A total of $121$ of neutral-current candidates are observed at the Far Detector, compared to a prediction of $122\pm11$(stat.)$\pm15$(syst.) assuming mixing between three active flavors. No evidence for $\barν_μ\rightarrow\barν_{s}$ oscillation is observed. Interpreting this result within a 3+1 model, constraints are placed on the mixing angles $θ_{24} < 25^{\circ}$ and $θ_{34} < 32^{\circ}$ at the 90% C.L. for $0.05$eV$^{2} \leq Δm^{2}_{41} \leq 0.5$eV$^{2}$, the range of mass splittings that produces no significant oscillations at the Near Detector. These are the first 3+1 confidence limits set using long-baseline accelerator antineutrinos. △ Less

Submitted 30 September, 2021; v1 submitted 8 June, 2021; originally announced June 2021.

Comments: 8 pages, 4 figures

Report number: FERMILAB-PUB-21-271-ND

Seasonal Variation of Multiple-Muon Cosmic Ray Air Showers Observed in the NOvA Detector on the Surface

Authors: M. A. Acero, P. Adamson, L. Aliaga, N. Anfimov, A. Antoshkin, E. Arrieta-Diaz, L. Asquith, A. Aurisano, A. Back, C. Backhouse, M. Baird, N. Balashov, P. Baldi, B. A. Bambah, S. Bashar, K. Bays, R. Bernstein, V. Bhatnagar, B. Bhuyan, J. Bian, J. Blair, A. C. Booth, R. Bowles, C. Bromberg, N. Buchanan , et al. (172 additional authors not shown)

Abstract: We report the rate of cosmic ray air showers with multiplicities exceeding 15 muon tracks recorded in the NOvA Far Detector between May 2016 and May 2018. The detector is located on the surface under an overburden of 3.6 meters water equivalent. We observe a seasonal dependence in the rate of multiple-muon showers, which varies in magnitude with multiplicity and zenith angle. During this period, t… ▽ More We report the rate of cosmic ray air showers with multiplicities exceeding 15 muon tracks recorded in the NOvA Far Detector between May 2016 and May 2018. The detector is located on the surface under an overburden of 3.6 meters water equivalent. We observe a seasonal dependence in the rate of multiple-muon showers, which varies in magnitude with multiplicity and zenith angle. During this period, the effective atmospheric temperature and surface pressure ranged between 210 K to 230 K and 940mbar to 990mbar, respectively; the shower rates are anti-correlated with the variation in the effective temperature. The variations are about 30% larger for the highest multiplicities than the lowest multiplicities and 20% larger for showers near the horizon than vertical showers. △ Less

Submitted 13 July, 2021; v1 submitted 9 May, 2021; originally announced May 2021.

Report number: FERMILAB-PUB-21-224-ND

Journal ref: Phys. Rev. D 104, 012014 (2021)

Free energy surface of two-step nucleation

Authors: Dean Eaton, Ivan Saika-Voivod, Richard K. Bowles, Peter H. Poole

Abstract: We test the theoretical free energy surface (FES) for two-step nucleation (TSN) proposed by Iwamatsu [J. Chem. Phys. 134, 164508 (2011)] by comparing the predictions of the theory to numerical results for the FES recently reported from Monte Carlo simulations of TSN in a simple lattice system [James, et al., J. Chem. Phys. 150, 074501 (2019)]. No adjustable parameters are used to make this compari… ▽ More We test the theoretical free energy surface (FES) for two-step nucleation (TSN) proposed by Iwamatsu [J. Chem. Phys. 134, 164508 (2011)] by comparing the predictions of the theory to numerical results for the FES recently reported from Monte Carlo simulations of TSN in a simple lattice system [James, et al., J. Chem. Phys. 150, 074501 (2019)]. No adjustable parameters are used to make this comparison. That is, all the parameters of the theory are evaluated directly for the model system, yielding a predicted FES which we then compare to the FES obtained from simulations. We find that the theoretical FES successfully predicts the numerically-evaluated FES over a range of thermodynamic conditions that spans distinct regimes of behavior associated with TSN. All the qualitative features of the FES are captured by the theory and the quantitative comparison is also very good. Our results demonstrate that Iwamatsu's extension of classical nucleation theory provides an excellent framework for understanding the thermodynamics of TSN. △ Less

Submitted 4 May, 2021; originally announced May 2021.

Comments: 15 pages, 15 figures

Journal ref: J. Chem. Phys. 154, 234507 (2021)

Developed liquid film passing a trailing edge: small-scale analysis and the `teapot effect' at large Reynolds numbers

Authors: Bernhard Scheichl, Robert I. Bowles, Georgios Pasias

Abstract: Recently, the authors considered a thin steady developed viscous liquid wall jet passing the sharp trailing edge of a horizontally aligned flat plate under surface tension and the weak action of gravity acting vertically in the asymptotic slender-layer limit (J. Fluid Mech. 850, pp. 924--953, 2018). We revisit the capillarity-driven short-scale viscous--inviscid interaction, on account of the inhe… ▽ More Recently, the authors considered a thin steady developed viscous liquid wall jet passing the sharp trailing edge of a horizontally aligned flat plate under surface tension and the weak action of gravity acting vertically in the asymptotic slender-layer limit (J. Fluid Mech. 850, pp. 924--953, 2018). We revisit the capillarity-driven short-scale viscous--inviscid interaction, on account of the inherent upstream influence, immediately downstream of the edge and scrutinise flow detachment on all smaller scales. We adhere to the assumption of a Froude number so large that choking at the plate edge is insignificant but envisage the variation of the relevant Weber number of O(1). The aspect in the main focus, tackled essentially analytically, is the continuation of the structure of the flow towards scales much smaller then the interactive ones and where it no longer can be treated as slender. As a remarkable phenomenon, this analysis predicts harmonic capillary ripples of Rayleigh type, prevalent on the free surface upstream of the trailing edge. They exhibit an increase of both the wavelength and amplitude as the characteristic Weber number decreases. Finally, the theory clarifies the actual detachment process, unprecedented in the rational theories of flow separation. At this stage, the wetting properties of the fluid and the microscopically wedge-shaped edge, viewed as infinitely thin on the larger scales, come into play. As this geometry typically models the exit of a spout, the predicted wetting of the wedge is related to what in the literature is referred to as the so-called teapot effect. △ Less

Submitted 18 May, 2021; v1 submitted 24 November, 2020; originally announced November 2020.

Comments: 45 pages, 12 figures

MSC Class: 76A20

Search for Slow Magnetic Monopoles with the NOvA Detector on the Surface

Authors: NOvA Collaboration, M. A. Acero, P. Adamson, L. Aliaga, T. Alion, V. Allakhverdian, N. Anfimov, A. Antoshkin, E. Arrieta-Diaz, L. Asquith, A. Aurisano, A. Back, C. Backhouse, M. Baird, N. Balashov, P. Baldi, B. A. Bambah, S. Bashar, K. Bays, S. Bending, R. Bernstein, V. Bhatnagar, B. Bhuyan, J. Bian, J. Blair , et al. (174 additional authors not shown)

Abstract: We report a search for a magnetic monopole component of the cosmic-ray flux in a 95-day exposure of the NOvA experiment's Far Detector, a 14 kt segmented liquid scintillator detector designed primarily to observe GeV-scale electron neutrinos. No events consistent with monopoles were observed, setting an upper limit on the flux of $2\times 10^{-14} \mathrm{cm^{-2}s^{-1}sr^{-1}}$ at 90% C.L. for mon… ▽ More We report a search for a magnetic monopole component of the cosmic-ray flux in a 95-day exposure of the NOvA experiment's Far Detector, a 14 kt segmented liquid scintillator detector designed primarily to observe GeV-scale electron neutrinos. No events consistent with monopoles were observed, setting an upper limit on the flux of $2\times 10^{-14} \mathrm{cm^{-2}s^{-1}sr^{-1}}$ at 90% C.L. for monopole speed $6\times 10^{-4} < β< 5\times 10^{-3}$ and mass greater than $5\times 10^{8}$ GeV. Because of NOvA's small overburden of 3 meters-water equivalent, this constraint covers a previously unexplored low-mass region. △ Less

Submitted 5 January, 2021; v1 submitted 10 September, 2020; originally announced September 2020.

Comments: 8 pages, 7 figures

Report number: FERMILAB-PUB-20-472-ND

Journal ref: Phys. Rev. D 103, 012007 (2021)

Adjusting Neutrino Interaction Models and Evaluating Uncertainties using NOvA Near Detector Data

Authors: NOvA Collaboration, M. A. Acero, P. Adamson, G. Agam, L. Aliaga, T. Alion, V. Allakhverdian, N. Anfimov, A. Antoshkin, L. Asquith, A. Aurisano, A. Back, C. Backhouse, M. Baird, N. Balashov, P. Baldi, B. A. Bambah, S. Bashar, K. Bays, S. Bending, R. Bernstein, V. Bhatnagar, B. Bhuyan, J. Bian, J. Blair , et al. (170 additional authors not shown)

Abstract: The two-detector design of the NOvA neutrino oscillation experiment, in which two functionally identical detectors are exposed to an intense neutrino beam, aids in canceling leading order effects of cross-section uncertainties. However, limited knowledge of neutrino interaction cross sections still gives rise to some of the largest systematic uncertainties in current oscillation measurements. We s… ▽ More The two-detector design of the NOvA neutrino oscillation experiment, in which two functionally identical detectors are exposed to an intense neutrino beam, aids in canceling leading order effects of cross-section uncertainties. However, limited knowledge of neutrino interaction cross sections still gives rise to some of the largest systematic uncertainties in current oscillation measurements. We show contemporary models of neutrino interactions to be discrepant with data from NOvA, consistent with discrepancies seen in other experiments. Adjustments to neutrino interaction models in GENIE that improve agreement with our data are presented. We also describe systematic uncertainties on these models, including uncertainties on multi-nucleon interactions from a newly developed procedure using NOvA near detector data. △ Less

Submitted 10 December, 2020; v1 submitted 15 June, 2020; originally announced June 2020.

Comments: Code implementing adjustments to GENIE 2.12.2 described in this paper is available at https://github.com/novaexperiment/NOvARwgt-public

Report number: FERMILAB-PUB-20-243-ND

Journal ref: Eur. Phys. J. C 80, 1119 (2020)

Supernova neutrino detection in NOvA

Authors: NOvA Collaboration, M. A. Acero, P. Adamson, G. Agam, L. Aliaga, T. Alion, V. Allakhverdian, N. Anfimov, A. Antoshkin, E. Arrieta-Diaz, L. Asquith, A. Aurisano, A. Back, C. Backhouse, M. Baird, N. Balashov, P. Baldi, B. A. Bambah, S. Bashar, K. Bays, S. Bending, R. Bernstein, V. Bhatnagar, B. Bhuyan, J. Bian , et al. (177 additional authors not shown)

Abstract: The NOvA long-baseline neutrino experiment uses a pair of large, segmented, liquid-scintillator calorimeters to study neutrino oscillations, using GeV-scale neutrinos from the Fermilab NuMI beam. These detectors are also sensitive to the flux of neutrinos which are emitted during a core-collapse supernova through inverse beta decay interactions on carbon at energies of… ▽ More The NOvA long-baseline neutrino experiment uses a pair of large, segmented, liquid-scintillator calorimeters to study neutrino oscillations, using GeV-scale neutrinos from the Fermilab NuMI beam. These detectors are also sensitive to the flux of neutrinos which are emitted during a core-collapse supernova through inverse beta decay interactions on carbon at energies of $\mathcal{O}(10~\text{MeV})$. This signature provides a means to study the dominant mode of energy release for a core-collapse supernova occurring in our galaxy. We describe the data-driven software trigger system developed and employed by the NOvA experiment to identify and record neutrino data from nearby galactic supernovae. This technique has been used by NOvA to self-trigger on potential core-collapse supernovae in our galaxy, with an estimated sensitivity reaching out to 10~kpc distance while achieving a detection efficiency of 23\% to 49\% for supernovae from progenitor stars with masses of 9.6M$_\odot$ to 27M$_\odot$, respectively. △ Less

Submitted 29 July, 2020; v1 submitted 14 May, 2020; originally announced May 2020.

Comments: 30 pages, 17 figures

Report number: FERMILAB-PUB-20-201-E

Journal ref: JCAP 10 (2020) 014

Mapping diffusivity of narrow channels into one-dimension

Authors: Mahdi Zarif, Richard K. Bowles

Abstract: The diffusion of particles trapped in long narrow channels occurs predominantly in one dimension. Here, molecular dynamics simulation is used to study the inertial dynamics of two-dimensional hard disks, confined to long, narrow, structureless channels with hard walls in the no-passing regime. We show that the diffusion coefficient obtained from the mean squared displacement can be mapped onto the… ▽ More The diffusion of particles trapped in long narrow channels occurs predominantly in one dimension. Here, molecular dynamics simulation is used to study the inertial dynamics of two-dimensional hard disks, confined to long, narrow, structureless channels with hard walls in the no-passing regime. We show that the diffusion coefficient obtained from the mean squared displacement can be mapped onto the exact results for the diffusion of the strictly one-dimensional hard rod system through an effective occupied volume fraction obtained from either the equation of state or a geometric projection of the particle interaction diameters. △ Less

Submitted 13 January, 2020; originally announced January 2020.

Comments: 7 pages, 10 figures

Journal ref: Phys. Rev. E 101, 012908 (2020)

The effect of soft repulsive interactions on the diffusion of particles in quasi-one-dimensional channels: A hopping time approach

Authors: Sheida Ahmadi, Marina Schmidt, Raymond J. Spiteri, Richard K. Bowles

Abstract: Fluids confined to quasi-one-dimensional channels exhibit a dynamic crossover from single file diffusion to normal diffusion as the channel becomes wide enough for particles to hop past each other. In the crossover regime, where hopping events are rare, the diffusion coefficient in the long time limit can be related to a hopping time that measures the average time it takes a particle to escape the… ▽ More Fluids confined to quasi-one-dimensional channels exhibit a dynamic crossover from single file diffusion to normal diffusion as the channel becomes wide enough for particles to hop past each other. In the crossover regime, where hopping events are rare, the diffusion coefficient in the long time limit can be related to a hopping time that measures the average time it takes a particle to escape the local cage formed by its neighbours. In this work, we show that a transition state theory that calculates the free energy barrier for two particles attempting to pass each other in the small system isobaric ensemble is able to quantitatively predict the hopping time in a system of two-dimensional soft repulsive discs [$U(r_{ij})=(σ/r_{ij})^α$] confined to a hard walled channel over a range of channel radii and degrees of particle softness measured in terms of $1/α$. The free energy barrier exhibits a maximum at intermediate values of $α$ that moves to smaller values of $1/α$ (harder particles) as the channel becomes narrower. However, the presence of the maximum is only observed in the hopping times for wide channels because the interaction potential dependence of the kinetic prefactor plays an increasingly important role for narrower channels. We also begin to explore how our transition state theory approach can be used to optimize and control dynamics in confined quasi-one-dimensional fluids. △ Less

Submitted 19 April, 2019; originally announced April 2019.

Comments: 9 pages, 10 figures, 1 table

Phase transitions in fluctuations and their role in two-step nucleation

Authors: Daniella James, Seamus Beairsto, Carmen Hartt, Oleksandr Zavalov, Ivan Saika-Voivod, Richard K. Bowles, Peter H. Poole

Abstract: We consider the thermodynamic behavior of local fluctuations occurring in a stable or metastable bulk phase. For a system with three or more phases, a simple analysis based on classical nucleation theory predicts that small fluctuations always resemble the phase having the lowest surface tension with the surrounding bulk phase, regardless of the relative chemical potentials of the phases. We also… ▽ More We consider the thermodynamic behavior of local fluctuations occurring in a stable or metastable bulk phase. For a system with three or more phases, a simple analysis based on classical nucleation theory predicts that small fluctuations always resemble the phase having the lowest surface tension with the surrounding bulk phase, regardless of the relative chemical potentials of the phases. We also identify the conditions at which a fluctuation may convert to a different phase as its size increases, referred to here as a "fluctuation phase transition" (FPT). We demonstrate these phenonena in simulations of a two dimensional lattice model by evaluating the free energy surface that describes the thermodynamic properties of a fluctuation as a function of its size and phase composition. We show that a FPT can occur in the fluctuations of either a stable or metastable bulk phase and that the transition is first-order. We also find that the FPT is bracketed by well-defined spinodals, which place limits on the size of fluctuations of distinct phases. Furthermore, when the FPT occurs in a metastable bulk phase, we show that the superposition of the FPT on the nucleation process results in two-step nucleation (TSN). We identify distinct regimes of TSN based on the nucleation pathway in the free energy surface, and correlate these regimes to the phase diagram of the bulk system. Our results clarify the origin of TSN, and elucidate a wide variety of phenomena associated with TSN, including the Ostwald step rule. △ Less

Submitted 5 September, 2018; originally announced September 2018.

Comments: 35 pages, 12 figures

Journal ref: J. Chem. Phys. 150, 074501 (2019)

"Swarm relaxation": Equilibrating a large ensemble of computer simulations

Authors: Shahrazad M. A. Malek, Richard K. Bowles, Ivan Saika-Voivod, Francesco Sciortino, Peter H. Poole

Abstract: It is common practice in molecular dynamics and Monte Carlo computer simulations to run multiple, separately-initialized simulations in order to improve the sampling of independent microstates. Here we examine the utility of an extreme case of this strategy, in which we run a large ensemble of $M$ independent simulations (a "swarm"), each of which is relaxed to equilibrium. We show that if $M$ is… ▽ More It is common practice in molecular dynamics and Monte Carlo computer simulations to run multiple, separately-initialized simulations in order to improve the sampling of independent microstates. Here we examine the utility of an extreme case of this strategy, in which we run a large ensemble of $M$ independent simulations (a "swarm"), each of which is relaxed to equilibrium. We show that if $M$ is of order $10^3$, we can monitor the swarm's relaxation to equilibrium, and confirm its attainment, within $\sim 10\barτ$, where $\barτ$ is the equilibrium relaxation time. As soon as a swarm of this size attains equilibrium, the ensemble of $M$ final microstates from each run is sufficient for the evaluation of most equilibrium properties without further sampling. This approach dramatically reduces the wall-clock time required, compared to a single long simulation, by a factor of several hundred, at the cost of an increase in the total computational effort by a small factor. It is also well-suited to modern computing systems having thousands of processors, and is a viable strategy for simulation studies that need to produce high-precision results in a minimum of wall-clock time. We present results obtained by applying this approach to several test cases. △ Less

Submitted 29 October, 2017; originally announced October 2017.

Comments: 12 pages. To appear in Eur. Phy. J. E, 2017

Journal ref: Eur. Phys. J. E 40, 98 (2017)

Vapor Condensed and Supercooled Glassy Nanoclusters

Authors: Weikai Qi, Richard K. Bowles

Abstract: We use molecular simulation to study the structural and dynamic properties of glassy nanoclusters formed both through the direct condensation of the vapor below the glass transition temperature, without the presence of a substrate, and \textit{via} the slow supercooling of unsupported liquid nanodroplets. An analysis of local structure using Voronoi polyhedra shows that the energetic stability of… ▽ More We use molecular simulation to study the structural and dynamic properties of glassy nanoclusters formed both through the direct condensation of the vapor below the glass transition temperature, without the presence of a substrate, and \textit{via} the slow supercooling of unsupported liquid nanodroplets. An analysis of local structure using Voronoi polyhedra shows that the energetic stability of the clusters is characterized by a large, increasing fraction of bicapped square antiprism motifs. We also show that nanoclusters with similar inherent structure energies are structurally similar, independent of their history, which suggests the supercooled clusters access the same low energy regions of the potential energy landscape as the vapor condensed clusters despite their different methods of formation. By measuring the intermediate scattering function at different radii from the cluster center, we find that the relaxation dynamics of the clusters are inhomogeneous, with the core becoming glassy above the glass transition temperature while the surface remains mobile at low temperatures. This helps the clusters sample the highly stable, low energy structures on the potential energy surface. Our work suggests the nanocluster systems are structurally more stable than the ultra-stable glassy thin films, formed through vapor deposition onto a cold substrate, but the nanoclusters do not exhibit the superheating effects characteristic of the ultra-stable glass states. △ Less

Submitted 16 February, 2017; originally announced February 2017.

Comments: 4 Figures

Journal ref: ACS Nano 10, 3416 (2016)

Diffusion in Quasi-One Dimensional Channels: A Small System $n,p,T,$ Transition State Theory for Hopping Times

Authors: Sheida Ahmadi, Richard K. Bowles

Abstract: Particles confined to a single file, in a narrow quasi-one dimensional channel, exhibit a dynamic crossover from single file diffusion to Fickian diffusion as the channel radius increases and the particles can begin to pass each other. The long time diffusion coefficient for a system in the crossover regime can be described in terms of a hopping time, which measures the time it takes for a particl… ▽ More Particles confined to a single file, in a narrow quasi-one dimensional channel, exhibit a dynamic crossover from single file diffusion to Fickian diffusion as the channel radius increases and the particles can begin to pass each other. The long time diffusion coefficient for a system in the crossover regime can be described in terms of a hopping time, which measures the time it takes for a particle to escape the cage formed by its neighbours. In this paper, we develop a transition state theory approach to the calculation of the hopping time, using the small system isobaric--isothermal ensemble to rigorously account for the volume fluctuations associated with the size of the cage. We also describe a Monte Carlo simulation scheme that can be used to calculate the free energy barrier for particle hopping. The theory and simulation method correctly predict the hopping times for a two dimensional confined ideal gas system and a system of confined hard discs over a range of channel radii, but the method breaks down for wide channels in the hard discs case, underestimating the height of the hopping barrier due to the neglect of interactions between the small system and its surroundings. △ Less

Submitted 26 April, 2017; v1 submitted 24 January, 2017; originally announced January 2017.

Comments: 13 Pages and & 7 Figures: Added Discussion corrected typos

Journal ref: Journal of Chemical Physics, 146, 154505 (2017)

Competitive Heterogeneous Nucleation onto a Microscopic Impurity in a Potts Model

Authors: Cletus C. Asuquo, Danielle McArthur, Richard K. Bowles

Abstract: Many metastable systems can nucleate to multiple competing stable or intermediate metastable states. In this work, a Potts model, subject to external fields, is used to study the competitive nucleation of two phases attempting to grow on a microscopic impurity. Monte Carlo simulations are used to calculate the free energy surfaces for the system under different conditions, where the relative stabi… ▽ More Many metastable systems can nucleate to multiple competing stable or intermediate metastable states. In this work, a Potts model, subject to external fields, is used to study the competitive nucleation of two phases attempting to grow on a microscopic impurity. Monte Carlo simulations are used to calculate the free energy surfaces for the system under different conditions, where the relative stability of the phases is adjusted by changing the interaction parameters, and the nucleation rates obtained using multicomponent transition state theory are compared with the rates measured using the survival probability method. We find that the two methods predict similar nucleation rates when the free energy barrier used in the transition state theory is defined as the work required to form a critical embryo from the metastable phase. An analysis of the free energy surfaces also reveals that the competition between the nucleating phases leads to an effective drying of the impurity which slows down the nucleation rate compared to the single phase case. △ Less

Submitted 17 August, 2016; originally announced August 2016.

Comments: 8 pages, 7 Figures

Journal ref: J. Chem. Phys. 145, 064511 (2016)

Free energy of formation of small ice nuclei near the Widom line in simulations of supercooled water

Authors: Connor R. C. Buhariwalla, Richard K. Bowles, Ivan Saika-Voivod, Francesco Sciortino, Peter H. Poole

Abstract: The ST2 interaction potential has been used in a large number of simulation studies to explore the possibility of a liquid-liquid phase transition (LLPT) in supercooled water. Using umbrella sampling Monte Carlo simulations of ST2 water, we evaluate the free energy of formation of small ice nuclei in the supercooled liquid in the vicinity of the Widom line, the region above the critical temperatur… ▽ More The ST2 interaction potential has been used in a large number of simulation studies to explore the possibility of a liquid-liquid phase transition (LLPT) in supercooled water. Using umbrella sampling Monte Carlo simulations of ST2 water, we evaluate the free energy of formation of small ice nuclei in the supercooled liquid in the vicinity of the Widom line, the region above the critical temperature of the LLPT where a number of thermodynamic anomalies occur. Our results show that in this region there is a substantial free energy cost for the formation of small ice nuclei, demonstrating that the thermodynamic anomalies associated with the Widom line in ST2 water occur in a well-defined metastable liquid phase. On passing through the Widom line, we identify changes in the free energy to form small ice nuclei that illustrate how the thermodynamic anomalies associated with the LLPT may influence the ice nucleation process. △ Less

Submitted 13 January, 2015; originally announced January 2015.

Comments: 11 pages, 15 figures

Journal ref: Eur. Phys. J. E 38, 39 (2015)

An orientationally ordered helical fluid phase in a quasi-one-dimensional system of confined hard spheres

Authors: Mahdi Zaeifi Yamchi, Richard K. Bowles

Abstract: We use a series of molecular dynamics simulations, and analytical theory, to demonstrate that a system of hard spheres confined to a narrow cylindrical channel exhibits a continuous phase transition from an isotropic fluid at low densities, to an orientationally ordered, but translationally disordered, helical fluid at high densities. The ordered fluid phase contains small sections of helix separa… ▽ More We use a series of molecular dynamics simulations, and analytical theory, to demonstrate that a system of hard spheres confined to a narrow cylindrical channel exhibits a continuous phase transition from an isotropic fluid at low densities, to an orientationally ordered, but translationally disordered, helical fluid at high densities. The ordered fluid phase contains small sections of helix separated by topological defects that change the direction of the twist, altering the local chirailty. The defects break up the translational order, but the fluid develops long range orientational order. An analysis of the particle packings show that the length separation between defects controls the geometrical properties of the helical sections, including the orientation, and that pairs of defects experience a weak, but long range attraction resulting from entropic free volumes effects. These collective long range interactions overcome the restrictions on quasi-one-dimensional transitions, even though the particle-particle interaction is short ranged. △ Less

Submitted 16 September, 2014; originally announced September 2014.

Comments: 5 pages, 6 figs

MIGRAINE: MRI Graph Reliability Analysis and Inference for Connectomics

Authors: William Gray Roncal, Zachary H. Koterba, Disa Mhembere, Dean M. Kleissas, Joshua T. Vogelstein, Randal Burns, Anita R. Bowles, Dimitrios K. Donavos, Sephira Ryman, Rex E. Jung, Lei Wu, Vince Calhoun, R. Jacob Vogelstein

Abstract: Currently, connectomes (e.g., functional or structural brain graphs) can be estimated in humans at $\approx 1~mm^3$ scale using a combination of diffusion weighted magnetic resonance imaging, functional magnetic resonance imaging and structural magnetic resonance imaging scans. This manuscript summarizes a novel, scalable implementation of open-source algorithms to rapidly estimate magnetic resona… ▽ More Currently, connectomes (e.g., functional or structural brain graphs) can be estimated in humans at $\approx 1~mm^3$ scale using a combination of diffusion weighted magnetic resonance imaging, functional magnetic resonance imaging and structural magnetic resonance imaging scans. This manuscript summarizes a novel, scalable implementation of open-source algorithms to rapidly estimate magnetic resonance connectomes, using both anatomical regions of interest (ROIs) and voxel-size vertices. To assess the reliability of our pipeline, we develop a novel nonparametric non-Euclidean reliability metric. Here we provide an overview of the methods used, demonstrate our implementation, and discuss available user extensions. We conclude with results showing the efficacy and reliability of the pipeline over previous state-of-the-art. △ Less

Submitted 17 December, 2013; originally announced December 2013.

Comments: Published as part of 2013 IEEE GlobalSIP conference

A Transition State Theory for Calculating Hopping Times and Diffusion in Highly Confined Fluids

Authors: Surajith N. Wanasundara, Raymond J. Spiteri, Richard K. Bowles

Abstract: Monte Carlo simulation is used to study the dynamical crossover from single file diffusion to normal diffusion in fluids confined to narrow channels. We show that the long time diffusion coefficients for a series of systems involving hard and soft interaction potentials can be described in terms of a hopping time that measures the time it takes for a particle to escape the cage formed by its neigh… ▽ More Monte Carlo simulation is used to study the dynamical crossover from single file diffusion to normal diffusion in fluids confined to narrow channels. We show that the long time diffusion coefficients for a series of systems involving hard and soft interaction potentials can be described in terms of a hopping time that measures the time it takes for a particle to escape the cage formed by its neighbors in the pore. Free energy barriers for the particle hopping process are calculated and used to show that transition state theory effectively describes the hopping time for all the systems studied, over a range of pore diameters. Our work suggests that the combination of hopping times and transition state theory offers a useful and general framework to describe the dynamics of these highly confined fluids. △ Less

Submitted 12 October, 2013; originally announced October 2013.

Comments: 6 figures

Vapor condensation onto a non-volatile liquid drop

Authors: Levent Inci, Richard. K. Bowles

Abstract: Molecular dynamics simulations of miscible and partially miscible binary Lennard--Jones mixtures are used to study the dynamics and thermodynamics of vapor condensation onto a non-volatile liquid drop in the canonical ensemble. When the system volume is large, the driving force for condensation is low and only a submonolayer of the solvent is adsorbed onto the liquid drop. A small degree of mixing… ▽ More Molecular dynamics simulations of miscible and partially miscible binary Lennard--Jones mixtures are used to study the dynamics and thermodynamics of vapor condensation onto a non-volatile liquid drop in the canonical ensemble. When the system volume is large, the driving force for condensation is low and only a submonolayer of the solvent is adsorbed onto the liquid drop. A small degree of mixing of the solvent phase into the core of the particles occurs for the miscible system. At smaller volumes complete film formation is observed and the dynamics of film growth are dominated by cluster-cluster coalescence. Mixing into the core of the droplet is also observed for partially miscible systems below an onset volume suggesting, the presence of a solubility transition. We also develop a non-volatile liquid drop model, based on the capillarity approximations, that exhibits a solubility transition between small and large drops for partially miscible mixtures and has a hysteresis loop similar to the one observed in the deliquescence of small soluble salt particles. The properties of the model are compared to our simulation results and the model is used to study the formulation of classical nucleation theory for systems with low free energy barriers. △ Less

Submitted 11 October, 2013; originally announced October 2013.

Comments: 12 Figures

Free energy surface of ST2 water near the liquid-liquid phase transition

Authors: Peter H. Poole, Richard K. Bowles, Ivan Saika-Voivod, Francesco Sciortino

Abstract: We carry out umbrella sampling Monte Carlo simulations to evaluate the free energy surface of the ST2 model of water as a function two order parameters, the density and a bond-orientational order parameter. We approximate the long-range electrostatic interactions of the ST2 model using the reaction-field method. We focus on state points in the vicinity of the liquid-liquid critical point proposed… ▽ More We carry out umbrella sampling Monte Carlo simulations to evaluate the free energy surface of the ST2 model of water as a function two order parameters, the density and a bond-orientational order parameter. We approximate the long-range electrostatic interactions of the ST2 model using the reaction-field method. We focus on state points in the vicinity of the liquid-liquid critical point proposed for this model in earlier work. At temperatures below the predicted critical temperature we find two basins in the free energy surface, both of which have liquid-like bond orientational order, but differing in density. The pressure and temperature dependence of the shape of the free energy surface is consistent with the assignment of these two basins to the distinct low density and high density liquid phases previously predicted to occur in ST2 water. △ Less

Submitted 7 December, 2012; originally announced December 2012.

Comments: 8 pages, 9 figures

Journal ref: J. Chem. Phys. 138, 034505 (2013)

A Fragile-Strong Fluid Crossover and Universal Relaxation Times in a Confined Hard Disc Fluid

Authors: Mahdi Zaeifi Yamchi, S. S. Ashwin, Richard K. Bowles

Abstract: We show that a system of hard discs confined to a narrow channel exhibits a fragile-strong fluid crossover located at the maximum of the isobaric heat capacity and that the relaxation times for different channel widths fall onto a single master curve when rescaled by the relaxation times and temperatures of the crossover. Calculations of the configurational entropy and the inherent structure equat… ▽ More We show that a system of hard discs confined to a narrow channel exhibits a fragile-strong fluid crossover located at the maximum of the isobaric heat capacity and that the relaxation times for different channel widths fall onto a single master curve when rescaled by the relaxation times and temperatures of the crossover. Calculations of the configurational entropy and the inherent structure equation of state find that the crossover is related to properties of the jamming landscape for the model but that the Adams-Gibbs relation does not predict the relaxation behavior. We also show that a facilitated dynamics description of the system, where kinetically excited regions are identified with local packing arrangements of the discs, successfully describes the fragile-strong crossover. △ Less

Submitted 22 October, 2012; v1 submitted 9 August, 2012; originally announced August 2012.

Comments: 5 pages 5 figs, updated text + Figs, Accepted PRL

The Inherent Structure Landscape Connection Between Liquids, Granular materials and the Jamming Phase Diagram

Authors: S. S. Ashwin, Mahdi Zaeifi Yamchi, Richard. K. Bowles

Abstract: We provide a comprehensive picture of the jamming phase diagram by connecting the athermal, granular ensemble of jammed states and the equilibrium fluid through the inherent structure paradigm for a system hard discs confined to a narrow channel. The J-line is shown to be divided into packings that are thermodynamically accessible from the equilibrium fluid and inaccessible packings. The J-point i… ▽ More We provide a comprehensive picture of the jamming phase diagram by connecting the athermal, granular ensemble of jammed states and the equilibrium fluid through the inherent structure paradigm for a system hard discs confined to a narrow channel. The J-line is shown to be divided into packings that are thermodynamically accessible from the equilibrium fluid and inaccessible packings. The J-point is found to occur at the transition between these two sets of packings and is located at the maximum the inherent structure distribution. A general thermodynamic argument suggests that the density of the states at the configurational entropy maximum represents a lower bound on the J-point density in hard sphere systems. Finally, we find that the granular and fluid systems only occupy the same set of inherent structures, under the same thermodynamic conditions, at two points, corresponding to zero and infinite pressures, where they sample the J-point states and the most dense packing respectively. △ Less

Submitted 7 August, 2012; originally announced August 2012.

Comments: 5 pages, 3 Figures

Heterogeneous Nucleation in the Low Barrier Regime

Authors: Benjamin Scheifele, Ivan Saika-Voivod, Richard K. Bowles, Peter H. Poole

Abstract: In simulations of the 2D Ising model, we examine heterogeneous nucleation induced by a small impurity consisting of a line of $l$ fixed spins. As $l$ increases, we identify a limit of stability beyond which the metastable phase is not defined. We evaluate the free energy barrier for nucleation of the stable phase, and show that contrary to expectation, the barrier does not vanish on approach to th… ▽ More In simulations of the 2D Ising model, we examine heterogeneous nucleation induced by a small impurity consisting of a line of $l$ fixed spins. As $l$ increases, we identify a limit of stability beyond which the metastable phase is not defined. We evaluate the free energy barrier for nucleation of the stable phase, and show that contrary to expectation, the barrier does not vanish on approach to the limit of stability. We also demonstrate that our values for the height of the barrier yield predictions for the nucleation time (from transition state theory) and the size of the critical cluster (from the nucleation theorem) that are in excellent agreement with direct measurements, even near the limit of stability. △ Less

Submitted 18 April, 2013; v1 submitted 5 August, 2012; originally announced August 2012.

Comments: 7 pages, 6 figures; to appear in Phys. Rev. E

Journal ref: Phys. Rev. E 87, 042407 (2013)

Heterogeneous condensation of the Lennard-Jones vapor onto a nanoscale seed particle

Authors: Levent Inci, Richard K. Bowles

Abstract: The heterogeneous condensation of a Lennard-Jones vapor onto a nanoscale seed particle is studied using molecular dynamics simulations. Measuring the nucleation rate and the height of the free energy barrier using the mean first passage time method shows that the presence of a weakly interacting seed has little effect on the work of forming very small cluster embryos but accelerates the rate by lo… ▽ More The heterogeneous condensation of a Lennard-Jones vapor onto a nanoscale seed particle is studied using molecular dynamics simulations. Measuring the nucleation rate and the height of the free energy barrier using the mean first passage time method shows that the presence of a weakly interacting seed has little effect on the work of forming very small cluster embryos but accelerates the rate by lowering the barrier for larger clusters. We suggest that this results from a competition between the energetic and entropic features of cluster formation in the bulk and at the heterogeneity. As the interaction is increased, the free energy of formation is reduced for all cluster sizes. We also develop a simple phenomenological model of film formation on a small seed that captures the general features of the nucleation process for small heterogeneities. A comparison of our simulation results with the model shows that heterogeneous classical nucleation theory provides a good estimate of the critical size of the film but significantly over-estimates the size of the barrier. △ Less

Submitted 4 March, 2011; originally announced March 2011.

Comments: 9 pages, 10 figures, In Print J. Chem. Phys

Edwards entropy and compactivity in a model of granular matter

Authors: Richard K. Bowles, S. S. Ashwin

Abstract: Formulating a statistical mechanics for granular matter remains a significant challenge, in part, due to the difficulty associated with a complete characterization of the systems under study. We present a fully characterized model of a granular material consisting of $N$ two-dimensional, frictionless, hard discs, confined between hard walls, including a complete enumeration of all possible jammed… ▽ More Formulating a statistical mechanics for granular matter remains a significant challenge, in part, due to the difficulty associated with a complete characterization of the systems under study. We present a fully characterized model of a granular material consisting of $N$ two-dimensional, frictionless, hard discs, confined between hard walls, including a complete enumeration of all possible jammed structures. We show the properties of the jammed packings are independent of the distribution of defects within the system and that all the packings are isostatic. This suggests the assumption of equal probability for states of equal volume, which provides one possible way of constructing the equivalent of a microcanonical ensemble, is likely to be vaild for our model. An application of the second law of thermodynamics involving two subsystems in contact shows that the expected spontaneous equilibration of defects between the two is accompanied by an increase in entropy and that the equilibirum, obtained by entropy maximization, is characterized by the equality of compactivities. Finally, we explore the properties of the equivalent to the canonical ensemble for this system. △ Less

Submitted 1 February, 2011; originally announced February 2011.

Comments: Accepted PRE

Crystal Nucleation in a Supercooled Liquid with Glassy Dynamics

Authors: Ivan Saika-Voivod, Richard K. Bowles, Peter H. Poole

Abstract: In simulations of supercooled, high-density liquid silica we study a range of temperature T in which we find both crystal nucleation, as well as the characteristic dynamics of a glass forming liquid, including a breakdown of the Stokes-Einstein relation. We find that the liquid cannot be observed below a homogeneous nucleation limit (HNL) at which the liquid crystallizes faster than it can equil… ▽ More In simulations of supercooled, high-density liquid silica we study a range of temperature T in which we find both crystal nucleation, as well as the characteristic dynamics of a glass forming liquid, including a breakdown of the Stokes-Einstein relation. We find that the liquid cannot be observed below a homogeneous nucleation limit (HNL) at which the liquid crystallizes faster than it can equilibrate. We show that the HNL would occur at lower T, and perhaps not at all, if the Stokes-Einstein relation were obeyed, and hence that glassy dynamics plays a central role in setting a crystallization limit on the liquid state in this case. We also explore the relation of the HNL to the Kauzmann temperature, and test for spinodal-like effects near the HNL. △ Less

Submitted 19 October, 2009; originally announced October 2009.

Comments: 4 pages, 4 figures

Journal ref: Phys. Rev. Lett. 103, 225701 (2009)

The Complete Jamming Landscape of Confined Hard Discs

Authors: S. S. Ashwin, Richard k. Bowles

Abstract: An exact description of the complete jamming landscape is developed for a system of hard discs of diameter $σ$, confined between two lines separated by a distance $1+\sqrt{3/4} < H/σ< 2$. By considering all possible local packing arrangements, the generalized ensemble partition function of jammed states is obtained using the transfer matrix method, which allows us to calculate the configurationa… ▽ More An exact description of the complete jamming landscape is developed for a system of hard discs of diameter $σ$, confined between two lines separated by a distance $1+\sqrt{3/4} < H/σ< 2$. By considering all possible local packing arrangements, the generalized ensemble partition function of jammed states is obtained using the transfer matrix method, which allows us to calculate the configurational entropy and the equation of state for the packings. Exploring the relationship between structural order and packing density, we find that the geometric frustration between local packing environments plays an important role in determining the density distribution of jammed states and that structural "randomness" is a non-monotonic function of packing density. Molecular dynamics simulations show that the properties of the equilibrium liquid are closely related to those of the landscape. △ Less

Submitted 25 March, 2009; originally announced March 2009.

Comments: 5 Pages, 4 figures

A Tiling Approach to Counting Inherent Structures in Hard Potential Systems

Authors: S. S. Ashwin, Richard K Bowles

Abstract: The number of distinguishable inherent structures of a liquid is the key component to understanding the thermodynamics of glass formers. In the case of hard potential systems such as hard discs, spheres and ellipsoids, an inherent structure corresponds to a collectively jammed configuration. This work develops a tiling based approach to counting inherent structures that constructs packings by co… ▽ More The number of distinguishable inherent structures of a liquid is the key component to understanding the thermodynamics of glass formers. In the case of hard potential systems such as hard discs, spheres and ellipsoids, an inherent structure corresponds to a collectively jammed configuration. This work develops a tiling based approach to counting inherent structures that constructs packings by combining sets of elementary locally jammed structures but eliminates those final packings that either, do not tile space, or are not collectively jammed, through the use of tile incompatibility rules. The resulting theory contains a number of geometric quantities, such as the number of growth sites on a tile and the number of tile compatibilities that provide insight into the number of inherent structures in certain limits. We also show that these geometric quantities become quite simple in a system of highly confined hard discs. △ Less

Submitted 10 March, 2009; originally announced March 2009.

Comments: 12 pages, 2 figures. Accepted for Publication : Journal of Non-Crystalline Solids

A Limit of Stability in Supercooled Liquid Clusters

Authors: Eduardo Mendez-Villuendas Ivan Saika-Voivod Richard K. Bowles

Abstract: We examine the metastable liquid phase of a supercooled gold nanocluster by studying the free energy landscape using the largest solid-like embryo as an order parameter. Just below freezing, the free energy exhibits a local minimum at small embryo sizes and a maximum at larger embryo sizes which denotes the critical embryo size. At T=660K the free energy becomes a monotonically decreasing functi… ▽ More We examine the metastable liquid phase of a supercooled gold nanocluster by studying the free energy landscape using the largest solid-like embryo as an order parameter. Just below freezing, the free energy exhibits a local minimum at small embryo sizes and a maximum at larger embryo sizes which denotes the critical embryo size. At T=660K the free energy becomes a monotonically decreasing function of the order parameter as the liquid phase becomes unstable, indicating we have reached a spinodal. In contrast to the usual mean-field theory predictions, the size of the critical embryo remains finite as the spinodal is approached. We also calculate the rate of nucleation, independently from our free energy calculations, and observe a rapid increase in its temperature dependence when the free energy barrier is in the order of $kT$. This supports the idea that freezing becomes a barrierless process around the spinodal temperature. △ Less

Submitted 14 May, 2007; originally announced May 2007.

Comments: 5 pages, 4 figures

Surface Nucleation in the Freezing of Gold Nanoparticles

Authors: Eduardo Mendez-Villuendas Richard K. Bowles

Abstract: We use molecular simulation to calculate the nucleation free energy barrier for the freezing of a 456 atom gold cluster over a range of temperatures. The results show that the embryo of the solid cluster grows at the vapor-surface interface for all temperatures studied and that the usual classical nucleation model, with the embryo growing in the core of the cluster, is unable to predict the shap… ▽ More We use molecular simulation to calculate the nucleation free energy barrier for the freezing of a 456 atom gold cluster over a range of temperatures. The results show that the embryo of the solid cluster grows at the vapor-surface interface for all temperatures studied and that the usual classical nucleation model, with the embryo growing in the core of the cluster, is unable to predict the shape of the free energy barrier. We use a simple partial wetting model that treats the crystal as a lens shaped nucleus at the liquid-vapor interface and find that the line tension plays an important role in the freezing of gold nanoparticles. △ Less

Submitted 3 April, 2007; v1 submitted 26 February, 2007; originally announced February 2007.

Comments: 4 pages, 5 figures. Accepted Phys. Rev. Lett.

The role of FCC tetrahedral subunits in the phase behavior of medium sized Lennard-Jones clusters

Authors: Ivan Saika-Voivod, Louis Poon, Richard K. Bowles

Abstract: The free energy of a 600-atom Lennard-Jones cluster is calculated as a function of surface and bulk crystallinity in order to study the structural transformations that occur in the core of medium sized clusters. Within the order parameter range studied, we find the existence of two free energy minima at temperatures near freezing. One minimum, at low values of both bulk and surface order, belongs… ▽ More The free energy of a 600-atom Lennard-Jones cluster is calculated as a function of surface and bulk crystallinity in order to study the structural transformations that occur in the core of medium sized clusters. Within the order parameter range studied, we find the existence of two free energy minima at temperatures near freezing. One minimum, at low values of both bulk and surface order, belongs to the liquid phase. The second minimum exhibits a highly ordered core with a disordered surface and is related to structures containing a single FCC-tetrahedral subunit, with an edge length of seven atoms (l=7), located in the particle core. At lower temperatures, a third minimum appears at intermediate values of the bulk order parameter which is shown to be related to the formation of multiple l=6 tetrahedra in the core of the cluster. We also use molecular dynamics simulations to follow a series of nucleation events and find that the clusters freeze to structures containing l=5,6,7 and 8 sized tetrahedra as well as those containing no tetrahedral units. The structural correlations between bulk and surface order with the size of the tetrahedral units in the cluster core are examined. Finally, the relationships between the formation of FCC tetrahedral subunits in the core, the phase behavior of medium sized clusters and the nucleation of noncrystalline global structures such as icosahedra and decahedra are discussed. △ Less

Submitted 26 May, 2010; v1 submitted 19 April, 2006; originally announced April 2006.

Comments: 8 pages, 7 figures

Test of classical nucleation theory on deeply supercooled high-pressure simulated silica

Authors: Ivan Saika-Voivod, Peter H. Poole, Richard K. Bowles

Abstract: We test classical nucleation theory (CNT) in the case of simulations of deeply supercooled, high density liquid silica, as modelled by the BKS potential. We find that at density $ρ=4.38$~g/cm$^3$, spontaneous nucleation of crystalline stishovite occurs in conventional molecular dynamics simulations at temperature T=3000 K, and we evaluate the nucleation rate J directly at this T via "brute force… ▽ More We test classical nucleation theory (CNT) in the case of simulations of deeply supercooled, high density liquid silica, as modelled by the BKS potential. We find that at density $ρ=4.38$~g/cm$^3$, spontaneous nucleation of crystalline stishovite occurs in conventional molecular dynamics simulations at temperature T=3000 K, and we evaluate the nucleation rate J directly at this T via "brute force" sampling of nucleation events. We then use parallel, constrained Monte Carlo simulations to evaluate $ΔG(n)$, the free energy to form a crystalline embryo containing n silicon atoms, at T=3000, 3100, 3200 and 3300 K. We find that the prediction of CNT for the n-dependence of $ΔG(n)$ fits reasonably well to the data at all T studied, and at 3300 K yields a chemical potential difference between liquid and stishovite that matches independent calculation. We find that $n^*$, the size of the critical nucleus, is approximately 10 silicon atoms at T=3300 K. At 3000 K, $n^*$ decreases to approximately 3, and at such small sizes methodological challenges arise in the evaluation of $ΔG(n)$ when using standard techniques; indeed even the thermodynamic stability of the supercooled liquid comes into question under these conditions. We therefore present a modified approach that permits an estimation of $ΔG(n)$ at 3000 K. Finally, we directly evaluate at T=3000 K the kinetic prefactors in the CNT expression for J, and find physically reasonable values; e.g. the diffusion length that Si atoms must travel in order to move from the liquid to the crystal embryo is approximately 0.2 nm. We are thereby able to compare the results for J at 3000 K obtained both directly and based on CNT, and find that they agree within an order of magnitude. △ Less

Submitted 18 April, 2006; v1 submitted 25 January, 2006; originally announced January 2006.

Comments: corrected calculation, new figure, accepted in JCP

Landscapes, dynamic heterogeneity and kinetic facilitation in a simple off-lattice model

Authors: R. K. Bowles, I. Saika-Voivod

Abstract: We present a simple off-lattice hard-disc model that exhibits glassy dynamics. The inherent structures are enumerated exactly, transitions between metabasins are well understood, and the particle configurations that act to facilitate dynamics are easily identified. The model readily maps to a coarse grained dynamic facilitation description. We present a simple off-lattice hard-disc model that exhibits glassy dynamics. The inherent structures are enumerated exactly, transitions between metabasins are well understood, and the particle configurations that act to facilitate dynamics are easily identified. The model readily maps to a coarse grained dynamic facilitation description. △ Less

Submitted 15 November, 2005; originally announced November 2005.

Comments: 5 pages, 5 figures, submitted to PRE

Osmotic force resisting chain insertion in a colloidal suspension

Authors: M. Castelnovo, R. K. Bowles, H. Reiss, W. M. Gelbart

Abstract: We consider the problem of inserting a stiff chain into a colloidal suspension of particles that interact with it through excluded volume forces. The free energy of insertion is associated with the work of creating a cavity devoid of colloid and sufficiently large to accomodate the chain. The corresponding work per unit length is the force that resists the entry of the chain into the colloidal s… ▽ More We consider the problem of inserting a stiff chain into a colloidal suspension of particles that interact with it through excluded volume forces. The free energy of insertion is associated with the work of creating a cavity devoid of colloid and sufficiently large to accomodate the chain. The corresponding work per unit length is the force that resists the entry of the chain into the colloidal suspension. In the case of a hard sphere fluid, this work can be calculated straightforwardly within the scaled particle theory; for solutions of flexible polymers, on the other hand, we employ simple scaling arguments. The forces computed in these ways are shown, for nanometer chain and colloid diameters, to be of the order of tens of pN for solution volume fraction for biophysical processes such as the ejection of DNA from viral capsids into the cell cytoplasm. △ Less

Submitted 3 March, 2003; originally announced March 2003.

Comments: 16 pages,3 figures. Accepted for publication in European Physical Journal E