Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Displacement”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Vertical instability forecasting and controllability assessment of multi-device tokamak plasmas in DECAF with data-driven optimization

Abstract Reliable vertical position control will be an essential element of any future tokamak-based fusion power plant in order to reduce disruptions and maximize performance. We investigate methods to improve vertical controllability boundary determination in plasma operational space and demonstrate a data-driven approach based on direct pseudoinversion of operational space data that is rigorously quantitative, applicable in real-time plasma control systems, and physically intuitive to interpret. Applied to historical shot data from entire run campaigns on the MAST-U, KSTAR, and NSTX tokamaks, this approach, implemented in DECAF, improves vertical displacement event identification accuracy to 98.9%–100%. Further, we explore the application of a physics-based vertical stability metric as an early warning forecaster for vertical displacement events. The development of a linear surrogate model for the plasma current density profile, with a coefficient of determination of 0.992 on the training dataset, enables potential employment of this forecaster in real-time. The application of this approach on historical data from the MAST-U MU02 campaign yields a forecaster with 62.6% accuracy, indicating promise for this method when further refined and potentially coupled with other stability metrics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Accurate predictions from small boxes: variance suppression via the Zel'dovich approximation

Simulations have become an indispensable tool for accurate modelling of observables measured in galaxy surveys, but can be expensive if very large dynamic range in scale is required. We describe how to combine Lagrangian perturbation theory models with N-body simulations to reduce the effects of finite computational volume in the prediction of ensemble average properties in the simulations within the context of control variates. In particular we use the fact that Zel'dovich displacements, computed during initial condition generation for any simulation, correlate strongly with the final density field. Since all the correlators of biased tracers can be computed with arbitrary precision for these displacements, pairing the Zel'dovich `simulation' with the N-body realization allows hundredfold reductions in sample variance for power spectrum or correlation function estimation. Additionally, Zel'dovich control variates can accurately extend matter or tracer field emulators to larger scales than previously possible, as well as improving measurements of statistics in simulations which are inherently limited to small volumes, such as hydrodynamical simulations of galaxy formation and reionization.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

An atom probe tomography and inventory calculation examination of second phase precipitates in neutron irradiated single crystal tungsten

Tungsten is the prime candidate material for the plasma facing divertor system and first wall armor of future nuclear fusion reactors. However, the understanding of the microstructural and chemical evolution of pure tungsten under neutron irradiation is relatively unknown, in part due to a lack of experimental data on this topic. Here, single crystal tungsten has been irradiated in the High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory at a nominal temperature of 700–800 °C to damage levels of 0.1 and 1.8 displacements per atom (dpa). Inventory calculations of the neutron irradiation experiments have been used to track isotope generation and decay to inform atom probe tomography (APT) results in the determination of the transmutation-induced precipitate compositions. Furthermore, APT crystallography has been used to show the relationship between precipitates and matrix. The composition of the precipitates is shown to progress towards that of the σ-phase at the highest dose studied, with those precipitates lying along crystallographic planes similar to those of displacement-induced dislocations. This work also sets the framework for APT studies of materials that contain isotopic ratios far from those observed in the natural state.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Influence of triangularity on the plasma response to resonant magnetic perturbations

The effect of the axisymmetric plasma shape on the non-axisymmetric plasma response to resonant magnetic perturbations is investigated in experiment and modeling for the DIII-D, EAST and ASDEX Upgrade tokamaks. Systematically modeling the effect of the triangularity, whilst keeping other equilibrium quantities largely unchanged, reveals that the plasma response is strongly suppressed at high triangularity compared to that at low triangularity. This is validated through targeted comparison with experiments at DIII-D, EAST, and ASDEX Upgrade. DIII-D and EAST magnetic measurements are used to validate simulations, while at ASDEX Upgrade the plasma edge displacement is measured. Both experiments and modeling find a reduced magnetic plasma response on the high-field side at high triangularity across devices. Multi-modal analysis of the simulation results extracts the mode structure and applied spectrum dependence of each mode. The amplitude of the dominant mode reveals similar trends with the edge resonance and radial displacement near the X-point, which suggests that the multi-mode plasma response provides another way to understand the edge localized mode (ELM) control physics. The plasma response is strongly reduced at high triangularity compared to that at low triangularity, which implies different ELM control effects as shaping is varied. Furthermore, these findings indicate that the plasma shape should be taken into consideration when designing an RMP-ELM control strategy in experiment, and that predictive plasma response calculations can be used to maximize access to RMP-ELM control in future devices by maximizing the coupling between coils and the plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Sensitivity of magnetic islands in permanent magnet stellarators using the gradient and Hessian methods

Stellarator plasmas are known to be very sensitive to perturbations in the magnetic field. The permanent magnet stellarator was in part developed as a solution to high machining tolerances placed on the shape properties of electromagnetic coils in traditional stellarators. However, as a consequence of this high sensitivity to the field structure, sensitivities of permanent magnet stellarator plasmas to perturbations of permanent magnet properties must necessarily be well-understood. The gradient and Hessian matrix methods have been previously demonstrated to be useful sensitivity analysis methods for modular coils. We apply these two methods to the study of island width sensitivities in both the MUSE and PM4STELL permanent magnet stellarator projects. These sensitivity methods were used to determine the relative impacts of permanent magnet parameter perturbations on island widths in the vacuum field approximation of both stellarator equilibria. The square of resonant magnetic field perturbation is used here as a proxy for island width. In particular, gradients of magnetizations of individual magnets were examined in MUSE, as well as gradients of magnet group displacements informed by device design. Three different forms of permanent magnet magnetization perturbations are investigated for MUSE, and the flux surface response to perturbations is demonstrated. The Hessian matrix method is applied to PM4STELL, illustrating the sensitivity of dominant island widths to displacements of toroidal wedge structures. These methods allow for selective direction of experimental resources toward regions of heightened sensitivity, while constraints on less impactful permanent magnet parameters can be relaxed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Error analysis of low-fidelity models for wake steering based on field measurements

The observations collected by two scanning lidars deployed on the roof of a 2.8-MW turbine undergoing a series of imposed yaw offsets are analyzed. The wake lateral displacement detected by the rear-facing lidar correlates well with the yaw offset sensed by the forward-facing lidar. We find that the high-frequency part of the yaw offset signal is connected to wake meandering, whereas the low frequency component is a good predictor for wake displacement due to yaw misalignment. Conditionally averaged wake velocity data for different yaw offsets are used as benchmarks for the validation of a linearized Reynolds-averaged Navier-Stokes and an empirical wake model. A mean error as low as 2% and a good prediction of the wake trajectory are achieved, provided that the wake recovery rate matches the observations.

17 WIND ENERGY↗
DOE OSTI2020

Measurement of space charge effects in the MicroBooNE LArTPC using cosmic muons

Large liquid argon time projection chambers (LArTPCs), especially those operating near the surface, are susceptible to space charge effects. In the context of LArTPCs, the space charge effect is the build-up of slow-moving positive ions in the detector primarily due to ionization from cosmic rays, leading to a distortion of the electric field within the detector. This effect leads to a displacement in the reconstructed position of signal ionization electrons in LArTPC detectors ("spatial distortions"), as well as to variations in the amount of electron-ion recombination experienced by ionization throughout the volume of the TPC. We present techniques that can be used to measure and correct for space charge effects in large LArTPCs by making use of cosmic muons, including the use of track pairs to unambiguously pin down spatial distortions in three dimensions. Finally, the performance of these calibration techniques are studied using both Monte Carlo simulation and MicroBooNE data, utilizing a UV laser system as a means to estimate the systematic bias associated with the calibration methodology.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Characterization of proton-induced damage in thick, p-channel skipper-CCDs

Here, in this work, we characterize the radiation-induced damage in two thick, p-channel skipper-CCDs irradiated unbiased and at room temperature with 217-MeV protons. We evaluate the overall performance of the sensors and demonstrate their single-electron/single-photon sensitivity after receiving a fluence on the order of 10 10 protons/cm 2 . Using the pocket-pumping technique, we quantify and characterize the proton-induced defects from displacement damage. We report an overall trap density of 0.134 traps/pixel for a displacement damage dose of 2.3 × 10 7 MeV/g. Three main proton-induced trap species were identified, V 2 , C i O i and V n O m , and their characteristic trap energies and cross sections were extracted. We found that while divacancies are the most common proton-induced defects, C i O i defects have a greater impact on charge integrity at typical operating temperatures because their emission-time constants are comparable or larger than typical readout times. To estimate ionization damage, we measure the characteristic output transistor curves. We found no threshold voltage shifts after irradiation. Our results highlight the potential of skipper-CCDs for applications requiring high-radiation tolerance and can be used to find the operating conditions in which effects of radiation-induced damage are mitigated.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Stability of trapped solutions of a nonlinear Schrödinger equation with a nonlocal nonlinear self-interaction potential

This work focuses on the study of the stability of trapped soliton-like solutions of a (1 + 1)-dimensional nonlinear Schrödinger equation (NLSE) in a nonlocal, nonlinear, self-interaction potential of the form [|Ψ(x,t)| 2 +|Ψ(-x,t)| 2 ] κ where κ is an arbitrary nonlinearity parameter. Although the system with κ = 1 (i.e. fully integrable case) was first reported by Yang (2018 Phys. Rev. E 98 042202), here in the present work, we extend this model to the one in which κ is arbitrary. This allows us to compare the stability properties of the now trapped solutions to previously found solutions of the more usual NLSE with κ ≠ 1 which are moving soliton solutions. We show that there is a simple, one-component, nonlocal Lagrangian and corresponding action governing the dynamics of the system. Using a collective coordinate method derived from the action as well as assuming the validity of Derrick's theorem, we find that these trapped solutions are stable for 0 < κ < 2 and unstable when κ > 2. At the critical value of κ, i.e. κ = 2, the solution can either collapse or blowup linearly in time when q 0 = 0, where q 0 is the center of the initial density ρ(x, t = 0) = ψ*ψ of the solution. For q 0 ≠ 0 the displaced solution collapses. When κ > 2 initial small displacements from the origin also lead to collapse of the wave function. This phenomenon is not seen in the usual NLSE.

collective coordinates↗
DOE OSTI2020

A deep neural network to search for new long-lived particles decaying to jets

A tagging algorithm to identify jets that are significantly displaced from the proton-proton (pp) collision region in the CMS detector at the LHC is presented. Displaced jets can arise from the decays of long-lived particles (LLPs), which are predicted by several theoretical extensions of the standard model. The tagger is a multiclass classifier based on a deep neural network, which is parameterised according to the proper decay length $\mathrm{c}\tau_0$ of the LLP. A novel scheme is defined to reliably label jets from LLP decays for supervised learning. Samples of pp collision data, recorded by the CMS detector at a centre-of-mass energy of 13 TeV, and simulated events are used to train the neural network. Domain adaptation by backward propagation is performed to improve the simulation modelling of the jet class probability distributions observed in pp collision data. The potential performance of the tagger is demonstrated with a search for long-lived gluinos, a manifestation of split supersymmetric models. The tagger provides a rejection factor of 10 000 for jets from standard model processes, while maintaining an LLP jet tagging efficiency of 30-80% for gluinos with 1 mm $\leq$ $c\tau_0$ $\leq$ 10 m. The expected coverage of the parameter space for split supersymmetry is presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Combining different 3-D global and regional seismic wave propagation solvers towards box tomography in the deep Earth

SUMMARY In previous publications, we presented a general framework, which we called ‘box tomography’, that allows the coupling of any two different numerical seismic wave propagation solvers, respectively outside and inside a target region, or ‘box’. The goal of such hybrid wavefield computations is to reduce the cost of computations in the context of full-waveform inversion for structure within the target region, when sources and/or receivers are located at large distances from the box. Previously, we had demonstrated this approach with sources and receivers outside the target region in a 2-D acoustic spherical earth model, and demonstrated and applied this methodology in the 3-D spherical elastic Earth in a continental scale inversion in which all stations were inside the target region. Here we extend the implementation of the approach to the case of a 3-D global elastic earth model in the case where both sources and stations are outside the box. We couple a global 3-D solver, SPECFEM3D_GLOBE, for the computation of the wavefield and Green’s functions in a reference 3-D model, with a regional 3-D solver, RegSEM, for the computation of the wavefield within the box, by means of time-reversal mirrors. We briefly review key theoretical aspects, showing in particular how only the displacement is needed to be stored at the boundary of the box. We provide details of the practical implementation, including the geometrical design of the mirrors, how we deal with different sizes of meshes in the two solvers, and how we address memory-saving through the use of B-spline compression of the recorded wavefield on the mirror. The proposed approach is numerically efficient but also versatile, since adapting it to other solvers is straightforward and does not require any changes in the solver codes themselves, as long as the displacement can be recovered at any point in time and space. We present benchmarks of the hybrid computations against direct computations of the wavefield between a source and an array of stations in a realistic geometry centred in the Yellowstone region, with and without a hypothetical plume within the ‘box’, and with a 1-D or a 3-D background model, down to a period of 20 s. The ultimate goal of this development is for applications in the context of imaging of remote target regions in the deep mantle, such as, for example, Ultra Low Velocity Zones.

Geochemistry & Geophysics↗
DOE OSTI2020

Reconstruction with velocities

ABSTRACT Reconstruction is becoming a crucial procedure of galaxy clustering analysis for future spectroscopic redshift surveys to obtain subper cent level measurement of the baryon acoustic oscillation scale. Most reconstruction algorithms rely on an estimation of the displacement field from the observed galaxy distribution. However, the displacement reconstruction degrades near the survey boundary due to incomplete data and the boundary effects extend to ${\sim}100\, \mathrm{Mpc}/h$ within the interior of the survey volume. We study the possibility of using radial velocities measured from the cosmic microwave background observation through the kinematic Sunyaev–Zeldovich effect to improve performance near the boundary. We find that the boundary effect can be reduced to ${\sim}30-40\, \mathrm{Mpc}/h$ with the velocity information from Simons Observatory. This is especially helpful for dense low redshift surveys where the volume is relatively small and a large fraction of total volume is affected by the boundary.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Achievable Rates for Concatenated Square Gottesman-Kitaev-Preskill Codes

The Gottesman-Kitaev-Preskill (GKP) codes are known to achieve optimal rates under displacement noise and pure-loss channels, which establishes theoretical foundations for its optimality. However, such optimal rates are only known to be achieved at a discrete set of noise strengths with the current self-dual symplectic lattice construction. In this work, we develop a new coding strategy using concatenated continuous variable-discrete variable encodings to go beyond past results and establish GKP’s optimal rate over all noise strengths. In particular, for displacement noise, the rate is obtained through a constructive approach by concatenating GKP codes with a quantum polar code and analog decoding. For a pure-loss channel, we prove the existence of capacity-achieving GKP codes through a random coding approach. These results highlight the capability of concatenation-based GKP codes and provides new methods for constructing good GKP lattices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Importance of finite-size corrections for accurate ab initio modeling of carrier capture at semiconductor defects: A case study of substitutional C N in GaN

In ab initio studies of carrier-capture processes in defective semiconductor materials, the single-effective-mode formalism and the static-coupling approximation have become the predominant theoretical approaches for determining carrier-capture coefficients. The single-mode formalism relies on accurate nonequilibrium defect energies obtained from density-functional theory (DFT), where required inputs are a series of configurationally displaced, defect-containing supercells obtained using an interpolative ansatz, and where the DFT outputs are corresponding total energies that have traditionally been postprocessed using a long-established ground-state formulation of finite-size corrections and defect-formation energies. This formulation remains commonly used even though the defects that form a configuration-coordinate (CC) diagram typically exist as structures that are displaced from the ground state. To remedy this inconsistency, Kumagai has recently proposed novel methods for implementing finite-size corrections specifically intended for DFT calculations of the defect energies used to construct CC diagrams and implement the single-mode formalism [Y. Kumagai, Phys. Rev. B 107, L220101 (2023)]. Kumagai's approach builds on the latest finite-size-correction methods introduced to describe vertical charge-state transitions for charge-localizing point defects in semiconductors and insulators [T. Gake et al., Phys. Rev. B 101, 020102 (2020); S. Falletta et al., Phys. Rev. B 102, 041115 (2020)]. The newly identified finite-size artifact treated in these studies is the polarization charge induced on a configurationally frozen defect and its subsequent interaction with a vertical transition in charge state. In this work, we evaluate Kumagai's proposed methodology by applying it in a high-precision DFT study of carrier capture by substitutional C N in GaN, a well-characterized and technologically relevant defect and material. We have rigorously calculated C N defect energies across various supercell sizes for each defect configuration and charge state on the hole-capture CC diagram of C N (𝑞=−1), enabling a direct comparison of the slopes of the defect energies versus inverse cell size with those predicted by Kumagai. The most consequential prediction of Kumagai's method is that these slopes distinctly vary as the square of the linear-interpolation parameter used to construct the nonequilibrium defect configurations. Our results quantitatively support this prediction. Moreover, with these new finite-size corrections and multiple-cell-size DFT calculations in place, we find that the classical energy barrier for hole capture by C N (𝑞=−1) in GaN decreases to 0.092–0.127 eV. This finding confirms the recent ≈ 0.1 eV prediction of Reshchikov based on the weak temperature dependence for hole capture observed in photoluminescence experiments [M. A. Reshchikov, J. Appl. Phys. 129, 121101 (2021)]. These results stand in stark contrast to previously calculated barriers of 0.486 and 0.73 eV, which also used the single-mode formalism but were obtained by instead using ground-state-based finite-size corrections. Our reduced classical barrier for capture increases the temperature-dependent hole-capture coefficient of a C N (𝑞=−1) defect by more than two to four orders of magnitude for temperatures of 100–600 K, compared to the previous 0.486 eV results. While other defects may not be as dramatically affected as here, we suggest that incorporating proper finite-size corrections for the vertical-transition-like states embedded within CC diagrams is an essential, yet previously unrecognized, component of accurate modeling of carrier-capture when using the single-effective-mode formalism.

dielectric properties↗
DOE OSTI2021

Mass enhancement in 3 d and s - p perovskites from symmetry breaking

In some d-electron oxides, the measured effective mass m* exptl has long been known to be significantly larger than the model effective mass m* model deduced from mean-field band theory, i.e., m* exptl =βm* model , where β>1 is the “mass-enhancement” or “mass-renormalization” factor. Previous applications of density functional theory (DFT), based on a symmetry-restricted structure with the smallest number of possible magnetic, orbital, and structural degrees of freedom, missed such mass enhancement. This fact has been taken as evidence of strong electronic correlation, often described via the symmetry-restricted dynamic mean-field approach of the many-body theory, being the exclusive enabling physics. This paper uses instead a static density functional approach that does not restrict positional or spin degrees of freedom (symmetry-broken structures). This approach analyzes the contributions of different symmetry-broken modalities to mass enhancement for a few nominally highly correlated d-electron perovskites as well as the nominally uncorrelated, closed-shell s-p bonding perovskites. It shows that the energy-lowering symmetry-broken spin effects (e.g., nonzero local moment in the paramagnetic phase) and structural effects (e.g., atomic displacement) as described in mean-field DFT already manifest mass enhancement for both electrons and holes in a range of d-electron perovskites SrVO 3 , SrTiO 3 , BaTiO 3 , and LaMnO 3 , as well as p-electron perovskites CsPbI 3 and SrBiO 3 , including both metals (SrVO 3 ) and insulators (the rest). This is revealed only when enlarged unit cells of the same parent global symmetry, which are large enough to allow for symmetry-breaking distortions and concomitant variations in spin order, are explored for their ability to lower the total energy. Positional symmetry breaking that leads to mass enhancement includes octahedral rotation in halide perovskites such as CsPbI 3 , Jahn-Teller-like Q${+}\atop{2}$ distortion in LaMnO 3 , and bond disproportionation in SrBiO 3 , while magnetic symmetry breaking resulting in mass enhancement includes the formation of a distribution of local moments in SrVO3 that averages to zero in the paramagnetic phase. Not all symmetry breaking leads to significant mass enhancement, e.g., the rather small octahedral rotations in the nearly perfectly cubic SrTiO 3 cause negligible mass enhancement, as do the paraelectric displacements in BaTiO 3 . In principle, by ergodicity, the two descriptions, i.e., the symmetry-restricted dynamic approach with a single time-fluctuating magnetic moment and the symmetry-broken mean-field approach with a static distribution of spatially fluctuated local moments, are not mutually exclusive but are a choice of representation and consequently, a choice of computational efficiency. In approximate implementations, the symmetry-broken mean-field approach appears to remove much of what was strong correlation in dynamically correlated symmetry-restricted solutions, leaving smaller (“weak”) residual correlation with respect to the exact solution.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Prevalence of pretransition disordering in the rutile-to- CaCl 2 phase transition of GeO 2

WE report the ability to tailor a material's electronic properties using density driven disordering has emerged as a powerful route to materials design. The observation of anomalous structural and electronic behavior in the rutile to CaCl 2 phase transition in SnO 2 led to the prediction that such behavior is inherent to all oxides experiencing such a phase transition sequence [Smith et al., J. Phys. Chem. Lett. 10, 5351 (2019)]. Here, the ultrawide band gap semiconductor GeO 2 is confirmed to exhibit anomalous behavior during the rutile to CaCl 2 phase transition. A phase pure rutile GeO 2 sample synthesized under high-pressure, high-temperature conditions is probed using synchrotron diffraction and x-ray and optical spectroscopy under high pressure conditions. Density functional theory calculations show that the enthalpic barrier to displacing an oxygen along the B 1g librational mode decreases with pressure leading up to the rutile to CaCl 2 phase transition. The band structure of the distorted state shows that such oxygen displacements form small polarons.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Reentrant superconductivity through a quantum Lifshitz transition in twisted trilayer graphene

A series of recent experiments has demonstrated robust superconductivity in magic-angle twisted trilayer graphene (TTG). In particular, a recent work by Cao et al. [Nature (London) 595, 526 (2021)] studies the behavior of the superconductor in an in-plane magnetic field and an out-of-plane displacement field, finding that the superconductor is unlikely to have purely spin-singlet pairing. Herein, this work also finds that at high magnetic fields and a smaller range of dopings and displacement fields, the superconductor undergoes a transition to a distinct field-induced superconducting state. Inspired by these results, we develop an understanding of the superconductivity in TTG using a combination of phenomenological reasoning and microscopic theory. We describe the role that an in-plane field plays in TTG, and we use this understanding to argue that the reentrant transition may be associated with a quantum Lifshitz phase transition, with the high-field phase possessing finite-momentum pairing. We argue that the superconductor is likely to involve a superposition of singlet and triplet pairing, and we describe the structure of the normal state. We also draw lessons for twisted bilayer graphene (TBG), and we explain the differences in the phenomenology with TTG despite their close microscopic relationship. We propose that a singlet-triplet superposition is realized in the TBG superconductor as well, and that the ν = –2 correlated insulator may be a time-reversal protected Z 2 topological insulator obtained through spontaneous spin symmetry breaking.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Optically and thermally driven huge lattice orbital and spin angular momenta from spinning fullerenes

Lattice vibration in solids may carry angular momentum. But unlike the intrinsic spin of electrons, the lattice vibration is rarely rotational. To induce angular momentum, one needs to find a material that can accommodate a twisted normal mode, two orthogonal modes, or excitation of magnons. If excitation is too strong, one may exceed the Lindemann limit, so the material melts. Therefore these methods are not ideal. Here, we theoretically propose a route to phonon angular momentum in a molecular crystal C 60 . We find that a single laser pulse is able to inject a significant amount of angular momentum to C 60 , and the momentum transfer is helicity dependent. Changing from right-circularly polarized light to left-circularly polarized light switches the direction of phonon angular momentum. On the ultrafast timescale, the orbital angular momentum change closely resembles the displacive excitation of coherent phonons, with a cosine-function dependence on time, different from the spin counterpart. Atomic displacements, even under strong laser excitation, remain far below the Lindemann criterion. Under thermal excitation, spinning C 60 even at room temperature generates a huge angular momentum close to several hundred h¯. Our finding opens the door to a large group of fullerenes, from C 60 , C 70 to their endohedral derivatives, where angular momentum can be generated through light or temperature. Furthermore, this paves the way to the phononic control electronic spin and harvesting thermal energy through phonon angular momentum.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗