-
Validating Automated Resonance Evaluation with Synthetic Data
Authors:
Oleksii Zivenko,
Noah A. W. Walton,
William Fritsch,
Jacob Forbes,
Amanda M. Lewis,
Aaron Clark,
Jesse M. Brown,
Vladimir Sobes
Abstract:
The integrity and precision of nuclear data are crucial for a broad spectrum of applications, from national security and nuclear reactor design to medical diagnostics, where the associated uncertainties can significantly impact outcomes. A substantial portion of uncertainty in nuclear data originates from the subjective biases in the evaluation process, a crucial phase in the nuclear data producti…
▽ More
The integrity and precision of nuclear data are crucial for a broad spectrum of applications, from national security and nuclear reactor design to medical diagnostics, where the associated uncertainties can significantly impact outcomes. A substantial portion of uncertainty in nuclear data originates from the subjective biases in the evaluation process, a crucial phase in the nuclear data production pipeline. Recent advancements indicate that automation of certain routines can mitigate these biases, thereby standardizing the evaluation process, reducing uncertainty and enhancing reproducibility. This article contributes to developing a framework for automated evaluation techniques testing, emphasizing automated fitting methods that do not require the user to provide any prior information. This approach simplifies the process and reduces the manual effort needed in the initial evaluation stage. It highlights the capability of the framework to validate and optimize subroutines, targeting the performance analysis and optimization of the fitting procedure using high-fidelity synthetic data (labeled experimental data) and the concept of a fully controlled computational experiment. An error metric is introduced to provide a clear and intuitive measure of the fitting quality by quantifying the accuracy and performance across the specified energy. This metric sets a scale for comparison and optimization of routines or hyperparameter selection, improving the entire evaluation process methodology and increasing reproducibility and objectivity.
△ Less
Submitted 18 June, 2024; v1 submitted 3 June, 2024;
originally announced June 2024.
-
Predicting the Electronic Density Response of Condensed-Phase Systems to Electric Field Perturbations
Authors:
Alan M Lewis,
Paolo Lazzaroni,
Mariana Rossi
Abstract:
We present a local and transferable machine learning approach capable of predicting the real-space density response of both molecules and periodic systems to external homogeneous electric fields. The new method, SALTER, builds on the Symmetry-Adapted Gaussian Process Regression SALTED framework. SALTER requires only a small, but necessary, modification to the descriptors used to represent the atom…
▽ More
We present a local and transferable machine learning approach capable of predicting the real-space density response of both molecules and periodic systems to external homogeneous electric fields. The new method, SALTER, builds on the Symmetry-Adapted Gaussian Process Regression SALTED framework. SALTER requires only a small, but necessary, modification to the descriptors used to represent the atomic environments. We present the performance of the method on isolated water molecules, bulk water and a naphthalene crystal. Root mean square errors of the predicted density response lie at or below 10% with barely more than 100 training structures. Derived quantities, such as polarizability tensors and even Raman spectra further derived from these tensors show a good agreement with those calculated directly from quantum mechanical methods. Therefore, SALTER shows excellent performance when predicting derived quantities, while retaining all of the information contained in the full electronic response. This method is thus capable of learning vector fields in a chemical context and serves as a landmark for further developments.
△ Less
Submitted 18 April, 2023;
originally announced April 2023.
-
Electronic-structure properties from atom-centered predictions of the electron density
Authors:
Andrea Grisafi,
Alan M. Lewis,
Mariana Rossi,
Michele Ceriotti
Abstract:
The electron density of a molecule or material has recently received major attention as a target quantity of machine-learning models. A natural choice to construct a model that yields transferable and linear-scaling predictions is to represent the scalar field using a multi-centered atomic basis analogous to that routinely used in density fitting approximations. However, the non-orthogonality of t…
▽ More
The electron density of a molecule or material has recently received major attention as a target quantity of machine-learning models. A natural choice to construct a model that yields transferable and linear-scaling predictions is to represent the scalar field using a multi-centered atomic basis analogous to that routinely used in density fitting approximations. However, the non-orthogonality of the basis poses challenges for the learning exercise, as it requires accounting for all the atomic density components at once. We devise a gradient-based approach to directly minimize the loss function of the regression problem in an optimized and highly sparse feature space. In so doing, we overcome the limitations associated with adopting an atom-centered model to learn the electron density over arbitrarily complex datasets, obtaining extremely accurate predictions. The enhanced framework is tested on 32-molecule periodic cells of liquid water, presenting enough complexity to require an optimal balance between accuracy and computational efficiency. We show that starting from the predicted density a single Kohn-Sham diagonalization step can be performed to access total energy components that carry an error of just 0.1 meV/atom with respect to the reference density functional calculations. Finally, we test our method on the highly heterogeneous QM9 benchmark dataset, showing that a small fraction of the training data is enough to derive ground-state total energies within chemical accuracy.
△ Less
Submitted 28 June, 2022;
originally announced June 2022.
-
Quantum communications infrastructure architecture: theoretical background, network structure and technologies. A review of recent studies from a European public infrastructure perspective
Authors:
Adam M. Lewis,
Petra F. Scudo
Abstract:
Progress in the development of techniques for the construction of multiuser quantum communications networks is reviewed in light of the plans for an EU quantum communications infrastructure (EU QCI). Quantum key distribution networks may be classified as trusted node; quantum repeater-entanglement-swapping based; point-to-multipoint based on passive beam splitters, active optical switchers or WDM;…
▽ More
Progress in the development of techniques for the construction of multiuser quantum communications networks is reviewed in light of the plans for an EU quantum communications infrastructure (EU QCI). Quantum key distribution networks may be classified as trusted node; quantum repeater-entanglement-swapping based; point-to-multipoint based on passive beam splitters, active optical switchers or WDM; high dimensional-multipartite entanglement and flexible reconfigurable multipartite networks. Both satellite and terrestrial implementations are possible and ore both planned for EU QCI; their relative merits are discussed. Current technology falls short in several respects of what is required to address use cases well. Bit rates as a function of distance are currently limited by the characteristics of available devices and are an order of magnitude or more below the theoretical bounds. Non-ideal device behaviour can open loopholes for side-channel attacks. Limited range constrains network geometry. New techniques are being developed to overcome these problems. The more sophisticated schemes depend on ongoing advances in the generation and distribution of entanglement. Particularly promising are the measurement device independent protocol, which eliminates photon-detector related security weaknesses, and the twin-field protocol, which uses similar techniques to extend range, in principle even beyond previously accepted theoretical bounds. Agreement of keys between more than two users can be accomplished classically after pairwise quantum key exchange but direct quantum multipartite agreement using may offer advantages. The same applies to secret sharing. Practical quantum repeaters, to eliminate the need for trusted nodes, are closer to being realised due to recent progress in optical memories.
△ Less
Submitted 5 January, 2022; v1 submitted 13 October, 2021;
originally announced October 2021.
-
Evaluation of $^{86}$Kr Cross Sections For Use in Fusion Diagnostics
Authors:
M. Vorabbi,
G. P. A. Nobre,
D. A. Brown,
A. M. Lewis,
E. Rubino,
S. Mughabghab
Abstract:
The National Ignition Facility at Lawrence Livermore National Laboratory uses $^{86}$Kr as a diagnostic tool to measure the neutron flux produced by fusion reactions. As krypton is chemically inert, it can be implanted directly into the fuel capsule, and the reaction products can be measured to determine the flux of fusion neutrons. $^{86}$Kr cross sections also provide model constraints for the…
▽ More
The National Ignition Facility at Lawrence Livermore National Laboratory uses $^{86}$Kr as a diagnostic tool to measure the neutron flux produced by fusion reactions. As krypton is chemically inert, it can be implanted directly into the fuel capsule, and the reaction products can be measured to determine the flux of fusion neutrons. $^{86}$Kr cross sections also provide model constraints for the $^{85}$Kr branching point in the s-process and the neutron flux in stars. In this work, experimental data on the neutron production, radiative capture, inelastic scattering, and total cross sections of $^{86}$Kr were used in conjunction with the fast region nuclear reaction code EMPIRE and a new resonance-region evaluation to produce a new evaluation of neutron-induced reactions on $^{86}$Kr. For the EMPIRE calculations, we fitted the optical model potential up to 12 MeV to simultaneously reproduce the experimental data for the total cross section and the main inelastic gamma transition from the $2^+$ state to the $0^+$ ground state. For energies above 12 MeV, due to large fluctuations and uncertainties in the total cross section data, we preferred to adopt the Koning-Delaroche global spherical optical model potential. With these models and corrections to the structure of $^{86}$Kr, the evaluated cross sections matched the experimental data. The new evaluation has been submitted for incorporation in the next release of the ENDF/B nuclear reaction library.
△ Less
Submitted 16 September, 2021;
originally announced September 2021.
-
Learning electron densities in the condensed phase
Authors:
Alan M. Lewis,
Andrea Grisafi,
Michele Ceriotti,
Mariana Rossi
Abstract:
We introduce a local machine-learning method for predicting the electron densities of periodic systems. The framework is based on a numerical, atom-centred auxiliary basis, which enables an accurate expansion of the all-electron density in a form suitable for learning isolated and periodic systems alike. We show that using this formulation the electron densities of metals, semiconductors and molec…
▽ More
We introduce a local machine-learning method for predicting the electron densities of periodic systems. The framework is based on a numerical, atom-centred auxiliary basis, which enables an accurate expansion of the all-electron density in a form suitable for learning isolated and periodic systems alike. We show that using this formulation the electron densities of metals, semiconductors and molecular crystals can all be accurately predicted using symmetry-adapted Gaussian process regression models, properly adjusted for the non-orthogonal nature of the basis. These predicted densities enable the efficient calculation of electronic properties which present errors on the order of tens of meV/atom when compared to ab initio density-functional calculations. We demonstrate the key power of this approach by using a model trained on ice unit cells containing only 4 water molecules to predict the electron densities of cells containing up to 512 molecules, and see no increase in the magnitude of the errors of derived electronic properties when increasing the system size. Indeed, we find that these extrapolated derived energies are more accurate than those predicted using a direct machine-learning model. Finally, on heterogeneous datasets SALTED can predict electron densities with errors below 4%.
△ Less
Submitted 9 November, 2021; v1 submitted 9 June, 2021;
originally announced June 2021.
-
Measurement and Modeling of Proton-Induced Reactions on Arsenic from 35 to 200 MeV
Authors:
Morgan B. Fox,
Andrew S. Voyles,
Jonathan T. Morrell,
Lee A. Bernstein,
Jon C. Batchelder,
Eva R. Birnbaum,
Cathy S. Cutler,
Arjan J. Koning,
Amanda M. Lewis,
Dmitri G. Medvedev,
Francois M. Nortier,
Ellen M. O'Brien,
Christiaan Vermeulen
Abstract:
$^{72}$As is a promising positron emitter for diagnostic imaging that can be employed locally using a $^{72}$Se generator. However, current reaction pathways to $^{72}…
▽ More
$^{72}$As is a promising positron emitter for diagnostic imaging that can be employed locally using a $^{72}$Se generator. However, current reaction pathways to $^{72}$Se have insufficient nuclear data for efficient production using regional 100-200 MeV high-intensity proton accelerators. In order to address this deficiency, stacked-target irradiations were performed at LBNL, LANL, and BNL to measure the production of the $^{72}$Se/$^{72}$As PET generator system via $^{75}$As(p,x) between 35 and 200 MeV. This work provides the most well-characterized excitation function for $^{75}$As(p,4n)$^{72}$Se starting from threshold. Additional focus was given to report the first measurements of $^{75}$As(p,x)$^{68}$Ge and bolster an already robust production capability for the highly valuable $^{68}$Ge/$^{68}$Ga PET generator. Thick target yield comparisons with prior established formation routes to both generators are made. In total, high-energy proton-induced cross sections are reported for 55 measured residual products from $^{75}$As, Cu, and Ti targets, where the latter two materials were present as monitor foils. These results were compared with literature data as well as the default theoretical calculations of the nuclear model codes TALYS, CoH, EMPIRE, and ALICE. Reaction modeling at these energies is typically unsatisfactory due to few prior published data and many interacting physics models. Therefore, a detailed assessment of the TALYS code was performed with simultaneous parameter adjustments applied according to a standardized procedure. Particular attention was paid to the formulation of the two-component exciton model in the transition between the compound and pre-equilibrium regions, with a linked investigation of level density models for nuclei off of stability and their impact on modeling predictive power.
△ Less
Submitted 11 June, 2021; v1 submitted 8 June, 2021;
originally announced June 2021.
-
Investigating High-Energy Proton-Induced Reactions on Spherical Nuclei: Implications for the Pre-Equilibrium Exciton Model
Authors:
Morgan B. Fox,
Andrew S. Voyles,
Jonathan T. Morrell,
Lee A. Bernstein,
Amanda M. Lewis,
Arjan J. Koning,
Jon C. Batchelder,
Eva R. Birnbaum,
Cathy S. Cutler,
Dmitri G. Medvedev,
Francois M. Nortier,
Ellen M. O'Brien,
Christiaan Vermeulen
Abstract:
A number of accelerator-based isotope production facilities utilize 100- to 200-MeV proton beams due to the high production rates enabled by high-intensity beam capabilities and the greater diversity of isotope production brought on by the long range of high-energy protons. However, nuclear reaction modeling at these energies can be challenging because of the interplay between different reaction m…
▽ More
A number of accelerator-based isotope production facilities utilize 100- to 200-MeV proton beams due to the high production rates enabled by high-intensity beam capabilities and the greater diversity of isotope production brought on by the long range of high-energy protons. However, nuclear reaction modeling at these energies can be challenging because of the interplay between different reaction modes and a lack of existing guiding cross section data. A Tri-lab collaboration has been formed among the Lawrence Berkeley, Los Alamos, and Brookhaven National Laboratories to address these complexities by characterizing charged-particle nuclear reactions relevant to the production of established and novel radioisotopes. In the inaugural collaboration experiments, stacked-targets of niobium foils were irradiated at the Brookhaven Linac Isotope Producer (E$_p$=200 MeV) and the Los Alamos Isotope Production Facility (E$_p$=100 MeV) to measure $^{93}$Nb(p,x) cross sections between 50 and 200 MeV. The measured cross-section results were compared with literature data as well as the default calculations of the nuclear model codes TALYS, CoH, EMPIRE, and ALICE. We developed a standardized procedure that determines the reaction model parameters that best reproduce the most prominent reaction channels in a physically justifiable manner. The primary focus of the procedure was to determine the best parametrization for the pre-equilibrium two-component exciton model. This modeling study revealed a trend toward a relative decrease for internal transition rates at intermediate proton energies (E$_p$=20-60 MeV) in the current exciton model as compared to the default values. The results of this work are instrumental for the planning, execution, and analysis essential to isotope production.
△ Less
Submitted 5 April, 2021; v1 submitted 24 October, 2020;
originally announced October 2020.
-
The Baghdad Atlas: A relational database of inelastic neutron-scattering $(n,n'γ)$ data
Authors:
A. M. Hurst,
L. A. Bernstein,
T. Kawano,
A. M. Lewis,
K. Song
Abstract:
A relational database has been developed based on the original ($n,n'γ$) work carried out by A. M. Demidov $et$ $al$., at the Nuclear Research Institute in Baghdad, Iraq [$"Atlas$ $of$ $Gamma$-$Ray$ $Spectra$ $from$ $the$ $Inelastic$ $Scattering$ $of$ $Reactor$ $Fast$ $Neutrons"$, Nuclear Research Institute, Baghdad, Iraq (Moscow, Atomizdat 1978)] for 105 independent measurements comprising 76 ele…
▽ More
A relational database has been developed based on the original ($n,n'γ$) work carried out by A. M. Demidov $et$ $al$., at the Nuclear Research Institute in Baghdad, Iraq [$"Atlas$ $of$ $Gamma$-$Ray$ $Spectra$ $from$ $the$ $Inelastic$ $Scattering$ $of$ $Reactor$ $Fast$ $Neutrons"$, Nuclear Research Institute, Baghdad, Iraq (Moscow, Atomizdat 1978)] for 105 independent measurements comprising 76 elemental samples of natural composition and 29 isotopically-enriched samples. The information from this Atlas includes: $γ$-ray energies and relative intensities; nuclide and level data corresponding to the residual nucleus and meta data associated with the target sample that allows for the extraction of the flux-weighted ($n,n'γ$) cross sections for a given transition relative to a defined value. The optimized angular-distribution-corrected fast-neutron flux-weighted partial $γ$-ray cross section for the production of the 846.8-keV $2^{+}_{1} \rightarrow 0^{+}_{\rm gs}$ $γ$-ray transition in $^{56}$Fe, determined to be $\langle σ_γ \rangle = 143(29)$ mb, is used for this purpose. However, different values for the adopted cross section can be readily implemented to accommodate user preference based on revised determinations of this quantity. The Atlas ($n,n'γ$) data has been compiled into a series of CSV-style ASCII data sets and a suite of Python scripts have been developed to build and install the database locally. The database can then be accessed directly through the SQLite engine, or using alternative methods such as the Jupyter Notebook Python-browser interface. Several examples exploiting different interaction methodologies are distributed with the complete software package.
△ Less
Submitted 15 January, 2021; v1 submitted 29 January, 2020;
originally announced January 2020.
-
Proton-induced reactions on Fe, Cu, & Ti from threshold to 55 MeV
Authors:
Andrew S. Voyles,
Amanda M. Lewis,
Jonathan T. Morrell,
M. Shamsuzzoha Basunia,
Lee A. Bernstein,
Jonathan W. Engle,
Stephen A. Graves,
Eric F. Matthews
Abstract:
Theoretical models often differ significantly from measured data in their predictions of the magnitude of nuclear reactions that produce radionuclides for medical, research, and national security applications. In this paper, we compare a priori predictions from several state-of-the-art reaction modeling packages (CoH, EMPIRE, TALYS, and ALICE) to cross sections measured using the stacked-target ac…
▽ More
Theoretical models often differ significantly from measured data in their predictions of the magnitude of nuclear reactions that produce radionuclides for medical, research, and national security applications. In this paper, we compare a priori predictions from several state-of-the-art reaction modeling packages (CoH, EMPIRE, TALYS, and ALICE) to cross sections measured using the stacked-target activation method. The experiment was performed using the LBNL 88-Inch Cyclotron with beams of 25 and 55 MeV protons on a stack of iron, copper, and titanium foils. 34 excitation functions were measured for 4 < Ep < 55 MeV, including the first measurement of the independent cross sections for natFe(p,x) 49,51Cr, 51,52m,52g,56Mn, and 58m,58gCo. All of the models failed to reproduce the isomer-to-ground state ratio for reaction channels at compound and pre-compound energies, suggesting issues in modeling the deposition or distribution of angular momentum in these residual nuclei.
△ Less
Submitted 22 October, 2019;
originally announced October 2019.
-
Ab Initio Linear and Pump-Probe Spectroscopy of Excitons in Molecular Crystals
Authors:
Alan M. Lewis,
Timothy C. Berkelbach
Abstract:
Linear and non-linear spectroscopies are powerful tools used to investigate the energetics and dynamics of electronic excited states of both molecules and crystals. While highly accurate \emph{ab initio} calculations of molecular spectra can be performed relatively routinely, extending these calculations to periodic systems is challenging. Here, we present calculations of the linear absorption spe…
▽ More
Linear and non-linear spectroscopies are powerful tools used to investigate the energetics and dynamics of electronic excited states of both molecules and crystals. While highly accurate \emph{ab initio} calculations of molecular spectra can be performed relatively routinely, extending these calculations to periodic systems is challenging. Here, we present calculations of the linear absorption spectrum and pump-probe two-photon photoemission spectra of the naphthalene crystal using equation-of-motion coupled-cluster theory with single and double excitations (EOM-CCSD). Molecular acene crystals are of interest due to the low-energy multi-exciton singlet states they exhibit, which have been studied extensively as intermediates involved in singlet fission. Our linear absorption spectrum is in good agreement with experiment, predicting a first exciton absorption peak at 4.4 eV, and our two-photon photoemission spectra capture the behavior of multi-exciton states, whose double-excitation character cannot be captured by current methods. The simulated pump-probe spectra provide support for existing interpretations of two-photon photoemission in closely-related acene crystals such as tetracene and pentacene.
△ Less
Submitted 24 September, 2019;
originally announced September 2019.
-
Ab Initio Lifetime and Concomitant Double-Excitation Character of Plasmons at Metallic Densities
Authors:
Alan M. Lewis,
Timothy C. Berkelbach
Abstract:
The accurate calculation of excited state properties of interacting electrons in the condensed phase is an immense challenge in computational physics. Here, we use state-of-the-art equation-of-motion coupled-cluster theory with single and double excitations (EOM-CCSD) to calculate the dynamic structure factor, which can be experimentally measured by inelastic x-ray and electron scattering. Our cal…
▽ More
The accurate calculation of excited state properties of interacting electrons in the condensed phase is an immense challenge in computational physics. Here, we use state-of-the-art equation-of-motion coupled-cluster theory with single and double excitations (EOM-CCSD) to calculate the dynamic structure factor, which can be experimentally measured by inelastic x-ray and electron scattering. Our calculations are performed on the uniform electron gas at densities corresponding to Wigner-Seitz radii of $r_s=5$, 4, and 3 corresponding to the valence electron densities of common metals. We compare our results to those obtained using the random-phase approximation, which is known to provide a reasonable description of the collective plasmon excitation and which resums only a small subset of the polarizability diagrams included in EOM-CCSD. We find that EOM-CCSD, instead of providing a perturbative improvement on the RPA plasmon, predicts a many-state plasmon resonance, where each contributing state has a double-excitation character of 80\% or more. This finding amounts to an ab initio treatment of the plasmon linewidth, which is in good quantitative agreement with previous diagrammatic calculations, and highlights the strongly correlated nature of lifetime effects in condensed-phase electronic structure theory.
△ Less
Submitted 8 July, 2019; v1 submitted 15 October, 2018;
originally announced October 2018.
-
Vertex corrections to the polarizability do not improve the GW approximation for the ionization potential of molecules
Authors:
Alan M. Lewis,
Timothy C. Berkelbach
Abstract:
The $GW$ approximation is based on the neglect of vertex corrections, which appear in the exact self-energy and the exact polarizability. Here, we investigate the importance of vertex corrections in the polarizability only. We calculate the polarizability with equation-of-motion coupled-cluster theory with single and double excitations (EOM-CCSD), which rigorously includes a large class of diagram…
▽ More
The $GW$ approximation is based on the neglect of vertex corrections, which appear in the exact self-energy and the exact polarizability. Here, we investigate the importance of vertex corrections in the polarizability only. We calculate the polarizability with equation-of-motion coupled-cluster theory with single and double excitations (EOM-CCSD), which rigorously includes a large class of diagrammatically-defined vertex corrections beyond the random phase approximation (RPA). As is well-known, the frequency-dependent polarizability predicted by EOM-CCSD is quite different and generally more accurate than that predicted by the RPA. We evaluate the effect of these vertex corrections on a test set of 20 atoms and molecules. When using a Hartree-Fock reference, ionization potentials predicted by the $GW$ approximation with the RPA polarizability are typically overestimated with a mean absolute error of 0.3 eV. However, those predicted with a vertex-corrected polarizability are typically underestimated with an increased mean absolute error of 0.5 eV. This result suggests that vertex corrections in the self-energy cannot be neglected, at least for molecules. We also assess the behavior of eigenvalue self-consistency in vertex-corrected $GW$ calculations, finding a further worsening of the predicted ionization potentials.
△ Less
Submitted 12 April, 2019; v1 submitted 30 September, 2018;
originally announced October 2018.
-
On the low magnetic field effect in radical pair reactions
Authors:
Alan M. Lewis,
Thomas P. Fay,
David E. Manolopoulos,
Christian Kerpal,
Sabine Richert,
Christiane R. Timmel
Abstract:
Radical pair recombination reactions are known to be sensitive to the application of both low and high magnetic fields. The application of a weak magnetic field reduces the singlet yield of a singlet-born radical pair, whereas the application of a strong magnetic field increases the singlet yield. The high field effect arises from energy conservation: when the magnetic field is stronger than the s…
▽ More
Radical pair recombination reactions are known to be sensitive to the application of both low and high magnetic fields. The application of a weak magnetic field reduces the singlet yield of a singlet-born radical pair, whereas the application of a strong magnetic field increases the singlet yield. The high field effect arises from energy conservation: when the magnetic field is stronger than the sum of the hyperfine fields in the two radicals, ${\rm S}\to {\rm T}_{\pm}$ transitions become energetically forbidden, thereby reducing the number of pathways for singlet to triplet interconversion. The low field effect arises from symmetry breaking: the application of a weak magnetic field lifts degeneracies among the zero field eigenstates and increases the number of pathways for singlet to triplet interconversion. However, the details of this effect are more subtle, and have not previously been properly explained. Here we present a complete analysis of the low field effect in a radical pair containing a single proton, and in a radical pair in which one of the radicals contains a large number of hyperfine-coupled nuclear spins. We find that the new transitions that occur when the field is switched on are between ${\rm S}$ and ${\rm T}_0$ in both cases, and not between ${\rm S}$ and ${\rm T}_{\pm}$ as has previously been claimed. We then illustrate this result by using it in conjunction with semiclassical spin dynamics simulations to account for the observation of a biphasic--triphasic--biphasic transition with increasing magnetic field strength in the magnetic field effect on the time-dependent survival probability of a photoexcited carotenoid-porphyrin-fullerene radical pair.
△ Less
Submitted 17 July, 2018;
originally announced July 2018.
-
Excitation functions for (p,x) reactions of niobium in the energy range of E$_{\text{p}}$ = 40-90 MeV
Authors:
Andrew S. Voyles,
Lee A. Bernstein,
Eva R. Birnbaum,
Jonathan W. Engle,
Stephen A. Graves,
Toshihiko Kawano,
Amanda M. Lewis,
Francois M. Nortier
Abstract:
A stack of thin Nb foils was irradiated with the 100 MeV proton beam at Los Alamos National Laboratory's Isotope Production Facility, to investigate the $^{93}$Nb(p,4n)$^{90}$Mo nuclear reaction as a monitor for intermediate energy proton experiments and to benchmark state-of-the-art reaction model codes. A set of 38 measured cross sections for $^{\text{nat}}$Nb(p,x) and $^{\text{nat}}$Cu(p,x) rea…
▽ More
A stack of thin Nb foils was irradiated with the 100 MeV proton beam at Los Alamos National Laboratory's Isotope Production Facility, to investigate the $^{93}$Nb(p,4n)$^{90}$Mo nuclear reaction as a monitor for intermediate energy proton experiments and to benchmark state-of-the-art reaction model codes. A set of 38 measured cross sections for $^{\text{nat}}$Nb(p,x) and $^{\text{nat}}$Cu(p,x) reactions between 40-90 MeV, as well as 5 independent measurements of isomer branching ratios, are reported. These are useful in medical and basic science radionuclide productions at intermediate energies. The $^{\text{nat}}$Cu(p,x)$^{56}$Co, $^{\text{nat}}$Cu(p,x)$^{62}$Zn, and $^{\text{nat}}$Cu(p,x)$^{65}$Zn reactions were used to determine proton fluence, and all activities were quantified using HPGe spectrometry. Variance minimization techniques were employed to reduce systematic uncertainties in proton energy and fluence, improving the reliability of these measurements. The measured cross sections are shown to be in excellent agreement with literature values, and have been measured with improved precision compared with previous measurements. This work also reports the first measurement of the $^{\text{nat}}$Nb(p,x)$^{82\text{m}}$Rb reaction, and of the independent cross sections for $^{\text{nat}}$Cu(p,x)$^{52\text{g}}$Mn and $^{\text{nat}}$Nb(p,x)$^{85\text{g}}$Y in the 40-90 MeV region. The effects of $^{\text{nat}}$Si(p,x)$^{22,24}$Na contamination, arising from silicone adhesive in the Kapton tape used to encapsulate the aluminum monitor foils, is also discussed as a cautionary note to future stacked-target cross section measurements. \emph{A priori} predictions of the reaction modeling codes CoH, EMPIRE, and TALYS are compared with experimentally measured values and used to explore the differences between codes for the $^{\text{nat}}$Nb(p,x) and $^{\text{nat}}$Cu(p,x) reactions.
△ Less
Submitted 21 June, 2018; v1 submitted 18 April, 2018;
originally announced April 2018.
-
Spin-dependent charge recombination along para-phenylene molecular wires
Authors:
Thomas P. Fay,
Alan M. Lewis,
David E. Manolopoulos
Abstract:
We have used an efficient new quantum mechanical method for radical pair recombination reactions to study the spin-dependent charge recombination along PTZ$^{\bullet+}$--Ph$_n$--PDI$^{\bullet-}$ molecular wires. By comparing our results to the experimental data of E. Weiss {\em et al.} [J. Am. Chem. Soc. {\bf 126}, 5577 (2004)], we are able to extract the spin-dependent (singlet and triplet) charg…
▽ More
We have used an efficient new quantum mechanical method for radical pair recombination reactions to study the spin-dependent charge recombination along PTZ$^{\bullet+}$--Ph$_n$--PDI$^{\bullet-}$ molecular wires. By comparing our results to the experimental data of E. Weiss {\em et al.} [J. Am. Chem. Soc. {\bf 126}, 5577 (2004)], we are able to extract the spin-dependent (singlet and triplet) charge recombination rate constants for wires with $n=2-5$. These spin-dependent rate constants have not been extracted previously from the experimental data because they require fitting its magnetic field-dependence to the results of quantum spin dynamics simulations. We find that the triplet recombination rate constant decreases exponentially with the length of the wire, consistent with the superexchange mechanism of charge recombination. However, the singlet recombination rate constant is nearly independent of the length of the wire, suggesting that the singlet pathway is dominated by an incoherent hopping mechanism. A simple qualitative explanation for the different behaviours of the two spin-selective charge recombination pathways is provided in terms of Marcus theory. We also find evidence for a magnetic field-independent background contribution to the triplet yield of the charge recombination reaction, and suggest several possible explanations for it. Since none of these explanations is especially compelling given the available experimental evidence, and since the result appears to apply more generally to other molecular wires, we hope that this aspect of our study will stimulate further experimental work.
△ Less
Submitted 30 October, 2017;
originally announced October 2017.
-
An efficient quantum mechanical method for radical pair recombination reactions
Authors:
Alan M. Lewis,
Thomas P. Fay,
David E. Manolopoulos
Abstract:
The standard quantum mechanical expressions for the singlet and triplet survival probabilities and product yields of a radical pair recombination reaction involve a trace over the states in a combined electronic and nuclear spin Hilbert space. If this trace is evaluated deterministically, by performing a separate time-dependent wavepacket calculation for each initial state in the Hilbert space, th…
▽ More
The standard quantum mechanical expressions for the singlet and triplet survival probabilities and product yields of a radical pair recombination reaction involve a trace over the states in a combined electronic and nuclear spin Hilbert space. If this trace is evaluated deterministically, by performing a separate time-dependent wavepacket calculation for each initial state in the Hilbert space, the computational effort scales as $O(Z^2\log Z)$, where $Z$ is the total number of nuclear spin states. Here we show that the trace can also be evaluated stochastically, by exploiting the properties of spin coherent states. This results in a computational effort of $O(MZ\log Z)$, where $M$ is the number of Monte Carlo samples needed for convergence. Example calculations on a strongly-coupled radical pair with $Z>10^6$ show that the singlet yield can be converged to graphical accuracy using just $M=200$ samples, resulting in a speed up by a factor of $>5000$ over a standard deterministic calculation. We expect that this factor will greatly facilitate future quantum mechanical simulations of a wide variety of radical pairs of interest in chemistry and biology.
△ Less
Submitted 30 December, 2016;
originally announced December 2016.
-
Magnetoelectroluminescence in organic light emitting diodes
Authors:
Joseph E. Lawrence,
Alan M. Lewis,
David E. Manolopoulos,
P. J. Hore
Abstract:
The magnetoelectroluminescence of conjugated organic polymer films is widely accepted to arise from a polaron pair mechanism, but their magnetoconductance is less well understood. Here we derive a new relationship between the experimentally measurable magnetoelectroluminescence and magnetoconductance and the theoretically calculable singlet yield of the polaron pair recombination reaction. This re…
▽ More
The magnetoelectroluminescence of conjugated organic polymer films is widely accepted to arise from a polaron pair mechanism, but their magnetoconductance is less well understood. Here we derive a new relationship between the experimentally measurable magnetoelectroluminescence and magnetoconductance and the theoretically calculable singlet yield of the polaron pair recombination reaction. This relationship is expected to be valid regardless of the mechanism of the magnetoconductance, provided the mobilities of the free polarons are independent of the applied magnetic field (i.e., provided one discounts the possibility of spin-dependent transport). We also discuss the semiclassical calculation of the singlet yield of the polaron pair recombination reaction for materials such as poly(2,5-dioctyloxy-paraphenylene vinylene) (DOO-PPV), the hyperfine fields in the polarons of which can be extracted from light-induced electron spin resonance measurements. The resulting theory is shown to give good agreement with experimental data for both normal (H-) and deuterated (D-) DOO-PPV over a wide range of magnetic field strengths once singlet-triplet dephasing is taken into account. Without this effect, which has not been included in any previous simulation of magnetoelectroluminescence, it is not possible to reproduce the experimental data for both isotopologues in a consistent fashion. Our results also indicate that the magnetoconductance of DOO-PPV cannot be solely due to the effect of the magnetic field on the dissociation of polaron pairs.
△ Less
Submitted 19 May, 2016;
originally announced May 2016.
-
Asymmetric recombination and electron spin relaxation in the semiclassical theory of radical pair reactions
Authors:
Alan M. Lewis,
David E. Manolopoulos,
P. J. Hore
Abstract:
We describe how the semiclassical theory of radical pair recombination reactions recently introduced by two of us [D. E. Manolopoulos and P. J. Hore, J. Chem. Phys. 139, 124106 (2013)] can be generalised to allow for different singlet and triplet recombination rates. This is a non-trivial generalisation because when the recombination rates are different the recombination process is dynamically cou…
▽ More
We describe how the semiclassical theory of radical pair recombination reactions recently introduced by two of us [D. E. Manolopoulos and P. J. Hore, J. Chem. Phys. 139, 124106 (2013)] can be generalised to allow for different singlet and triplet recombination rates. This is a non-trivial generalisation because when the recombination rates are different the recombination process is dynamically coupled to the coherent electron spin dynamics of the radical pair. Furthermore, because the recombination operator is a two-electron operator, it is no longer sufficient simply to consider the two electrons as classical vectors: one has to consider the complete set of 16 two-electron spin operators as independent classical variables. The resulting semiclassical theory is first validated by comparison with exact quantum mechanical results for a model radical pair containing 12 nuclear spins. It is then used to shed light on the spin dynamics of a carotenoid-porphyrin-fullerene (CPF) triad containing considerably more nuclear spins which has recently been used to establish a 'proof of principle' for the operation of a chemical compass [K. Maeda et al., Nature 453, 387 (2008)]. We find in particular that the intriguing biphasic behaviour that has been observed in the effect of an Earth-strength magnetic field on the time-dependent survival probability of the photo-excited C+PF- radical pair arises from a delicate balance between its asymmetric recombination and the relaxation of the electron spin in the carotenoid radical.
△ Less
Submitted 8 July, 2014;
originally announced July 2014.
-
Implications of the Cosmic Background Imager Polarization Data
Authors:
J. L. Sievers,
C. Achermann,
J. R. Bond,
L. Bronfman,
R. Bustos,
C. R. Contaldi,
C. Dickinson,
P. G. Ferreira,
M. E. Jones,
A. M. Lewis,
B. S. Mason,
J. May,
S. T. Myers,
S. Padin,
T. J. Pearson,
M. Pospieszalski,
A. C. S. Readhead,
R. Reeves,
A. C. Taylor,
S. Torres
Abstract:
We present new measurements of the power spectra of the E-mode of CMB polarization, the temperature T, the cross-correlation of E and T, and upper limits on the B-mode from 2.5 years of dedicated Cosmic Background Imager (CBI) observations. Both raw maps and optimal signal images in the uv-plane and real space show strong detections of the E-mode (11.7 sigma for the EE power spectrum overall) an…
▽ More
We present new measurements of the power spectra of the E-mode of CMB polarization, the temperature T, the cross-correlation of E and T, and upper limits on the B-mode from 2.5 years of dedicated Cosmic Background Imager (CBI) observations. Both raw maps and optimal signal images in the uv-plane and real space show strong detections of the E-mode (11.7 sigma for the EE power spectrum overall) and no detection of the B-mode. The power spectra are used to constrain parameters of the flat tilted adiabatic Lambda-CDM models: those determined from EE and TE bandpowers agree with those from TT, a powerful consistency check. There is little tolerance for shifting polarization peaks from the TT-forecast locations, as measured by the angular sound crossing scale theta = 100 ell_s = 1.03 +/- 0.02 from EE and TE cf. 1.044 +/- 0.005 with the TT data included. The scope for extra out-of-phase peaks from subdominant isocurvature modes is also curtailed. The EE and TE measurements of CBI, DASI and BOOMERANG are mutually consistent, and, taken together rather than singly, give enhanced leverage for these tests.
△ Less
Submitted 9 November, 2006; v1 submitted 8 September, 2005;
originally announced September 2005.
-
Cosmological Parameters from the 2003 flight of BOOMERANG
Authors:
C. J. MacTavish,
P. A. R. Ade,
J. J. Bock,
J. R. Bond,
J. Borrill,
A. Boscaleri,
P. Cabella,
C. R. Contaldi,
B. P. Crill,
P. de Bernardis,
G. De Gasperis,
A. de Oliveira-Costa,
G. De Troia,
G. Di Stefano,
E. Hivon,
A. H. Jaffe,
W. C. Jones,
T. S. Kisner,
A. E. Lange,
A. M. Lewis,
S. Masi,
P. D. Mauskopf,
A. Melchiorri,
T. E. Montroy,
P. Natoli
, et al. (13 additional authors not shown)
Abstract:
We present the cosmological parameters from the CMB intensity and polarization power spectra of the 2003 Antarctic flight of the BOOMERANG telescope. The BOOMERANG data alone constrains the parameters of the $Λ$CDM model remarkably well and is consistent with constraints from a multi-experiment combined CMB data set. We add LSS data from the 2dF and SDSS redshift surveys to the combined CMB data…
▽ More
We present the cosmological parameters from the CMB intensity and polarization power spectra of the 2003 Antarctic flight of the BOOMERANG telescope. The BOOMERANG data alone constrains the parameters of the $Λ$CDM model remarkably well and is consistent with constraints from a multi-experiment combined CMB data set. We add LSS data from the 2dF and SDSS redshift surveys to the combined CMB data set and test several extensions to the standard model including: running of the spectral index, curvature, tensor modes, the effect of massive neutrinos, and an effective equation of state for dark energy. We also include an analysis of constraints to a model which allows a CDM isocurvature admixture.
△ Less
Submitted 21 July, 2005;
originally announced July 2005.
-
The Cosmic Microwave Background and Inflation Parameters
Authors:
J. R. Bond,
C. R. Contaldi,
A. M. Lewis,
D. Pogosyan
Abstract:
We review the currrent cosmic parameter determinations of relevance to inflation using the WMAP-1year, Boomerang, CBI, Acbar and other CMB data. The basic steps in the pipelines which determine the bandpowers from the raw data from which these estimations are made are summarized. We forecast how the precision is likely to improve with more years of WMAP in combination with future ground-based ex…
▽ More
We review the currrent cosmic parameter determinations of relevance to inflation using the WMAP-1year, Boomerang, CBI, Acbar and other CMB data. The basic steps in the pipelines which determine the bandpowers from the raw data from which these estimations are made are summarized. We forecast how the precision is likely to improve with more years of WMAP in combination with future ground-based experiments and with Planck. We address whether the current data indicates strong breaking from uniform acceleration through the relatively small region of the inflaton potential that the CMB probes, manifest in the much-discussed running spectral index or in even more radical braking/breaking scenarios. Although some weak ``anomalies'' appear in the current data, the statistical case is not there. However increased precision, at the high multipole end and with polarization measurements, will significantly curtail current freedom.
△ Less
Submitted 7 June, 2004;
originally announced June 2004.
-
Large Scale Cosmic Microwave Background Anisotropies and Dark Energy
Authors:
J. Weller,
A. M. Lewis
Abstract:
In this note we investigate the effects of perturbations in a dark energy component with a constant equation of state on large scale cosmic microwave background anisotropies. The inclusion of perturbations increases the large scale power. We investigate more speculative dark energy models with w<-1 and find the opposite behaviour. Overall the inclusion of perturbations in the dark energy compone…
▽ More
In this note we investigate the effects of perturbations in a dark energy component with a constant equation of state on large scale cosmic microwave background anisotropies. The inclusion of perturbations increases the large scale power. We investigate more speculative dark energy models with w<-1 and find the opposite behaviour. Overall the inclusion of perturbations in the dark energy component increases the degeneracies. We generalise the parameterization of the dark energy fluctuations to allow for an arbitrary const ant sound speeds and show how constraints from cosmic microwave background experiments change if this is included. Combining cosmic microwave background with large scale structure, Hubble parameter and Supernovae observations we obtain w=-1.02+-0.16 (1 sigma) as a constraint on the equation of state, which is almost independent of the sound speed chosen. With the presented analysis we find no significant constraint on the constant speed of sound of the dark energy component.
△ Less
Submitted 29 August, 2003; v1 submitted 4 July, 2003;
originally announced July 2003.
-
Reconstructing the primordial power spectrum
Authors:
S. L. Bridle,
A. M. Lewis,
J. Weller,
G. Efstathiou
Abstract:
We reconstruct the shape of the primordial power spectrum from the latest cosmic microwave background data, including the new results from the Wilkinson Microwave Anisotropy Probe (WMAP), and large scale structure data from the two degree field galaxy redshift survey (2dFGRS). We tested four parameterizations taking into account the uncertainties in four cosmological parameters. First we paramet…
▽ More
We reconstruct the shape of the primordial power spectrum from the latest cosmic microwave background data, including the new results from the Wilkinson Microwave Anisotropy Probe (WMAP), and large scale structure data from the two degree field galaxy redshift survey (2dFGRS). We tested four parameterizations taking into account the uncertainties in four cosmological parameters. First we parameterize the initial spectrum by a tilt and a running spectral index, finding marginal evidence for a running spectral index only if the first three WMAP multipoles (ell=2,3,4) are included in the analysis. Secondly, to investigate further the low CMB large scale power, we modify the conventional power-law spectrum by introducing a scale above which there is no power. We find a preferred position of the cut at k_c ~ 3e-4 /Mpc although k_c=0 (no cut) is not ruled out. Thirdly we use a model independent parameterization, with 16 bands in wavenumber, and find no obvious sign of deviation from a power law spectrum on the scales investigated. Furthermore the values of the other cosmological parameters defining the model remain relatively well constrained despite the freedom in the shape of the initial power spectrum. Finally we investigate a model motivated by double inflation, in which the power spectrum has a break between two characteristic wavenumbers. We find that if a break is required to be in the range 0.01 < k Mpc <0.1 then the ratio of amplitudes across the break is constrained to be 1.23 +/- 0.14. Our results are consistent with a power law spectrum that is featureless and close to scale invariant over the wavenumber range 0.005 ~< k Mpc ~< 0.15, with a hint of a decrease in power on the largest scales.
△ Less
Submitted 16 May, 2003; v1 submitted 17 February, 2003;
originally announced February 2003.
-
Deep inelastic events containing two forward jets at HERA
Authors:
J. Kwiecinski,
C. A. M. Lewis,
A. D. Martin
Abstract:
We use the BFKL equation to calculate the rate of deep inelastic scattering events containing two forward jets (adjacent to the proton remnants) at HERA. We compare the production of two forward jets with that of only one forward jet (the "Mueller" process). We obtain a stable prediction for this two to one jet ratio, which may serve as a measure of the BFKL vertex function.
We use the BFKL equation to calculate the rate of deep inelastic scattering events containing two forward jets (adjacent to the proton remnants) at HERA. We compare the production of two forward jets with that of only one forward jet (the "Mueller" process). We obtain a stable prediction for this two to one jet ratio, which may serve as a measure of the BFKL vertex function.
△ Less
Submitted 17 July, 1997;
originally announced July 1997.
-
Observable jets from the BFKL chain
Authors:
J. Kwiecinski,
C. A. M. Lewis,
A. D. Martin
Abstract:
We derive a modified form of the BFKL equation which enables the structure of the gluon emissions to be studied in small $x$ deep inelastic scattering. The equation incorporates the resummation of the virtual and unresolved real gluon emissions. We solve the equation to calculate the number of small $x$ deep-inelastic events containing 0,1,2 ...resolved gluon jets, that is jets with transverse m…
▽ More
We derive a modified form of the BFKL equation which enables the structure of the gluon emissions to be studied in small $x$ deep inelastic scattering. The equation incorporates the resummation of the virtual and unresolved real gluon emissions. We solve the equation to calculate the number of small $x$ deep-inelastic events containing 0,1,2 ...resolved gluon jets, that is jets with transverse momenta $q_{T} > μ$. We study the jet decomposition for different choices of the jet resolution parameter $μ$.
△ Less
Submitted 19 June, 1996;
originally announced June 1996.