Source author record

J. Engel

J. Engel appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

32works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

32 published item(s)

preprint2022arXiv

Ab initio studies of double Gamow-Teller transition and its correlation with neutrinoless double beta decay

We use chiral interactions and several {\em ab initio} methods to compute the nuclear matrix elements (NMEs) for ground-state to ground-state double Gamow-Teller transitions in a range of isotopes, and explore the correlation of these NMEs with those for neutrinoless double beta decay produced by the exchange of a light Majorana neutrino. When all the NMEs of both isospin-conserving and isospin-changing transitions from the {\em ab initio} calculations are considered, the correlation is strong. For the experimentally relevant isospin-changing transitions by themselves, however, the correlation is weaker and may not be helpful for reducing the uncertainty in the NMEs for neutrinoless double-beta decay.

preprint2022arXiv

Solving Nuclear Structure Problems with the Adaptive Variational Quantum Algorithm

We use the Lipkin-Meshkov-Glick (LMG) model and the valence-space nuclear shell model to examine the likely performance of variational quantum eigensolvers in nuclear-structure theory. The LMG model exhibits both a phase transition and spontaneous symmetry breaking at the mean-field level in one of the phases, features that characterize collective dynamics in medium-mass and heavy nuclei. We show that with appropriate modifications, the ADAPT-VQE algorithm, a particularly flexible and accurate variational approach, is not troubled by these complications. We treat up to 12 particles and show that the number of quantum operations needed to approach the ground-state energy scales linearly with the number of qubits. We find similar scaling when the algorithm is applied to the nuclear shell model with realistic interactions in the $sd$ and $pf$ shells. Although most of these simulations contain no noise, we use a noise model from real IBM hardware to show that for the LMG model with four particles, weak noise has no effect on the efficiency of the algorithm.

preprint2021arXiv

Finite-temperature electron-capture rates for neutron-rich nuclei around N=50 and effects on core-collapse supernovae simulations

The temperature dependence of stellar electron-capture (EC) rates is investigated, with a focus on nuclei around $N=50$, just above $Z=28$, which play an important role during the collapse phase of core-collapse supernovae (CCSN). Two new microscopic calculations of stellar EC rates are obtained from a relativistic and a non-relativistic finite-temperature quasiparticle random-phase approximation approaches, for a conventional grid of temperatures and densities. In both approaches, EC rates due to Gamow-Teller transitions are included. In the relativistic calculation contributions from first-forbidden transitions are also included, and add strongly to the EC rates. The new EC rates are compared with large-scale shell model calculations for the specific case of $^{86}$Kr, providing insight into the finite-temperature effects on the EC rates. At relevant thermodynamic conditions for core-collapse, the discrepancies between the different calculations of this work are within about one order of magnitude. Numerical simulations of CCSN are performed with the spherically-symmetric GR1D simulation code to quantify the impact of such differences on the dynamics of the collapse. These simulations also include EC rates based on two parametrized approximations. A comparison of the neutrino luminosities and enclosed mass at core bounce shows that differences between simulations with different sets of EC rates are relatively small ($\approx 5\%$), suggesting that the EC rates used as inputs for these simulations have become well constrained.

preprint2021arXiv

Two-body weak currents in heavy nuclei

In light and medium-mass nuclei, two-body weak currents from chiral effective field theory account for a significant portion of the phenomenological quenching of Gamow-Teller transition matrix elements. Here we examine the systematic effects of two-body axial currents on Gamow-Teller strength and $β$-decay rates in heavy nuclei within energy-density functional theory. Using a Skyrme functional and the charge-changing finite amplitude method, we add the contributions of two-body currents to the usual one-body linear response in the Gamow-Teller channel, both exactly and though a density-matrix expansion. The two-body currents, as expected, usually quench both summed Gamow-Teller strength and decay rates, but by an amount that decreases as the neutron excess grows. In addition, they can enhance individual low-lying transitions, leading to decay rates that are quite different from those that an energy-independent quenching would produce, particularly in neutron-rich nuclei. We show that both these unexpected effects are related to changes in the total nucleon density as the number of neutrons increases.

preprint2020arXiv

Ab Initio Treatment of Collective Correlations and the Neutrinoless Double Beta Decay of $^{48}$Ca

Working with Hamiltonians from chiral effective field theory, we develop a novel framework for describing arbitrary deformed medium-mass nuclei by combining the in-medium similarity renormalization group with the generator coordinate method. The approach leverages the ability of the first method to capture dynamic correlations and the second to include collective correlations without violating symmetries. We use our scheme to compute the matrix element that governs the neutrinoless double beta decay of $^{48}$Ca to $^{48}$Ti, and find it to have the value $0.61$, near or below the predictions of most phenomenological methods. The result opens the door to ab initio calculations of the matrix elements for the decay of heavier nuclei such as $^{76}$Ge, $^{130}$Te, and $^{136}$Xe.

preprint2020arXiv

Benchmark neutrinoless double-beta decay matrix elements in a light nucleus

We compute nuclear matrix elements of neutrinoless double-beta decay mediated by light Majorana-neutrino exchange in the A = 6 system. The goal is to benchmark two many-body approaches, the No-Core Shell Model and the Multi-Reference In-Medium Similarity Renormalization Group. We use the SRG-evolved chiral N3LO-EM500 potential for the nuclear interaction, and make the approximation that isospin is conserved. We compare the results of the two approaches as a function of the cutoff on the many-body basis space. Although differences are seen in the predicted nuclear radii, the ground-state energies and neutrinoless double-beta decay matrix elements produced by the two approaches show significant agreement. We discuss the implications for calculations in heavier nuclei.

preprint2016arXiv

Beta decay of deformed r-process nuclei near A = 80 and A= 160, including odd-A and odd-odd nuclei, with the Skyrme finite-amplitude method

After identifying the nuclei in the regions near A =80 and A = 160 for which beta-decay rates have the greatest effect on weak and main r-process abundance patterns, we apply the finite-amplitude method (FAM) with Skyrme energy-density functionals (EDFs) to calculate beta-decay half-lives of those nuclei in the quasiparticle random-phase approximation (QRPA). We use the equal filling approximation to extend our implementation of the charge-changing FAM, which incorporates pairing correlations and allows axially symmetric deformation, to odd-A and odd-odd nuclei. Within this framework we find differences of up to a factor of seven between our calculated beta-decay half-lives and those of previous efforts. Repeated calculations with nuclei near A = 160 and multiple EDFs show a spread of two to four in beta-decay half-lives, with differences in calculated Q values playing an important role. We investigate the implications of these results for r-process simulations.

preprint2016arXiv

Octupole correlations in low-lying states of 150Nd and 150Sm and their impact on neutrinoless double-beta decay

We present a generator-coordinate calculation, based on a relativistic energy-density functional, of the low-lying spectra in the isotopes $^{150}$Nd and $^{150}$Sm and of the nuclear matrix element that governs the neutrinoless double-beta decay of the first isotope to the second. We carefully examine the impact of octupole correlations on both nuclear structure and the double-beta decay matrix element. Octupole correlations turn out to reduce quadrupole collectivity in both nuclei. Shape fluctuations, however, dilute the effects of octupole deformation on the double-beta decay matrix element, so that the overall octupole-induced quenching is only about 7\%.

preprint2016arXiv

Testing the importance of collective correlations in neutrinoless $ββ$ decay

We investigate the extent to which theories of collective motion can capture the physics that determines the nuclear matrix elements governing neutrinoless double-beta decay. To that end we calculate the matrix elements for a series of isotopes in the full $pf$ shell, omitting no spin-orbit partners. With the inclusion of isoscalar pairing, a separable collective Hamiltonian that is derived from the shell model effective interaction reproduces the full shell-model matrix elements with good accuracy. A version of the generator coordinate method that includes the isoscalar pairing amplitude as a coordinate also reproduces the shell model results well, an encouraging result for theories of collective motion, which can include more single-particle orbitals than the shell model. We briefly examine heavier nuclei relevant for experimental double-beta decay searches, in which shell-model calculations with all spin-orbit partners are not feasible; our estimates suggest that isoscalar pairing also plays a significant role in these nuclei, though one we are less able to quantify precisely.

preprint2015arXiv

Global description of beta-minus decay in even-even nuclei with the axially-deformed Skyrme finite amplitude method

We use the finite amplitude method for computing charge-changing Skyrme-QRPA transition strengths in axially-deformed nuclei together with a modern Skyrme energy-density functional to fit several previously unconstrained parameters in the charge-changing time-odd part of the functional. With the modified functional we then calculate rates of beta-minus decay for all medium-mass and heavy even-even nuclei between the valley of stability and the neutron drip line. We fit the Skyrme parameters to a limited set of beta-decay rates, a set of Gamow-Teller resonance energies, and a set of spin-dipole resonance energies, in both spherical and deformed nuclei. Comparison to available experimental beta-decay rates shows agreement at roughly the same level as in other global QRPA calculations. We estimate the uncertainty in our rates all the way to the neutron drip line through a construction that extrapolates the errors of known beta-decay rates in nuclei with intermediate Q values to less stable isotopes with higher Q values.

preprint2015arXiv

Subtraction method in the second random--phase approximation: first applications with a Skyrme energy functional

We make use of a subtraction procedure, introduced to overcome double--counting problems in beyond--mean--field theories, in the second random--phase--approximation (SRPA) for the first time. This procedure guarantees the stability of SRPA (so that all excitation energies are real). We show that the method fits perfectly into nuclear density--functional theory. We illustrate applications to the monopole and quadrupole response and to low--lying $0^+$ and $2^+$ states in the nucleus $^{16}$O. We show that the subtraction procedure leads to: (i) results that are weakly cutoff dependent; (ii) a considerable reduction of the SRPA downwards shift with respect to the random--phase approximation (RPA) spectra (systematically found in all previous applications). This implementation of the SRPA model will allow a reliable analysis of the effects of 2 particle--2 hole configurations ($2p2h$) on the excitation spectra of medium--mass and heavy nuclei.

preprint2014arXiv

Ab-initio coupled-cluster effective interactions for the shell model: Application to neutron-rich oxygen and carbon isotopes

We derive and compute effective valence-space shell-model interactions from ab-initio coupled-cluster theory and apply them to open-shell and neutron-rich oxygen and carbon isotopes. Our shell-model interactions are based on nucleon-nucleon and three-nucleon forces from chiral effective-field theory. We compute the energies of ground and low-lying states, and find good agreement with experiment. In particular our calculations are consistent with the N=14, 16 shell closures in oxygen-22 and oxygen-24, while for carbon-20 the corresponding N=14 closure is weaker. We find good agreement between our coupled-cluster effective-interaction results with those obtained from standard single-reference coupled-cluster calculations for up to eight valence neutrons.

preprint2014arXiv

Chiral Two-Body Currents and Neutrinoless Double-Beta Decay in the QRPA

We test the effects of an approximate treatment of two-body contributions to the axial-vector current on the QRPA matrix elements for neutrinoless double-beta decay in a range of isotopes. The form and strength of the two-body terms come from chiral effective-field theory. The two-body currents typically reduce the matrix elements by about 20%, not as much as in shell-model calculations. One reason for the difference is that standard practice in the QRPA is to adjust the strength of the isoscalar pairing interaction to reproduce two-neutrino double-beta decay lifetimes. Another may be the larger QRPA single-particle space. Whatever the reasons, the effects on neutrinoless decay are significantly less than those on two-neutrino decay, both in the shell model and the QRPA.

preprint2014arXiv

Finite Amplitude Method for Charge-Changing Transitions in Axially-Deformed Nuclei

We describe and apply a version of the finite amplitude method for obtaining the charge-changing nuclear response in the quasiparticle random phase approximation. The method is suitable for calculating strength functions and beta-decay rates, both allowed and forbidden, in axially-deformed open-shell nuclei. We demonstrate the speed and versatility of the code through a preliminary examination of the effects of tensor terms in Skyrme functionals on beta decay in a set of spherical and deformed open-shell nuclei. Like the isoscalar pairing interaction, the tensor terms systematically increase allowed beta-decay rates. This finding generalizes previous work in semimagic nuclei and points to the need for a comprehensive study of time-odd terms in nuclear density functionals.

preprint2013arXiv

Large-Scale Calculations of the Double-Beta Decay of 76Ge, 130Te, 136Xe, and 150Nd in the Deformed Self-Consistent Skyrme Quasiparticle Random-Phase Approximation

We use the axially-deformed Skyrme Quasiparticle Random-Phase Approximation (QRPA) together with the SkM* energy-density functional, both as originally presented and with the time-odd part adjusted to reproduce the Gamow-Teller resonance energy in 208Pb, to calculate the matrix elements governing the neutrinoless double-beta decay of 76Ge, 130Te, 136Xe, and 150Nd. Our matrix elements in 130Te and 136Xe are significantly smaller than those of previous QRPA calculations, primarily because of the difference in pairing or deformation between the initial and final nuclei. In 76Ge and 150Nd our results are similar to those of less computationally intensive QRPA calculations. We suspect the 76Ge result, however, because we are forced to use a spherical ground-state, even though the HFB indicates a deformed minimum.

preprint2013arXiv

New determination of double-beta-decay properties in 48Ca: high-precision Q-value measurement and improved nuclear matrix element calculations

We report a direct measurement of the Q-value of the neutrinoless double-beta-decay candidate 48Ca at the TITAN Penning-trap mass spectrometer, with the result that Q = 4267.98(32) keV. We measured the masses of both the mother and daughter nuclides, and in the latter case found a 1 keV deviation from the literature value. In addition to the Q-value, we also present results of a new calculation of the neutrinoless double-beta-decay nuclear matrix element of 48Ca. Using diagrammatic many-body perturbation theory to second order to account for physics outside the valence space, we constructed an effective shell-model double-beta-decay operator, which increased the nuclear matrix element by about 75% compared with that produced by the bare operator. The new Q-value and matrix element strengthen the case for a 48Ca double-beta-decay experiment.

preprint2011arXiv

Jastrow functions in double-beta decay

We use simple analytic considerations and a Monte Carlo calculation of nucleons in a box to argue that the use of Jastrow functions as short-range correlators in the commonly employed two-body-cluster approximation causes significant errors in the matrix elements for double-beta decay. The Jastrow approach appears to agree with others, however, if many-body clusters are included. A careful treatment of the charge-changing analog of the nuclear pair density shows, in addition, that differences between Unitary Correlator Operator Method and Brueckner methods for treating short-range correlations in double-beta are less significant than suggested by previous work.

preprint2011arXiv

Testing Skyrme energy-density functionals with the QRPA in low-lying vibrational states of rare-earth nuclei

Although nuclear energy density functionals are determined primarily by fitting to ground state properties, they are often applied in nuclear astrophysics to excited states, usually through the quasiparticle random phase approximation (QRPA). Here we test the Skyrme functionals SkM* and SLy4 along with the self-consistent QRPA by calculating properties of low-lying vibrational states in a large number of well-deformed even-even rare-earth nuclei. We reproduce trends in energies and transition probabilities associated with gamma-vibrational states, but our results are not perfect and indicate the presences of multi-particle-hole correlations that are not included in the QRPA. The Skyrme functional SkM* performs noticeably better than SLy4. In a few nuclei, changes in the treatment of the pairing energy functional have a significant effect. The QRPA is less successful with "beta-vibrational" states than with the gamma-vibrational states.

preprint2011arXiv

Towards a Microscopic Reaction Description Based on Energy-Density-Functional Structure Models

A microscopic calculation of reaction cross sections for nucleon-nucleus scattering has been performed by explicitly coupling the elastic channel to all particle-hole excitations in the target and one-nucleon pickup channels. The particle-hole states may be regarded as doorway states through which the flux flows to more complicated configurations, and subsequently to long-lived compound nucleus resonances. Target excitations for $^{40,48}$Ca, $^{58}$Ni, $^{90}$Zr and $^{144}$Sm were described in a random-phase framework using a Skyrme functional. Reaction cross sections obtained agree very well with experimental data and predictions of a state-of-the-art fitted optical potential. Couplings between inelastic states were found to be negligible, while the pickup channels contribute significantly. The effect of resonances from higher-order channels was assessed. Elastic angular distributions were also calculated within the same method, achieving good agreement with experimental data. For the first time observed absorptions are completely accounted for by explicit channel coupling, for incident energies between 10 and 70 MeV, with consistent angular distribution results.

preprint2010arXiv

Coupled-channels calculations of nonelastic cross sections using a density-functional structure model

A microscopic calculation of the reaction cross-section for nucleon-nucleus scattering has been performed by explicitly coupling the elastic channel to all particle-hole (p-h) excitation states in the target and to all one-nucleon pickup channels. The p-h states may be regarded as doorway states through which the flux flows to more complicated configurations, and subsequently to long-lived compound nucleus resonances. Target excitations for 40,48Ca, 58Ni, 90Zr and 144Sm were described in a QRPA framework using a Skyrme functional. Reaction cross sections calculated in this approach were compared to predictions of a fitted optical potential and to experimental data, reaching very good agreement. Couplings between inelastic states were found to be negligible, while the couplings to pickup channels contribute significantly. For the first time observed reaction cross-sections are completely accounted for by explicit channel coupling, for incident energies between 10 and 40 MeV.

preprint2010arXiv

Self-consistent Skyrme QRPA for use in axially-symmetric nuclei of arbitrary mass

We describe a new implementation of the quasiparticle random phase approximation (QRPA) in axially-symmetric deformed nuclei with Skyrme and volume-pairing energy-density functionals. After using a variety of tests to demonstrate the accuracy of the code in ^{24,26}Mg and ^{16}O, we report the first fully self-consistent application of the Skyrme QRPA to a heavy deformed nucleus, calculating strength distributions for several K^pi in ^{172}Yb. We present energy-weighted sums, properties of gamma-vibrational and low-energy K^pi=0^+ states, and the complete isovector E1 strength function. The QRPA calculation reproduces the properties of the low-lying 2^+ states as well or better than it typically does in spherical nuclei.

preprint2006arXiv

Multi-layer scintillation detector for the MOON double beta decay experiment: Scintillation photon responses studied by a prototype detector MOON-1

An ensemble of multi-layer scintillators is discussed as an option of the high-sensitivity detector Mo Observatory Of Neutrinos (MOON) for spectroscopic measurements of neutrino-less double beta decays. A prototype detector MOON-1, which consists of 6 layer plastic-scintillator plates, was built to study the sensitivity of the MOON-type detector. The scintillation photon collection and the energy resolution, which are key elements for the high-sensitivity experiments, are found to be 1835+/-30 photo-electrons for 976 keV electrons and sigma = 2.9+/-0.1% (dE/E = 6.8+/-0.3 % in FWHM) at the Qbb ~ 3 MeV region, respectively. The multi-layer plastic-scintillator structure with good energy resolution as well as good background suppression of beta-gamma rays is crucial for the MOON-type detector to achieve the inverted hierarchy neutrino mass sensitivity.

preprint2005arXiv

Time-Reversal-Violating Schiff Moment of 199Hg

We calculate the Schiff moment of the nucleus 199Hg, created by pi-N-N vertices that are odd under parity (P) and time-reversal (T). Our approach, formulated in diagrammatic perturbation theory with important core-polarization diagrams summed to all orders, gives a close approximation to the expectation value of the Schiff operator in the odd-A Hartree-Fock-Bogoliubov ground state generated by a Skyrme interaction and a weak P- and T-odd pion-exchange potential. To assess the uncertainty in the results, we carry out the calculation with several Skyrme interactions (the quality of which we test by checking predictions for the isoscalar-E1 strength distribution in 208Pb) and estimate most of the important diagrams we omit.

preprint2004arXiv

Self-Consistent Description of Multipole Strength in Exotic Nuclei I: Method

We use the canonical Hartree-Fock-Bogoliubov basis to implement a completely self-consistent quasiparticle-random-phase approximation with arbitrary Skyrme energy density functionals and density-dependent pairing functionals. The point of the approach is to accurately describe multipole strength functions in spherical even-even nuclei, including weakly-bound drip-line systems. We describe the method and carefully test its accuracy, particularly in handling spurious modes. To illustrate our approach, we calculate isoscalar and isovector monopole, dipole, and quadrupole strength functions in several Sn isotopes, both in the stable region and at the drip lines.

preprint2002arXiv

What Can Be Learned with a Lead-Based Supernova-Neutrino Detector?

We examine the prospects for using lead as a supernova-neutrino detector by considering the spectrum of electrons produced, and the number of one- and two-neutron events. We show that the electron energy spectrum from charged-current reactions can be used to extract information about the high-temperature component of the neutrino spectrum. Some degree of electron neutrino oscillation is expected in the supernova envelope. We examine the prospects for untangling the signatures of various oscillation scenarios, including, e.g. normal or inverted hierarchies, and different values for the small mixing angle, theta_13.

preprint1999arXiv

Nuclear Octupole Correlations and the Enhancement of Atomic Time-Reversal Violation

We examine the time-reversal-violating nuclear ``Schiff moment'' that induces electric dipole moments in atoms. After presenting a self-contained derivation of the form of the Schiff operator, we show that the distribution of Schiff strength, an important ingredient in the ground-state Schiff moment, is very different from the electric-dipole-strength distribution, with the Schiff moment receiving no strength from the giant dipole resonance in the Goldhaber-Teller model. We then present shell-model calculations in light nuclei that confirm the negligible role of the dipole resonance and show the Schiff strength to be strongly correlated with low-lying octupole strength. Next, we turn to heavy nuclei, examining recent arguments for the strong enhancement of Schiff moments in octupole-deformed nuclei over that of 199Hg, for example. We concur that there is a significant enhancement while pointing to effects neglected in previous work (both in the octupole-deformed nuclides and 199Hg) that may reduce it somewhat, and emphasizing the need for microscopic calculations to resolve the issue. Finally, we show that static octupole deformation is not essential for the development of collective Schiff moments; nuclei with strong octupole vibrations have them as well, and some could be exploited by experiment.

preprint1997arXiv

R-Process Freezeout, Nuclear Deformation, and the Rare-Earth Element Peak

We use network calculations of r-process nucleosynthesis to explore the origin of the peak in the solar r-process abundance distribution near nuclear mass number A = 160. The peak is due to a subtle interplay of nuclear deformation and beta decay, and forms not in the steady phase of the r-process, but only just prior to freezeout, as the free neutrons rapidly disappear. Its existence should therefore help constrain the conditions under which the r-process occurs and freezes out.

preprint1995arXiv

Structure of Pairs in Heavy Weakly-Bound Nuclei

We study the structure of nucleon pairs within a simple model consisting of a square well in three dimensions and a delta-function residual interaction between two weakly-bound particles at the Fermi surface. We include the continuum by enclosing the entire system in a large spherical box. To a good approximation, the continuum can be replaced by a small set of optimally-determined resonance states, suggesting that in many nuclei far from stability it may be possible to incorporate continuum effects within traditional shell-model based approximations.

preprint1995arXiv

What Can Be Learned with an Iodine Solar-Neutrino Detector?

We study the potential benefits of an iodine-based solar-neutrino detector for testing hypotheses that involve neutrino oscillations. We argue that such a detector will have a good chance of distinguishing the two allowed regions of $Δm^2$ -- $\sin^22θ$ parameter space if neutrino conversion is occurring in the sun. It should also be able to detect seasonal variations in the signal due to vacuum oscillations and might be sensitive enough to detect day/night variations due to MSW transitions in the earth. Although it would need to be calibrated, a working iodine detector could be completed before more ambitious projects that seek to accomplish the same things.

preprint1994arXiv

Scatttering of Weakly Interacting Massive Particles from Ge73

We use a "hybrid" method, mixing variationally-determined triaxial nuclear Slater determinants, to calculate the response of Ge73 to hypothetical dark-matter particles such as neutralinos. The method is a hybrid in that rotational invariance is approximately restored prior to variation and then fully restored before the mixing of intrinsic states. We discuss such features of Ge73 as shape coexistence and triaxiality, and their effects on spin-dependent neutralino cross sections. Our calculations yield a satisfactory quadrupole moment and an accurate magnetic moment in this very complicated nucleus, suggesting that the spin structure and thus the axial--vector response to dark matter particles is modeled well.

preprint1993arXiv

Quasielastic neutrino scattering from oxygen and the atmospheric neutrino problem

We examine several phenomena beyond the scope of Fermi-gas models that affect the quasielastic scattering (from oxygen) of neutrinos in the 0.1 -- 3.0 GeV range. These include Coulomb interactions of outgoing protons and leptons, a realistic finite-volume mean field, and the residual nucleon-nucleon interaction. None of these effects are accurately represented in the Monte Carlo simulations used to predict event rates due to $μ$ and $e$ neutrinos from cosmic-ray collisions in the atmosphere. We nevertheless conclude that the neglected physics cannot account for the anomalous $μ$ to $e$ ratio observed at Kamiokande and IMB, and is unlikely to change absolute event rates by more than 10--15\%. We briefly mention other phenomena, still to be investigated in detail, that may produce larger changes.