Search NASA⌕ Search

SEARCH · Search NASA

Results for “ground motion modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

A Theory of Ultrafast Charge Transfer Relaxation with Non-Innocent Solvent Molecules

In this work we revisit the photodynamics of tetracyanoethylene-hexamethylbenzene (TCNE- HMB), the molecular complex studied by Hochstrasser et al. [J. Chem. Phys. 100, 4797–4810, 1994] that has long challenged the applicability of Marcus’s theory of elec- tron transfer for predicting photochemical reactions. Using a novel black-box electronic structure algorithm (time-dependent density functional theory with one double, TD- DFT-1D) to efficiently run molecular dynamics that can treat charge recombination, we run ab initio surface hopping molecular dynamics and confirm that, for a polar solvent, charge recombination rates can be incredibly fast (indeed faster than the sol- vent relaxation time); for non-polar solvents, the rate is much slower. We demonstrate that, although Marcus theory cannot be directly applied, these nonequilibrium (and sometimes incredibly fast) photoexcited dynamics can be effectively explained within a two-state model without any evidence of a transition through a conical intersection. Most importantly, for this paradigmatic model system, we are able to identify two nuclear coordinates of interest (rather than the single coordinate predicted by Marcus or a full set of internal quantum modes studied by Bixon and Jortner): the solvent relaxation in the first shell (that strongly modulates the energies of the charge trans- fer state and differentiates time scales for relaxation) and a nuclear displacement in the TCNE-HMB complex arising from a handful of vibrations that induces non-Born Oppenheimer motion and eventually facilitates an abrupt electronic transition to the ground state. Altogether, these findings suggest a tractable generalization of Marcus theory for future simulations of photochemistry with non-innocent solvent environ- ments in the spirit of a Hamiltonian suggested by Stuchebrukhov (J. Chem. Phys. 107, 3821, 1997).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Observations and seismoacoustic simulations of earthquake-generated infrasound waves in non-epicentral regions

We analysed infrasound waves associated with the Gyeongju earthquake (ML 5.8) that occurred on 2016 September 12, in the southeastern Korean Peninsula. For infrasound wave detection, the Progressive Multichannel Correlation method was applied to the infrasound data set recorded at seven arrays operating in South Korea at epicentral distances ranging from 178 to 472 km. Based on the back-projection method constrained by array-dependent celerity and azimuth deviation models, the source regions were identified in both the epicentral and non-epicentral regions. Remarkably, the non-epicentral secondary sources of this earthquake were located in regions with shallow water depths: (i) the western coastal area in the Yellow Sea and (ii) the shallow ocean basin and bank in the East Sea. The location results obtained from the earthquake could be corroborated through its foreshock (ML 5.1), yielding location results consistent with those of the main shock. The generation of infrasound waves over shallow water depths was fortuitously validated by direct recordings of dominant single-frequency (~0.3 Hz) infrasound waves at close range via temporary sensors near the ocean basin and bank. We interpreted that low-frequency infrasound signals could be generated from interactions among the ocean floor, shallow seawater and atmosphere. We performed numerical simulations of seismoacoustic fields to predict ground motions on the seafloor and acoustic transmission efficiency between the water and air interface. The simulations quantified the energy transfer through different media and clarified our observational results. We found that because this solid Earth‒water‒atmosphere coupled air wave has a relatively low frequency (~0.3 Hz), it can survive propagation over long distances compared with high-frequency infrasound waves generated in inland and mountain regions. In this study, we extend our understanding of water‒atmosphere coupling and the monitoring framework for earthquake-associated non-epicentral infrasound waves, encompassing not only inland ground shaking but also shallow sea regions located far from the epicentre.

Earthquake monitoring and test-ban treaty verifica↗

Linear Seismic Source Equivalents in 3D Nonlinear Models: Effects of Embedded Small-Scale, Near-Source Structures

Gaining a proper understanding of how Earth structure and other near-source properties affect estimates of explosion yield is important to the nonproliferation mission. The yields of explosion sources are often based on seismic moment or waveform amplitudes. Quantifying how the seismic waveforms or estimates of the source characteristics derived from those waveforms are influenced by natural or man-made structures within the near-source region, where the wavefield behaves nonlinearly, is required to understand the full range of uncertainty in those yield estimates. We simulate tamped chemical explosions using a nonlinear, shock physics code and couple the ground motions beyond the elastic radius to a linear elastic, full waveform seismic simulation algorithm through 3D media. In order to isolate the effects of simple small-scale 3D structures on the seismic wavefield and linear seismic source estimates, we embed spheres and cylinders close to the fully- tamped source location within an otherwise homogenous half-space. The 3 m diameters spheres, given their small size compared to the predominate wavelengths investigated, not surprisingly are virtually invisible with only negligible perturbations to the far-field waveforms and resultant seismic source time functions. Similarly, the 11 m diameter basalt sphere has a larger, but still relatively minor impact on the wavefield. However, the 11 m diameter air-filled sphere has the largest impact on both waveforms and the estimated seismic moment of any of the investigated cases with a reduction of ~25% compared to the tamped moment. This significant reduction is likely due in large part to the cavity collapsing from the shock instead of being solely due to diffraction effects . Although the cylinders have the same diameters as the 3 m spheres, their length of interaction with the wavefield produces noticeable changes to the seismic waveforms and estimated source terms with reductions in the peak seismic moment on the order of 10%. Both the cylinders and 11 m diameter spheres generate strong shear waves that appear to emanate from body force sources.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

On the meaning of Berry force for unrestricted systems treated with mean-field electronic structure

We show that the Berry force as computed by an approximate, mean-field electronic structure can be meaningful if properly interpreted. In particular, for a model Hamiltonian representing a molecular system with an even number of electrons interacting via a two-body (Hubbard) interaction and a spin–orbit coupling, we show that a meaningful nonzero Berry force emerges whenever there is spin unrestriction—even though the Hamiltonian is real-valued and formally the on-diagonal single-surface Berry force must be zero. Moreover, if properly applied, this mean-field Berry force yields roughly the correct asymptotic motion for scattering through an avoided crossing. That being said, within the context of a ground-state calculation, several nuances do arise as far interpreting the Berry force correctly, and as a practical matter, the Berry force diverges near the Coulson–Fischer point (which can lead to numerical instabilities). Here, we do not address magnetic fields here.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A three-step model of high harmonic generation using complex classical trajectories

Highlights: • Rigorous trajectory formulation of strong field physics and high harmonic generation. • The Coulomb ground state is represented by an ensemble of complex trajectories. • Trajectories transition seamlessly from the ground state to ionization or recollision. • The ionized wavepacket is reproduced well from the complex trajectories. • Classically allowed and forbidden processes are treated in a single framework. We present a new trajectory formulation of high harmonic generation that treats classically allowed and classically forbidden processes within a single dynamical framework. Complex trajectories orbit the nucleus, producing the stationary Coulomb ground state. When the field is turned on, these complex trajectories continue their motion in the field-dressed Coulomb potential and therefore tunnel ionization, unbound evolution and recollision are described within a single, seamless framework. The new formulation can bring mechanistic understanding to a broad range of strong field physics effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Nucleon-pair coupling scheme in Elliott's SU(3) model

Elliott's SU(3) model is at the basis of the shell-model description of rotational motion in atomic nuclei. Here we demonstrate that SU(3) symmetry can be realized in a truncated shell-model space if constructed in terms of a sufficient number of collective S, D, G,...pairs (i.e., with angular momentum zero, two, four,...) and if the structure of the pairs is optimally determined either by a conjugate-gradient minimization method or from a Hartree-Fock intrinsic state. We illustrate the procedure for six protons and six neutrons in the pf (sdg) shell and exactly reproduce the level energies and electric quadrupole properties of the ground-state rotational band with SDG (SDGI) pairs. The SD-pair approximation without significant renormalization, on the other hand, cannot describe the full SU(3) collectivity. A mapping from Elliott's fermionic SU(3) model to systems with s, d, g,... bosons provides insight into the existence of a decoupled collective subspace in terms of S, D, G,... pairs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Equation-of-motion internally contracted multireference unitary coupled-cluster theory

The accurate computation of excited states remains a challenge in electronic structure theory, especially for systems with a ground state that requires a multireference treatment. In this work, we introduce a novel equation-of-motion (EOM) extension of the internally contracted multireference unitary coupled-cluster framework (ic-MRUCC), termed EOM-ic-MRUCC. EOM-ic-MRUCC follows the transform-then-diagonalize approach, in analogy to its non-unitary counterpart. By employing a projective approach to optimize the ground state, the method retains additive separability and proper scaling with system size. We show that excitation energies are size-intensive if the EOM operator satisfies the “killer” and the projective conditions. Furthermore, we propose to represent changes in the reference state upon electron excitation via projected many-body operators that span the active orbitals and show that the EOM equations formulated in this way are invariant with respect to active orbital rotations. We test the EOM-ic-MRUCC method truncated to single and double excitations by computing the potential energy curves for several excited states of a BeH2 model system, the HF molecule, and water undergoing symmetric dissociation. Across these systems, our method delivers accurate excitation energies and potential energy curves within 5 mE h (∼0.14 eV) from full configuration interaction. Here, we find that truncating the Baker–Campbell–Hausdorff series to fourfold commutators contributes negligible errors (on the order of 10 −5 E h or less), offering a practical route to highly accurate excited-state calculations with reduced computational overhead.

74 ATOMIC AND MOLECULAR PHYSICS↗

Seismic waves in medium with poroelastic/elastic interfaces: a two-dimensional P-SV finite-difference modelling

SUMMARY We present a new methodology of the finite-difference (FD) modelling of seismic wave propagation in a strongly heterogeneous medium composed of poroelastic (P) and (strictly) elastic (E) parts. The medium can include P/P, P/E and E/E material interfaces of arbitrary shapes. The poroelastic part can be with (i) zero resistive friction, (ii) non-zero constant resistive friction or (iii) JKD model of the frequency-dependent permeability and resistive friction. Our FD scheme is capable of subcell resolution: a material interface can have an arbitrary position in the spatial grid. The scheme keeps computational efficiency of the scheme for a smoothly and weakly heterogeneous medium (medium without material interfaces). Numerical tests against independent analytical, semi-analytical and spectral-element methods prove the efficiency and accuracy of our FD modelling. In numerical examples, we indicate effect of the P/E interfaces for the poroelastic medium with a constant resistive friction and medium with the JKD model of the frequency-dependent permeability and resistive friction. We address the 2-D P-SV problem. The approach can be readily extended to the 3-D problem.

58 GEOSCIENCES↗

Electronic and structural properties of RbCe X 2 ( X 2 : O 2 , S 2 , SeS, Se 2 , TeSe, Te 2 )

We report triangular lattice delafossite compounds built from magnetic lanthanide ions are a topic of recent interest due to their frustrated magnetism and realization of quantum disordered magnetic ground states. Here we report the evolution of the structure and electronic ground states of RbCeX 2 compounds, built from a triangular lattice of Ce 3+ ions, upon varying their anion character (X 2 =O 2 , S 2 , SeS, Se 2 , TeSe, Te 2 ). This includes the discovery of a new member of this series, RbCeO 2 , that potentially realizes a quantum disordered ground state analogous to NaYbO 2 . Magnetization and susceptibility measurements reveal that all compounds manifest mean-field antiferromagnetic interactions and, with the exception of the oxide, possess signatures of magnetic correlations onset below 1 K. The crystalline electric field level scheme is explored via neutron scattering and ab initio calculations in order to model the intramultiplet splitting of the J=5/2 multiplet. In addition to the two excited doublets expected within the J=5/2 manifold, we observe one extra local mode present across the sample series. This added mode shifts downward in energy with increasing anion mass and decreasing crystal field strength, suggesting a long-lived anomalous mode endemic to anion motion about the Ce3 + sites.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Following excited-state chemical shifts in molecular ultrafast x-ray photoelectron spectroscopy

The conversion of photon energy into other energetic forms in molecules is accompanied by charge moving on ultrafast timescales. We directly observe the charge motion at a specific site in an electronically excited molecule using time-resolved x-ray photoelectron spectroscopy (TR-XPS). We extend the concept of static chemical shift from conventional XPS by the excited-state chemical shift (ESCS), which is connected to the charge in the framework of a potential model. This allows us to invert TR-XPS spectra to the dynamic charge at a specific atom. We demonstrate the power of TR-XPS by using sulphur 2p-core-electron-emission probing to study the UV-excited dynamics of 2-thiouracil. The method allows us to discover that a major part of the population relaxes to the molecular ground state within 220–250 fs. In addition, a 250-fs oscillation, visible in the kinetic energy of the TR-XPS, reveals a coherent exchange of population among electronic states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Stress-Drop and Source Scaling of the 2019 Ridgecrest, California, Earthquake Sequence

Stress drop, while difficult to measure reliably and at scale, is a key source parameter for understanding the earthquake rupture process and its relationship to strong ground motion. Here, we use a P-wave spectral decomposition approach, designed for large and densely sampled datasets, to measure earthquake stress drop in the region surrounding the 2019 Ridgecrest, California, earthquake sequence. With more than 11,000 measurements of earthquake stress drop in the 20-yr time period from 2000 through 2019, this dataset provides an opportunity to understand how coseismic stress changes and how other geophysical factors relate to the distribution of stress drop and its evolution in space and time. We observe a mild but persistent deviation from self-similar scaling, with larger events having systematically higher stress drops, though this trend depends on the assumption of an omega-square source spectral model. Earthquake stress drop increases with hypocentral depth in this study region, and the Ridgecrest aftershocks tend to have higher stress drops than the pre-event seismicity. This is in part due to their deeper hypocenters. Coherent spatial patterns of stress drop in the aftershock sequence correlate with the slip distribution of the M 7.1 mainshock, whose northwest rupture tip terminated in a long-lived zone of enervated stress drop. Although physical interpretation of these results is complicated by the trade-offs between the timing, depth, and location of these earthquakes, the observations provide new insight into the physics of the earthquake source in an area of renewed seismic activity in southern California.

58 GEOSCIENCES↗

Modeling Subsurface Explosions Recorded on a Distributed Fiber Optic Sensor

Fiber optic distributed acoustic sensors (DAS) are becoming a widely used tool for seismic sensing. Here we examine recordings of two subsurface chemical explosions, DAG-1 and DAG-3, each of which was about one metric ton (TNT equivalent), that were recorded from a helical fiber installed in two boreholes 80 m away from the source location. Several clear phases including the initial P wave, a weak S wave, and a surface reflected P wave are observed on the helical DAS data. We estimate a velocity model using arrival times measured from the fiber. The DAS waveform data were compared with colocated accelerometers at specific depths in both frequency and time domains. Furthermore, the spectra of the DAS data matched spectra estimated from the accelerometer records. Comparisons of observed waveform shape between the accelerometer records and the fiber measurements (strain-rate) show reasonable agreement except for the data near the event depth. The DAS data and the accelerometer agreed in relative amplitudes but we had difficulties in matching absolute amplitudes, possibly due to errors in metadata. Synthetic strain-rate waveforms were calculated using a 2D wavenumber algorithm and matched the waveform shape and relative amplitudes. In general, DAS is effective at recording strong ground motions at high spatial density. Comparison of the synthetic seismograms with observed data indicate that the waveforms are not consistent with a pure isotropic explosion source and that the observed S waves originate from very near the source region.

58 GEOSCIENCES↗

Urban Seismic Site Characterization by Fiber-Optic Seismology

ccurate ground motion prediction requires detailed site effect assessment, but in urban areas where such assessments are most important, geotechnical surveys are difficult to perform, limiting their availability. Distributed acoustic sensing (DAS) offers an appealing alternative by repurposing existing fiber-optic cables, normally employed for telecommunication, as an array of seismic sensors. We present a proof-of-concept demonstration by using DAS to produce high-resolution maps of the shallow subsurface with the Stanford DAS array, California. We describe new methods and their assumptions to assess H/V spectral ratio—a technique widely used to estimate the natural frequency of the soil—and to extract Rayleigh wave dispersion curves from ambient seismic field. These measurements are jointly inverted to provide models of shallow seismic velocities and sediment thicknesses above bedrock in central campus. The good agreement with an independent survey validates the methodology and demonstrates the power of DAS for microzonation.

Geochemistry & Geophysics↗

Monte Carlo Simulation and Reconstruction: Assessment of Myocardial Perfusion Imaging of Tracer Dynamics With Cardiac Motion Due to Deformation and Respiration Using Gamma Camera With Continuous Acquisition

Purpose: Myocardial perfusion imaging (MPI) with single photon emission computed tomography (SPECT) is routinely used for stress testing in nuclear medicine. Recently, our group extended its potential going from 3D visual qualitative image analysis to 4D spatiotemporal reconstruction of dynamically acquired data to capture the time variation of the radiotracer concentration and the estimated myocardial blood flow (MBF) and coronary flow reserve (CFR). However, the quality of reconstructed image is compromised due to cardiac deformation and respiration. The work presented here develops an algorithm that reconstructs the dynamic sequence of separate respiratory and cardiac phases and evaluates the algorithm with data simulated with a Monte Carlo simulation for the continuous image acquisition and processing with a slowly rotating SPECT camera. Methods: A clinically realistic Monte Carlo (MC) simulation is developed using the 4D Extended Cardiac Torso (XCAT) digital phantom with respiratory and cardiac motion to model continuous data acquisition of dynamic cardiac SPECT with slowly rotating gamma cameras by incorporating deformation and displacement of the myocardium due to cardiac and respiratory motion. We extended our previously developed 4D maximum-likelihood expectation-maximization (MLEM) reconstruction algorithm for a data set binned from a continuous list mode (LM) simulation with cardiac and respiratory information. Our spatiotemporal image reconstruction uses splines to explicitly model the temporal change of the tracer for each cardiac and respiratory gate that delineates the myocardial spatial position as the tracer washes in and out. Unlike in a fully list-mode data acquisition and reconstruction the accumulated photons are binned over a specific but very short time interval corresponding to each cardiac and respiratory gate. Reconstruction results are presented showing the dynamics of the tracer in the myocardium as it continuously deforms. These results are then compared with the conventional 4D spatiotemporal reconstruction method that models only the temporal changes of the tracer activity. Mean Stabilized Activity (MSA), signal to noise ratio (SNR) and Bias for the myocardium activities for three different target-to-background ratios (TBRs) are evaluated. Dynamic quantitative indices such as wash-in (K1) and wash-out (k2) rates at each gate were also estimated. Results: The MSA and SNR are higher with higher TBRs while biases were improved with higher TBRs to less than 10%. The correlation between exhalation-inhalation sequence with the ground truth during respiratory cycle was excellent. Our reconstruction method showed better resolved myocardial walls during diastole to systole as compared to the ungated 4D image. Estimated values of K1 and k2 were also consistent with the ground truth. Conclusion: The continuous image acquisition for dynamic scan using conventional two-head gamma cameras can provide valuable information for MPI. Our study demonstrated the viability of using a continuous image acquisition method on a widely used clinical two-head SPECT system. Our reconstruction method showed better resolved myocardial walls during diastole to systole as compared to the ungated 4D image. Precise implementation of reconstruction algorithms, better segmentation techniques by generating images of different tissue types and background activity would improve the feasibility of the method in real clinical environment.

60 APPLIED LIFE SCIENCES↗

Measuring Local Turbulence Along the Optical Path: Multi-Beam Optical Seeing Sensor

Deflection of light along the optical path is a major source of image degradation for ground-based telescopes. Methods have been developed to measure upper atmospheric seeing based on models of the turbulence in the atmosphere, but due to boundary conditions, transmission within telescope enclosures is more complex. The Multi-beam Optical Seeing Sensor (MOSS) directly measures the component of the image quality degradation from inhomogeneity of the index of refraction within the telescope dome. MOSS outputs four near-parallel beams of light that travel along the optical path and are imaged by the telescope’s detector, landing like starlight on the telescope’s focal plane. By using a strobed light source, we can ‘freeze’ the instantaneous index variations transverse to the optical path. This system captures both ‘dome’ and ‘mirror’ seeing. Through plotting the standard deviation of differential motion between pairs of beams, MOSS enables characterization of the length scale of turbulence within the dome. The temporal coherence of temperature gradients can be probed with different pulse lengths, and the spatial coherence by comparing pairs at different separations across the aperture of the telescope. Optical path turbulence measurements, alongside other telemetry metrics, will guide thermal and airflow management to optimize image quality. A MOSS prototype was installed in the 1.2[Formula: see text]m Auxiliary Telescope (AuxTel) at the Vera C. Rubin Observatory in Chile, and preliminary data constrain the optical path turbulence with a lower bound of 1.4 arcsec. The optical path turbulence varied throughout the night of observing.

Astronomical seeing↗

Dynamic Characteristics of Multistory Buildings

In order to evaluate unit stresses or to estimate possible damage in multistory buildings from earthquake or ground motion from underground nuclear explosion it is essential to determine how much of the dynamic response at any level is due to various types of freedom. The effects of shear deformation between floors, joint rotation, over-all flexure, and ground compliance are considered. Joint rotation and over-all flexure are evaluated over a wide range of building characteristics with simply determined indices. Further, a period synthesis concept and a pseudo-stiffness procedure are proposed to enable the determination of periods, mode shapes and stiffnesses with consideration of joint rotation, over-all flexure and ground compliance while performing simple labor-saving analyses for assumed rigid-floor shear buildings. The first three modes of vibration are considered in elastic free vibration. The buildings and models are symmetric in plan without torsional coupling.

42 ENGINEERING↗

Electromagnetic Signatures from the Tidal Tail of a Black Hole—Neutron Star Merger

Black hole - neutron star (BH-NS) mergers are a major target for ground-based gravitational wave observatories. A merger can also produce an electromagnetic counterpart (a kilonova) if it ejects neutron-rich matter that assembles into heavy elements through r-process nucleosynthesis. We study the kilonova signatures of the unbound dynamical ejecta of a BH-NS merger. We take as our initial state the results from a numerical relativity simulation and then use a general relativistic hydrodynamics code to study the evolution of the ejecta with parameterized r-process heating models. The unbound dynamical ejecta is initially a flattened, directed tidal tail largely confined to a plane. Heating from the r-process inflates the ejecta into a more spherical shape and smooths its small-scale structure, though the ejecta retains its bulk directed motion. We calculate the electromagnetic signatures using a 3D radiative transfer code and a parameterized opacity model for lanthanide-rich matter. The light curve varies with viewing angle because of two effects: asphericity results in brighter emission for orientations with larger projected areas, while Doppler boosting results in brighter emission for viewing angles more aligned with the direction of bulk motion. For typical r-process heating rates, the peak bolometric luminosity varies by a factor of ~3 with orientation while the peak in the optical bands varies by ~3 magnitudes. The spectrum is blueshifted at viewing angles along the bulk motion, which increases the V-band peak magnitude to ~-14 despite the lanthanide-rich composition.

79 ASTRONOMY AND ASTROPHYSICS↗

Long-lived Bell states in an array of optical clock qubits

The generation of long-lived entanglement in optical atomic clocks is one of the main goals of quantum metrology. Arrays of neutral atoms, where Rydberg-based interactions may generate entanglement between individually controlled and resolved atoms, constitute a promising quantum platform to achieve this. Here we leverage the programmable state preparation afforded by optical tweezers and the efficient strong confinement of a three-dimensional optical lattice to prepare an ensemble of strontium-atom pairs in their motional ground state. We engineer global single-qubit gates on the optical clock transition and two-qubit entangling gates via adiabatic Rydberg dressing, enabling the generation of Bell states with a state-preparation-and-measurement-corrected fidelity of 92.8(2.0)% (87.1(1.6)% without state-preparation-and-measurement correction). For use in quantum metrology, it is furthermore critical that the resulting entanglement be long lived; we find that the coherence of the Bell state has a lifetime of 4.2(6) s via parity correlations and simultaneous comparisons between entangled and unentangled ensembles. Such long-lived Bell states can be useful for enhancing metrological stability and bandwidth. In the future, atomic rearrangement will enable the implementation of many-qubit gates and cluster state generation, as well as explorations of the transverse field Ising model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗