Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linearity preservation”

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 73 records · Page 4

Non-invertible defects from the Conway SCFT to K3 sigma models. Part I. General results

We initiate the study of supersymmetry-preserving topological defect lines (TDLs) in the Conway moonshine module V f♮ . We show that the tensor category of such defects, under suitable assumptions, admits a surjective but non-injective ring homomorphism into the ring of $\mathbb{Z}$-linear maps of the Leech lattice into itself. This puts strong constraints on possible defects and their quantum dimensions. We describe a simple construction of non-invertible TDLs from orbifolds of holomorphic (super)vertex operator algebras, which yields non-trivial examples of TDLs satisfying our main theorem. We conjecture a correspondence between four-plane-preserving TDLs in V f♮ and supersymmetry-preserving TDLs in K3 non-linear sigma models, which extends the correspondence between symmetry groups to the level of tensor category symmetry. We establish evidence for this conjecture by constructing non-invertible TDLs in special K3 non-linear sigma models.

Conformal Field Models in String Theory↗

Elastocaloric signatures of symmetric and antisymmetric strain-tuning of quadrupolar and magnetic phases in DyB 2 C 2

The adiabatic elastocaloric effect measures the temperature change of a given system with strain and provides a thermodynamic probe of the entropic landscape in the temperature-strain space. Here, we demonstrate that the DC bias strain-dependence of AC elastocaloric effect allows decomposition of the latter into symmetric (rotation-symmetry-preserving) and antisymmetric (rotation-symmetry-breaking) strain channels, using a tetragonal f -electron intermetallic DyB 2 C 2 —whose antiferroquadrupolar order breaks local fourfold rotational symmetries while globally remaining tetragonal—as a showcase example. We capture the strain evolution of its quadrupolar and magnetic phase transitions using both singularities in the elastocaloric coefficient and its jumps at the transitions, and the latter we show follows a modified Ehrenfest relation. We find that antisymmetric strain couples to the underlying order parameter in a biquadratic (linear-quadratic) manner in the antiferroquadrupolar (canted antiferromagnetic) phase, which are attributed to a preserved (broken) global tetragonal symmetry, respectively. The broken tetragonal symmetry in the magnetic phase is further evidenced by elastocaloric strain-hysteresis and optical birefringence. Additionally, within the staggered quadrupolar order, the observed elastocaloric response reflects a quadratic increase of entropy with antisymmetric strain, analogous to the role magnetic field plays for Ising antiferromagnetic orders by promoting pseudospin flips. Our results demonstrate AC elastocaloric effect as a compact and incisive thermodynamic probe into the coupling between electronic degrees of freedom and strain in free energy, which holds the potential for investigating and understanding the symmetry of a wide variety of ordered phases in broader classes of quantum materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Nevada Higher Education Benefits from US Department of Energy (DOE) Environmental Management (EM) Nevada Program's Transfer of Geologic Samples - 20442

The Nevada National Security Site (NNSS), formerly the Nevada Test Site, was the location of 100 historic atmospheric and 828 historic underground nuclear tests from 1951 to 1992. Related to this historic nuclear testing, geologic and hydrologic studies of the site were conducted utilizing the skills and the expertise of the U.S. Geological Survey (USGS) and national laboratories (Los Alamos National Laboratory, Lawrence Livermore National Laboratory), who were principal in leading the development of weapons and were responsible for specific underground testing programs. The subsurface samples and technical/scientific data associated with NNSS geologic studies (past and present) are preserved and stored at the USGS Mercury Core Library and Data Center located at the NNSS. Currently, the facility stores over 2,000,000 linear feet of cores and cuttings from more than 2,600 drill holes that can be accessed for study. Most of the samples and historic geologic work was focused directly on subsurface geologic settings that relate directly to historic underground nuclear testing. Since 1992, the United States has observed a unilateral moratorium on full-scale nuclear testing, and the U.S. Department of Energy (DOE) Environmental Management (EM) Nevada Program is now responsible for hydrogeologic characterization of the potential impacts to the natural groundwater systems that may have resulted from historic underground nuclear tests. As part of this effort, the EM Nevada Program Underground Test Area (UGTA) Activity has drilled and completed over 58 deep (2,000 - 7,000 ft.) characterization wells, totaling in excess of 170,000 linear feet of cuttings and core samples. Cutting samples were collected as triplicate samples for each respective depth interval, to account for potential later nondestructive/destructive analysis and to preserve samples for regulatory purposes. Recently, it was recognized by the EM Nevada Program that opportunities may exist to reduce the cost and floor space required for the storage of geologic samples at the Mercury Core Library without impacting the integrity and representative nature of the samples necessary for project execution. An initiative was sponsored by EM Nevada Program to evaluate several options: 1) disposal of a portion of cuttings and cores in a land fill setting; 2) reduction in sample volume through skeletonizing core and cuttings from wells to a representative but much smaller number of samples; and 3) solicit potential interest in the academic community where other geoscientists could freely access the samples for studies. The Nevada state university system through the Nevada Bureau of Mines and Geology responded positively to the opportunity to receive these NNSS samples from the EM Nevada Program. In June 2019, approximately 17,000 geologic samples, representing greater than 170,000 linear feet of drilling and weighing over 20,000 pounds, were shipped from the NNSS to the Great Basin Science Sample and Records Library located in Reno, Nevada. The benefits resulting from this transfer were realized by the EM Nevada Program in terms of cost and space savings for geologic sample storage at the NNSS. The transfer did not impact the EM Nevada Program mission as access to representative geologic samples for regulatory and scientific purposes was preserved. Moreover, the Nevada Bureau of Mines and Geology acquired a significant resource of subsurface geologic samples and supporting technical data to support academic and scientific studies in a complex volcanic setting in southwest Nevada. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A flexible gyro-fluid system of equations

Gyro-fluid equations are velocity space moments of the gyrokinetic equations. Special gyro-Landau-fluid closures have been developed that include the damping due to kinetic resonances by fitting to the collisionless local plasma response functions. This damping allows for accurate linear eigenmodes to be computed with a relatively low number of velocity space moments compared to the number of velocity quadrature points in gyrokinetic codes. However, none of the published gyro-Landau-fluid closure schemes considers the Onsager symmetries of the resulting quasi-linear fluxes as a constraint. Onsager symmetry guarantees that the matrix of diffusivities is positive definite, an important property for the numerical stability of a transport solver. A two-parameter real closure for improving the accuracy of low-resolution gyro-fluid equations, which preserves the Onsager symmetry and allows higher velocity space moments, is presented in this paper. The new linear gyro-fluid system (GFS) is used to extend the TGLF quasi-linear transport model so that it can compute the energy and momentum fluxes due to parallel magnetic fluctuations, completing the transport matrix. The GFS equations do not use a bounce average approximation. The GFS equations are fully electromagnetic with general flux surface magnetic geometry, pitch angle scattering for electron collisions, and subsonic equilibrium toroidal rotation. Using GFS eigenmodes in the quasi-linear TGLF model will be shown to yield a more accurate match to fluxes computed by CGYRO turbulence simulations. In conclusion, prospects for future applications of a quasi-linear theory to new plasma transport regimes and magnetic confinement devices in addition to tokamaks are opened by the flexibility of the GFS eigensolver.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Six-dimensional matching of intense beam with linear accelerating structure

Beam matching is a common technique that is routinely employed in accelerator design with the aim of minimizing beam losses and preservation of beam brightness. Despite being widely used, a full theoretical understanding of beam matching in 6D remains elusive. In this work, we present an analytical treatment of 6D beam matching of a high-intensity beam onto an RF structure. We begin our analysis within the framework of a linear model, and apply the averaging method to a set of 3D beam envelope equations. Accordingly, we obtain a matched solution that is comprised of smoothed envelopes and periodic terms, describing envelope oscillations with the period of the focusing structure. We then consider the nonlinear regime, where the beam size is comparable with the separatrix size. Stating with a Hamiltonian analysis in 6D phase space, we attain a self-consistent beam profile and show that it is significantly different from the commonly used ellipsoidal shape. Subsequently, we analyze the special case of an equilibrium with equal space charge depression between all degrees of freedom. Comparison of beam dynamics for equipartitioned, equal space charge depression, and equal emittances beams is given. Finally, we present experimental results on beam matching in the LANSCE linac.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Is Mini Beam Ready for Human Trials? Results of Randomized Study of Treating De-Novo Brain Tumors in Canines Using Linear Accelerator Generated Mini Beams

The main challenge in treating malignant brain neoplasms lies in eradicating the tumor while minimizing treatmentrelated damage. Conventional radiation treatments are associated with considerable side effects. Synchrotron generated micro-beam radiation (SMBRT) has shown to preserve brain architecture while killing tumor cells, however physical characteristics and limited facility access restrict its use. We have created a new clinical device which produces mini beams on a linear accelerator, to provide a new type of treatment called mini-beam radiation therapy (MBRT). The objective of this study is to compare the treatment outcomes of linear accelerator based MBRT versus standard radiation treatment (SRT), to evaluate the tumor response and the treatment-related changes in the normal brain with respect to each treatment type. Pet dogs with de-novo brain tumors were accrued for treatment. Dogs were randomized between standard fractionated stereotactic (9 Gy in 3 fractions) radiation treatment vs. a single fraction of MBRT (26 Gy mean dose). Dogs were monitored after treatment for clinical assessment and imaging. When the dogs were euthanized, a veterinary pathologist assessed the radiation changes and tumor response. We accrued 16 dogs, 8 dogs in each treatment arm. In the MBRT arm, 71% dogs achieved complete pathological remission. The radiation-related changes were all confined to the target region. Structural damage was not observed in the beam path outside of the target region. In contrast, none of the dogs in control group achieved remission and the treatment related damage was more extensive. Therapeutic superiority was observed with MBRT, including both tumor control and the normal structural preservation. The MBRT findings are suggestive of an immune related mechanism which is absent in standard treatment. These findings together with the widespread availability of clinical linear accelerators make MBRT a promising research topic to explore further treatment and clinical trial opportunities.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Structure of Shock Waves and Inelasticity in Shock-Compressed Cemented Tungsten Carbides

In this study, shock wave experiments are conducted on General Carbide cemented tungsten carbide with 3.7wt.% cobalt binder to determine its shock-induced compression behavior up to 100 GPa. The measured wave profiles indicate the cemented tungsten carbide to undergo elastic-plastic deformation during shock compression. A three-stage particle velocity profile is observed in the experiments -- an initial elastic-rise to the Hugoniot Elastic Limit (HEL), an elastic-plastic ramp indicating substantial post-yield hardening, and finally a rise to the peak shocked Hugoniot state. The results of the experiments are used to determine the HEL, the shock velocity (U s ) vs. particle velocity (u p ) Hugoniot relation, and the longitudinal stress (σ x ) vs. specific volume (V) curve for the samples. The HEL of the material was determined to lie between 4.41 GPa and 4.58 GPa. The U s - u p relation was determined to be U s = 4.97 + 1. 457u p for particle velocities greater than 0.75 km/s. The measured plastic shock velocities for particle velocities less than 0.7 km/s were found to be larger than those predicted using the linear U s - u p Hugoniot relationship, indicating the cemented WC samples to preserve substantial shear strength in the post-yield deformation region. No phase transformation was observed up to 100 GPa.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Celestial amplitudes as AdS-Witten diagrams

Both celestial and momentum space amplitudes in four dimensions are beset by divergences resulting from spacetime translation and sometimes scale invariance. In this paper we consider a (linearized) marginal deformation of the celestial CFT for Yang-Mills theory which preserves 2D conformal invariance but breaks both spacetime translation and scale invariance and involves a chirally coupled massive scalar. The resulting MHV celestial amplitudes are completely finite (apart from the usual soft and collinear divergences and isolated poles in the sum of the weights) and take the canonical CFT form. Moreover, we show they can be simply rewritten in terms of AdS 3 -Witten contact diagrams which evaluate to the well-known D-functions, thereby forging a direct connection between flat and AdS holography.

79 ASTRONOMY AND ASTROPHYSICS↗

Circularly Polarized Attosecond Pulses Enabled by an Azimuthal Phase and Polarization Grating

High-harmonic generation (HHG) is an extreme nonlinear optical process that can map the properties of an infrared driving laser beam onto short wavelength attosecond pulse trains. However, current techniques for generating circularly polarized high harmonics for probing magnetic materials and chiral systems have limitations: two-color collinear counter-rotating driving lasers result in a low cutoff photon energy, while single-color non-collinear counter-rotating schemes suffer from low conversion efficiency. In this work, we generate circularly polarized attosecond pulse trains by using a structured laser driver which has a rotating polarization and phase grating along the azimuthal coordinate. Furthermore, our experimental and numerical results demonstrate the production of left and right circularly polarized harmonics, which naturally separate upon propagation. Our approach uses a single laser color in a collinear geometry, that can be scaled for high efficiency. Simulations show this scheme can extend into the soft x-ray region when driven by mid-infrared driving lasers, while preserving the same high phase-matching cutoff photon energy as for linearly-polarized high harmonics.

attosecond science↗

Shock response of cobalt-based cemented tungsten carbides at pressures up to 100 GPa

In the present study, shock wave experiments are conducted on General Carbide cemented tungsten carbides with 3.7 wt.% and 6.0 wt.% cobalt binder to determine their shock compression response up to 100 GPa. A three-stage particle velocity profile is observed in the experiments -- an initial elastic-rise to the Hugoniot elastic limit (HEL), an elastic-plastic ramp indicating post-yield hardening, and finally a rise to the peak shocked Hugoniot state. The results of the experiments are used to determine the HEL, shock velocity (U s ) vs. particle velocity (u p ) Hugoniot relationship, and the longitudinal stress vs. specific volume Hugoniot relationship for the two samples. Here, the HEL for the WC 3.7wt% Co and WC 6.0wt% Co samples was determined to be ~4.45 GPa and 3.72 GPa, respectively. The U s – u p relation was determined to be U s = 4.97(0.006)+1. 446 (0.018) u p for WC 3.7 wt.% Co and U s = 4.93(0.006)+1. 454(0.017)u p for the WC 6.0 wt.% Co sample at peak particle velocities>0.75 km/s. For both WC grades, in the particle velocity regime less than ~0.75 km/s, the measured shock wave velocities were found to be larger than those predicted by the linear U s –u p Hugoniot relationship, indicating the two WC samples to preserve substantial strength in the post-yield deformation regime. Both WC grades show a catastrophic drop in shear stress carrying capacity when shocked to longitudinal stresses greater than ~70 GPa.

36 MATERIALS SCIENCE↗

Optimization of transformer ratio and beam loading in a plasma wakefield accelerator with a structure-exploiting algorithm

Plasma-based acceleration has emerged as a promising candidate as an accelerator technology for a future linear collider or a next-generation light source. We consider the plasma wakefield accelerator (PWFA) concept where a plasma wave wake is excited by a particle beam and a trailing beam surfs on the wake. For a linear collider, the energy transfer efficiency from the drive beam to the wake and from the wake to the trailing beam must be large, while the emittance and energy spread of the trailing bunch must be preserved. One way to simultaneously achieve this when accelerating electrons is to use longitudinally shaped bunches and nonlinear wakes. In the linear regime, there is an analytical formalism to obtain the optimal shapes. In the nonlinear regime, however, the optimal shape of the driver to maximize the energy transfer efficiency cannot be precisely obtained because currently no theory describes the wake structure and excitation process for all degrees of nonlinearity. In addition, the ion channel radius is not well defined at the front of the wake where the plasma electrons are not fully blown out by the drive beam. We present results using a novel optimization method to effectively determine a current profile for the drive and trailing beam in PWFA that provides low energy spread, low emittance, and high acceleration efficiency. We parameterize the longitudinal beam current profile as a piecewise-linear function and define optimization objectives. For the trailing beam, the algorithm converges quickly to a nearly inverse trapezoidal trailing beam current profile similar to that predicted by the ultrarelativistic limit of the nonlinear wakefield theory. For the drive beam, the beam profile found by the optimization in the nonlinear regime that maximizes the transformer ratio also resembles that predicted by linear theory. Furthermore, the current profiles found from the optimization method provide higher transformer ratios compared with the linear ramp predicted by the relativistic limit of the nonlinear theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Geometry of Non-Diffusive Tracer Transport in Gridded Atmospheric Models

A first-order linear and numerically non-diffusive Eulerian transport algorithm, minVAR, was recently developed for preservation of correlations between interrelated tracers during advective transport. The present study extends this work by: (1) providing further investigation of several interesting geometric constructions found in contours of constant minVAR, short for minimum variance, through extension to three coordinate dimensions. These contours capture point-by-point representations of thousands of individual atmospheric aerosol and/or cloud particles as they evolve and are rendered on Eulerian grids at a level of sub-grid resolution limited only by numerical precision; and (2) exploration of geometric similarities between the Arakawa C-grid, used to obtain interpolated values of the wind field at grid scale and minVAR. In particular, we consider interpolation of the u and v horizontal components of wind velocity from grid to sub-grid scales. The last results are motivated by recent applications of the Weather Research and Forecasting (WRF) model applied in the coastal Houston region, where the recent TRacking Aerosol Convection Interactions ExpeRiment (TRACER) field campaign was organized. A unique and fully consistent mapping is obtained between particles moving along meteorological wind trajectories and the non-diffusive, non-dispersive representation of such trajectories on an Eulerian grid.

54 ENVIRONMENTAL SCIENCES↗

An implicit-explicit time splitting strategy for the far SOL plasma fluid model with DG-FEM discretization

We consider a far scrape-off layer (SOL) plasma fluid model of ions that is governed by a Braginskiitype model: a one-dimensional, nonlinear system of advection-diffusion equations coupled with a diffusion equation for neutral particles. Our motivation for studying this system arises from the coupling between the edge plasma and radio-frequency (RF) heating, where solving a far SOL plasma fluid model provides critical insights into edge plasma dynamics. Numerical simulations of plasma fluid models require advanced computational techniques to achieve both efficiency and accuracy, especially when resolving the boundary layer in magnetically confined plasmas. In this work, we propose an implicit-explicit time operator splitting strategy that allows for an efficient solution algorithm, where the diffusive terms are treated semi-implicitly requiring only a linear solve, while the advection part is handled explicitly using a strong-stability-preserving Runge-Kutta (SSP-RK3) scheme. This leads to a fully decoupled system in which the diffusion and advection sub-problems can be solved separately, simplifying the overall solution procedure and allowing for efficient parallelization, which is particularly relevant for exploring the impact of RF heating on the SOL plasma. The main challenge of the discretization is due to the strong coupling between diffusion and advection, particularly through the boundary conditions. This makes implementation of such a scheme in an accurate and stable manner nontrivial. We discuss in detail how to split the equations and manage boundary conditions to maintain stability and well-posedness for each subsystem. We also describe a spatial discretization approach, based on the discontinuous Galerkin finite element method (DG-FEM) and present numerical results for a one-dimensional system.

Burkovska, Olena [ORNL] (ORCID:0000000163101130)↗

Parton distributions from boosted fields in the Coulomb gauge

We propose a new method to calculate parton distribution functions (PDFs) from lattice correlations of boosted quarks and gluons in the Coulomb gauge. Compared to the widely used gauge-invariant Wilson-line operators, these correlations greatly simplify the renormalization thanks to the absence of linear power divergence. Additionally, they enable access to larger off-axis momenta under preserved 3D rotational symmetry, as well as enhanced long-range precision that facilitates the Fourier transform. We verify the factorization formula that relates this new observable to the quark PDF at one-loop order in perturbation theory. Moreover, through a lattice calculation of the pion valence quark PDF, we demonstrate the aforementioned advantage and features of the Coulomb gauge correlation and show that it yields consistent results with the gauge-invariant method. This opens the door to a more efficient way to calculate parton physics on the lattice.

Baryons↗

A structure-preserving machine learning framework for accurate prediction of structural dynamics for systems with isolated nonlinearities

The nonlinearities present in structural systems are often found in isolated regions within the structure, such as those containing joints or interfaces. However, despite the localized nature of these nonlinearities their presence serves to couple together the modes of the underlying linear system and significantly complicate the development of appropriate reduced-order models; the localized nonlinearities have a global effect on the dynamics of the system. Further, in the presence of evolving structural health the nonlinearities can arise from accumulating damage, with dynamics distinct from those observed in the healthy state. The present work develops a data-driven formulation to identify and include the contributions of the isolated nonlinearities on the dynamics of the underlying linear structure. A novel coordinate separation is developed that decomposes those nonlinearities restricted to the isolated subdomain from the known linear system defined over the entire domain, and the influence of the isolated nonlinearities is reintroduced as an appropriately identified traction at the boundary of the isolated subdomain, referred to as the deviatoric force. In the region exterior to the nonlinear subdomain the response of the ideal linear system recovers that of the original nonlinear system. In this work, the deviatoric force component is predicted using a structure-preserving multilayer perceptron, based only on measured responses at the boundary of the isolated subdomain. Therefore introduction of the perceptron is able to bypass the direct numerical simulation of the nonlinearities within the isolated subdomain. This approach is illustrated through a simple structural system in which an interior region contains cubic nonlinearities and hysteretic damping. Once trained, the machine learning system is able to accurately predict the deviatoric force so that the ideal system recovers the response of the original system in the region outside the isolated nonlinear subdomain. Moreover, the data-driven approach is able to accurately predict the response when the system is subject to differing initial conditions and external excitation without the need for retraining, so that the proposed approach provides a robust description of the structural dynamics of the overall system.

Machine learning↗

Magnetic brightening and nanoscale imaging of spin-polarized helical edge modes in ZrTe 5

Efficient charge transport remains a fundamental challenge for nanoelectronic devices, as their performance is constrained by high dissipation and the impedance mismatch between high-frequency signal sources and nanoscale circuit interfaces. Although topologically protected helical edge modes offer a dissipationless and backscattering-resilient alternative, achieving nanoscale control over these modes requires direct visualization of their spatial electrodynamics, especially under high-frequency operation. Here, by utilizing cryogenic magneto-infrared scattering-type scanning near-field optical microscopy (cm-IR-sSNOM), we image spin-polarized helical edge channels in ZrTe 5 at the nanoscale, revealing magnetic-field-induced brightening as a high-frequency electrodynamic signature of robust topological edge modes. Operating at 1.8 K and under a magnetic field of up to 5 T, we observe the emergence of edge-state polarizability at infrared frequencies that is notably resilient to the magnetic gaps that typically quench d.c. and microwave edge transport. Our results reveal a topological ‘two-lane’ spatial reorganization in which an external magnetic field induces a spin-population imbalance between counterpropagating edge modes. Favoured helical branches are confined against physical boundaries to activate a net infrared near-field contrast. This electrodynamic response scales linearly with the number of atomic layers, which confirms that individual layers in ZrTe 5 preserve their discrete quantum spin Hall identities. These findings may be useful for developing topological spintronic devices, as the magnetic infrared tunability of helical edge channels provides a pathway for low-loss nanoscale interconnects and high-speed information processing.

36 MATERIALS SCIENCE↗

Quantum Langevin theory for two coupled phase-conjugated electromagnetic waves

We provide a general macroscopic phenomenological formula of quantum Langevin equations for two coupled phase-conjugated electromagnetic fields with linear loss (gain) and complex nonlinear coupling coefficient. The macroscopic phenomenological formula is obtained from the coupling matrix to preserve the field commutation relations and correlations, which does not require knowing the microscopic details of light-matter interaction and internal atomic structures. To validate this phenomenological formula, we take spontaneous four-wave mixing in a double-Λ four-level atomic system as an example to numerically confirm that our macroscopic phenomenological result is consistent with that obtained from the microscopic Heisenberg-Langevin theory. We find that a complex-valued nonlinear coupling coefficient can lead to noises even without linear gain or loss. Lastly, we apply the quantum Langevin equations to study the effects of linear gain and loss, complex phase mismatching, as well as complex nonlinear coupling coefficient in entangled photon pair (biphoton) generation, particularly to their temporal quantum correlations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

One-Body Properties and Their Perturbative Accuracy with Aufbau Suppressed Coupled Cluster Theory

In this work we derived and implemented the calculation of the one-body reduced density matrix for Aufbau suppressed coupled cluster theory, from which excited state natural orbitals and one-body properties, like atomic populations and dipole moments, are obtained. We utilized the natural orbitals to refine the ASCC solution for simple valence and Rydberg systems, exploring the process of repeatedly solving the ASCC equations in successive natural orbital bases to achieve independence from the starting molecular orbitals. For dipole moments in small molecules where high-level comparison data is available, we find that the accuracy of ASCC essentially matches that of linear response and equation-of-motion coupled cluster as long as care is taken to preserve the response's perturbative completeness.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗