Search NASA⌕ Search

SEARCH · Search NASA

Results for “quadratic”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

The Generality of the GUGA MRCI Approach in COLUMBUS for Treating Complex Quantum Chemistry

The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g., investigations of π-conjugated biradicaloid compounds, calculations of multitudes of excited states, development of diabatization procedures, and furnishing the electronic structure information for on-the-fly surface nonadiabatic dynamics. With fully variational uncontracted spin-orbit MRCI, COLUMBUS provides a unique possibility of performing high-level calculations on compounds containing heavy atoms up to lanthanides and actinides. Crucial for carrying out all of these calculations effectively is the availability of an efficient parallel code for the CI step. Configuration spaces of several billion in size now can be treated quite routinely on standard parallel computer clusters. Emerging developments in COLUMBUS, including the all configuration mean energy multiconfiguration self-consistent field method and the graphically contracted function method, promise to allow practically unlimited configuration space dimensions. Spin density based on the GUGA approach, analytic spin-orbit energy gradients, possibilities for local electron correlation MR calculations, development of general interfaces for nonadiabatic dynamics, and MRCI linear vibronic coupling models conclude this overview.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Percolation and nanosecond fluctuators in V 2 O 3 films within the metal–insulator transition

Vanadium sesquioxide (V 2 O 3 ) exhibits a metal–insulator transition (MIT) at 160 K between a low temperature, monoclinic, antiferromagnetic Mott insulator and a high temperature, rhombohedral, paramagnetic, metallic phase. In thin films, a percolative transition takes place over a finite temperature range of phase coexistence. We study the fluctuating dynamics of this percolative MIT by measuring voltage noise spectra at both low frequencies (up to 100 kHz) and radio frequencies (between 10 MHz and 1 GHz). Noise intensity quadratic in bias is observed in the MIT region, as expected for resistive fluctuations probed nonperturbatively by the current. The low frequency noise resembles flicker-type 1/f β noise, often taking on the form of Lorentzian noise dominated by a small number of fluctuators as the volume fraction of the insulating phase dominates. Radio frequency noise intensity also quadratic in the bias current allows the identification of resistance fluctuations with lifetimes below 1 ns, approaching timescales seen in non-equilibrium pump–probe studies of the transition. We find quantitative consistency with a model for fluctuations in the percolative fraction. The thermodynamics of the MIT suggests that dominant fluctuations are ones that alter small volumes affecting the connectivity of domain boundaries. This noise serves as a sensitive and nonperturbative probe for the dynamics of switching phenomena in this system.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Observations on the fourth-power scaling of high-pressure shock waves in solids

This paper is concerned with the bounds and approximations that are found in the original derivation of the fourth-power scaling of shock waves in solids [D. Grady, Appl. Phys. Lett. 38, 825 (1981)]. The analysis is focused on the framework of the derivation and is independent of constitutive assumption, such as visco-plastic behavior. Results include an upper bound for the shock pressure and a restriction on the range of power coefficient for materials having a power-type function for the shock velocity–particle velocity relation. Relaxation of this restriction is proposed based on the idea that the rise time along the Rayleigh line is strain dependent. The idea led to the application of the Onsager relation to strain rate calculation that in turn resulted in a quadratic function for the strain rate–shock pressure relation. The fourth-power relation is obtained by generalizing the Onsager relation through the introduction of a quadratic dissipative potential in analogy to Rayleigh's dissipative potential.

Horie, Yasuyuki (ORCID:0000000245207246)↗
DOE OSTI2023

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↗
DOE OSTI2021

Static response and Love numbers of Schwarzschild black holes

We derive the quadratic action for the physical degrees of freedom of massless spin-0, spin-1, and spin-2 perturbations on a Schwarzschild-(A)dS background in arbitrary dimensions. Here, we then use these results to compute the static response of asymptotically flat Schwarzschild black holes to external fields. Our analysis reproduces known facts about black hole Love numbers — in particular that they vanish for all types of perturbation in four spacetime dimensions — but also leads to new results. For instance, we find that neutral Schwarzschild black holes polarize in the presence of an electromagnetic background in any number of spacetime dimensions except four. Moreover, we calculate for the first time black hole Love numbers for vector-type gravitational perturbations in higher dimensions and find that they generically do not vanish. Along the way, we shed some light on an apparent discrepancy between previous results in the literature, and clarify some aspects of the matching between perturbative calculations of static response on a Schwarzschild background and the point-particle effective theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Reconstructing cosmic polarization rotation with ResUNet-CMB

Cosmic polarization rotation, which may result from parity-violating new physics or the presence of primordial magnetic fields, converts E-mode polarization of the cosmic microwave background (CMB) into B-mode polarization. Anisotropic cosmic polarization rotation leads to statistical anisotropy in CMB polarization and can be reconstructed with quadratic estimator techniques similar to those designed for gravitational lensing of the CMB. At the sensitivity of upcoming CMB surveys, lensing-induced B-mode polarization will act as a limiting factor in the search for anisotropic cosmic polarization rotation, meaning that an analysis which incorporates some form of delensing will be required to improve constraints on the effect with future surveys. In this study we extend the ResUNet-CMB convolutional neural network to reconstruct anisotropic cosmic polarization rotation in the presence of gravitational lensing and patchy reionization, and we show that the network simultaneously reconstructs all three effects with variance that is lower than that from the standard quadratic estimator nearly matching the performance of an iterative reconstruction method.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Synchronizing the consistency relation

We study the N-point function of the density contrast to quadratic order in the squeezed limit during the matter-dominated (MD) and radiation-dominated (RD) eras in synchronous gauge. Since synchronous gauge follows the free-fall frame of observers, the equivalence principle dictates that in the gradient approximation for the long-wavelength mode there is only a single, manifestly time-independent consistency relation for the N-point function. This simple form is dictated by the initial mapping between synchronous and local coordinates, unlike Newtonian gauge and its correspondingly separate dilation and Newtonian consistency relations. Dynamical effects only appear at quadratic order in the squeezed limit and are again characterized by a change in the local background, also known as the separate universe approach. We show that for the 3-point function the compatibility between these squeezed-limit relations and second-order perturbation theory requires both the initial and dynamical contributions to match, as they do in single-field inflation. This clarifies the role of evolution or late-time projection effects in establishing the consistency relation for observable bispectra, which is especially important for radiation acoustic oscillations and for establishing consistency below the matter-radiation equality scale in the MD era. Defining an appropriate angle and time average of these oscillations is also important for making separate universe predictions of spatially varying local observables during the RD era, which can be useful for a wider range of cosmological predictions beyond N-point functions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Dynamical chemistry: non-equilibrium effective actions for reactive fluids

Abstract We present two approaches for describing chemical reactions taking place in fluid phase. The first method mirrors the usual derivation of the hydrodynamic equations of motion by relating conserved—or to account for chemical reactions, non-conserved—currents to local-equilibrium parameters. The second method involves a higher-brow approach in which we attack the same problem from the perspective of non-equilibrium effective field theory (EFT). Non-equilibrium effective actions are defined using the in–in formalism on the Schwinger–Keldysh contour and are therefore capable of describing thermal fluctuations and dissipation as well as quantum effects. The non-equilibrium EFT approach is especially powerful as all terms in the action are fully specified by the symmetries of the system; in particular the second law of thermodynamics does not need to be included by hand, but is instead derived from the action itself. We find that the equations of motion generated by both methods agree, but the EFT approach yields certain advantages. To demonstrate some of these advantages we construct a quadratic action that is valid to very small distance scales—much smaller than the scales at which ordinary hydrodynamic theories break down. Such an action captures the full thermodynamic and quantum behavior of reactions and diffusion at quadratic order. Finally, taking the low-frequency and low-wavenumber limit, we reproduce the linearized version of the well-known reaction–diffusion equations as a final coherence check.

Mechanics↗
DOE OSTI2024

Analysis and Benchmarking of feature reduction for classification under computational constraints

Abstract Machine learning is most often expensive in terms of computational and memory costs due to training with large volumes of data. Current computational limitations of many computing systems motivate us to investigate practical approaches, such as feature selection and reduction, to reduce the time and memory costs while not sacrificing the accuracy of classification algorithms. In this work, we carefully review, analyze, and identify the feature reduction methods that have low costs/overheads in terms of time and memory. Then, we evaluate the identified reduction methods in terms of their impact on the accuracy, precision, time, and memory costs of traditional classification algorithms. Specifically, we focus on the least resource intensive feature reduction methods that are available in Scikit-Learn library. Since our goal is to identify the best performing low-cost reduction methods, we do not consider complex expensive reduction algorithms in this study. In our evaluation, we find that at quadratic-scale feature reduction, the classification algorithms achieve the best trade-off among competitive performance metrics. Results show that the overall training times are reduced 61%, the model sizes are reduced 6×, and accuracy scores increase 25% compared to the baselines on average with quadratic scale reduction.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Sparse Cholesky factorization for solving nonlinear PDEs via Gaussian processes

In recent years, there has been widespread adoption of machine learning-based approaches to automate the solving of partial differential equations (PDEs). Among these approaches, Gaussian processes (GPs) and kernel methods have garnered considerable interest due to their flexibility, robust theoretical guarantees, and close ties to traditional methods. They can transform the solving of general nonlinear PDEs into solving quadratic optimization problems with nonlinear, PDE-induced constraints. However, the complexity bottleneck lies in computing with dense kernel matrices obtained from pointwise evaluations of the covariance kernel, and its partial derivatives, a result of the PDE constraint and for which fast algorithms are scarce. The primary goal of this paper is to provide a near-linear complexity algorithm for working with such kernel matrices. We present a sparse Cholesky factorization algorithm for these matrices based on the near-sparsity of the Cholesky factor under a novel ordering of pointwise and derivative measurements. The near-sparsity is rigorously justified by directly connecting the factor to GP regression and exponential decay of basis functions in numerical homogenization. We then employ the Vecchia approximation of GPs, which is optimal in the Kullback-Leibler divergence, to compute the approximate factor. This enables us to compute ϵ-approximate inverse Cholesky factors of the kernel matrices with complexity O(N log d (N/ϵ)) in space and O(N log 2d (N/ϵ)) in time. We integrate sparse Cholesky factorizations into optimization algorithms to obtain fast solvers of the nonlinear PDE. We numerically illustrate our algorithm’s near-linear space/time complexity for a broad class of nonlinear PDEs such as the nonlinear elliptic, Burgers, and Monge-Ampère equations. In summary, we provide a fast, scalable, and accurate method for solving general PDEs with GPs and kernel methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Bias-hardened estimators of patchy screening profiles

Detecting anisotropic screening of the cosmic microwave background (CMB) holds the promise of revealing the distribution of gas in the Universe, characterizing the complex processes of galaxy formation and feedback, and studying the epoch of reionization. Estimators for inhomogeneous screening, including some recently proposed small-scale (stacked) estimators, are quadratic or higher order in the CMB temperature or polarization fields and are therefore subject to contamination from CMB lensing. We review the origin of this lensing bias and show that, when stacking on unWISE galaxies, the expected lensing bias dominates the signal if left unmitigated. Hardening techniques that null the lensing bias have been proposed for standard quadratic estimators, whereas only approximate methods have been proposed for stacked estimators. In conclusion, we review these techniques and apply the former to stacked estimators, presenting several strategies (including the optimal strategy) to null lensing contamination when stacking on any large-scale structure tracer.

Cosmic microwave background↗
DOE OSTI2022

Entanglement perspective on the quantum approximate optimization algorithm

Many quantum algorithms seek to output a specific bitstring solving the problem of interest—or a few if the solution is degenerate. It is the case for the quantum approximate optimization algorithm (QAOA) in the limit of large circuit depth, which aims to solve quadratic unconstrained binary optimization problems. Hence, the expected final state for these algorithms is either a product state or a low-entangled superposition involving a few bitstrings. What happens in between the initial N -qubit product state | 0 〉 ⊗ N and the final one regarding entanglement? Here, we consider the QAOA algorithm for solving the paradigmatic MaxCut problem on different types of graphs. We study the entanglement growth and spread resulting from randomized and optimized QAOA circuits and find that there is a volume-law entanglement barrier between the initial and final states. We also investigate the entanglement spectrum in connection with random matrix theory. In addition, we compare the entanglement production with a quantum annealing protocol aiming to solve the same MaxCut problems. Finally, we discuss the implications of our results for the simulation of QAOA circuits with tensor network-based methods relying on low-entanglement for efficiency, such as matrix product states.

Dupont, Maxime↗
DOE OSTI2021

Finite and infinite matrix product states for Gutzwiller projected mean-field wave functions

Matrix product states (MPS) and “dressed” ground states of quadratic mean fields (e.g., Gutzwiller projected Slater determinants) are both important classes of variational wave functions. This latter class has played important roles in understanding superconductivity and quantum spin liquids. We present a method to obtain both the finite and infinite MPS (iMPS) representation of the ground state of an arbitrary fermionic quadratic mean-field Hamiltonian (which in the simplest case is a Slater determinant and in the most general case is a Pfaffian). We also show how to represent products of such states (e.g., determinants times Pfaffians). From this representation one can project to single occupancy and evaluate the entanglement spectra after Gutzwiller projection. We then obtain the MPS and iMPS representation of Gutzwiller projected mean-field states that arise from the variational slave-fermion approach to the S = 1 bilinear-biquadratic quantum spin chain. To accomplish this, we develop an approach to orthogonalize degenerate iMPS to find all the states in the degenerate ground-state manifold. We find the energies of the MPS and iMPS states match the variational energies closely, indicating the method is accurate and there is minimal loss due to truncation error. We then present an exploration of the entanglement spectra of projected slave-fermion states, exploring their qualitative features and finding good qualitative agreement with the respective exact ground-state spectra found from density matrix renormalization group.

1-dimensional spin chains↗
DOE OSTI2023

Nonlinear photoconductivities and quantum geometry of chiral multifold fermions

Here, chiral multifold fermions are quasiparticles that appear only in chiral crystals such as transition metal silicides in the cubic B20 structure (i.e., the CoSi family), and they may show exotic physical properties. Here we study the injection and shift photoconductivities and also the related geometrical quantities for several types of chiral multifold fermions, including spin-1/2 as well as pseudospin-1 and -3/2 fermions, dubbed as Kramers Weyl, triple-point, and Rarita-Schwinger-Weyl (RSW) fermions, respectively. We utilize the minimal symmorphic model to describe the triple-point fermions (TPF). We also consider the more realistic model Hamiltonian for the CoSi family including both linear and quadratic terms. We find that injection currents due to circularly polarized light are quantized as a result of the Chern numbers carried by the multifold fermions within the linear models. Surprisingly, we discover that in the TPF model, the linear shift conductivities, responsible for the shift current generation by linearly polarized light, are proportional to the pseudo spin-orbit coupling and independent of photon frequency. In contrast, for the RSW and Kramer Weyl fermions, the linear shift conductivity is linearly proportional to photon frequency. The numerical results agree with the power-counting analysis for quadratic Hamiltonians. The frequency independence of the linear shift conductivity could be attributed to the strong resonant symplectic Christoffel symbols of the flat bands. Moreover, the calculated symplectic Christoffel symbols show significant peaks at the nodes, suggesting that the shift currents are due to the strong geometrical response near the topological nodes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Electronic density of states of a U ( 1 ) quantum spin liquid with spinon Fermi surface. I. Orbital magnetic field effects

A quantum spin liquid with a spinon Fermi surface is an exotic insulator that hosts neutral Fermi surfaces inside the insulating gap. In an external magnetic field, it has been pointed out that the neutral Fermi surfaces are Landau quantized to form Landau levels due to the coupling to the induced emergent gauge magnetic field. In this paper, we calculate the electronic density of states (as observed in tunneling experiments) of the quantum spin liquid in an orbital magnetic field. Here, we find that the Landau levels from the neutral Fermi surfaces give rise to a set of steps emerging at the upper and lower Hubbard band edges. Each of the Hubbard band edge steps further develop into a band edge resonance peak when a weak gauge binding arising from the gauge field fluctuations is taken into account. Importantly, each Hubbard band edge step and its resulting resonance peak in the weak gauge binding are found to have a corresponding Landau level from the neutral Fermi surfaces, so the Hubbard band edge steps and the band edge resonance peaks provide signatures to the unique feature that characterizes the Landau quantization of the in-gap neutral Fermi surfaces in the spin liquid. We further consider the strong gauge-binding regime where the band edge resonance peaks move into the Mott gap and develop into true in-gap bound states. In the strong gauge binding regime, we solve the Landau-level spectrum of the in-gap bound states in an orbital magnetic field. For the in-gap bound state with a Mexican-hat-like band dispersion, we find that the envelop energy to have a state excited from the bound state Landau levels decreases quadratically with the magnetic field. The quadratic decrease behavior of the envelop energy is consistent with the intuition that applying magnetic field localizes the states and energetically promotes the in-gap bound states formation. Finally, we discuss the connection of our results to the electronic density of states spectra measured in the layered 1T-TaS 2 . We point out that a quantum spin liquid with a quasibound state in the upper Hubbard band can give the density of states spectra similar to the one measured in the experiment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

High-precision mass measurements of the isomeric and ground states of 44 V: Improving constraints on the isobaric multiplet mass equation parameters of the A=44, 0 + quintet

Background: The quadratic isobaric multiplet mass equation (IMME) has been very successful at predicting the masses of isobaric analog states in the same multiplet, while its coefficients are known to follow specific trends as functions of mass number. The Atomic Mass Evaluation 2016 [Chin. Phys. C 41, 030003 (2017)] 44 V mass value results in an anomalous negative c coefficient for the IMME quadratic term; a consequence of large uncertainty and an unresolved isomeric state. The b and c coefficients can provide useful constraints for construction of the isospin-nonconserving Hamiltonians for the pf shell. In addition, the excitation energy of the 0 + , T = 2 level in 44 V is currently unknown. This state can be used to constrain the mass of the more exotic 44 Cr. Purpose: The aim of the experimental campaign was to perform high-precision mass measurements to resolve the difference between 44 V isomeric and ground states, to test the IMME using the new ground state mass value and to provide necessary ingredients for the future identification of the 0 + , T = 2 state in 44 V. Method: High-precision Penning trap mass spectrometry was performed at LEBIT, located at the National Superconducting Cyclotron Laboratory, to measure the cyclotron frequency ratios of [ 44g,m VO] + versus [ 32 SCO] + , a well-known reference mass, to extract both the isomeric and ground state masses of 44 V. Results: The mass excess of the ground and isomeric states in 44 V were measured to be -23 804.9(80) keV/c 2 and -23 537.0(55) keV/c 2 , respectively. This yielded a new proton separation energy of S p = 1 773(10) keV. Conclusion: The new values of the ground state and isomeric state masses of 44 V have been used to deduce the IMME b and c coefficients of the lowest 2 + and 6 + triplets in A = 44. The 2 + c coefficient is now verified with the IMME trend for lowest multiplets and is in good agreement with the shell-model predictions using charge-dependent Hamiltonians. The mirror energy differences were determined between 44 V and 44 Sc, in line with isospin-symmetry for this multiplet. Finally, the new value of the proton separation energy determined, to an uncertainty of 10 keV, will be important for the determination of the 0 + , T = 2 state in 44 V and, consequently, for prediction of the mass excess of 44 Cr.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

First Penning trap mass measurement of 36 Ca

Background: Isobaric quintets provide the best test of the isobaric multiplet mass equation (IMME) and can uniquely identify higher order corrections suggestive of isospin symmetry breaking effects in the nuclear Hamilto nian. The Generalized IMME (GIMME) is a novel microscopic interaction theory that predicts an extension to the quadratic form of the IMME. Only the A = 20, 32 T = 2 quintets have the exotic T z = –2 member ground state mass determined to high-precision by Penning trap mass spectrometry. Purpose: To establish A = 36 as the third high-precision T = 2 isobaric quintet with the T z = –2 member ground state mass measured by Penning trap mass spectrometry and provide the first test of the predictive power of the GIMME. Method: Here, a radioactive beam of neutron-deficient 36 Ca was produced by projectile fragmentation at the National Superconducting Cyclotron Laboratory. The beam was thermalized and the mass of 36 Ca + and 36 Ca 2+ measured by the Time of Flight - Ion Cyclotron Resonance method in the LEBIT 9.4 T Penning trap. Results: We measure the mass excess of 36 Ca to be ME= –6483.6(56) keV, an improvement in precision by a factor of 6 over the literature value. The new datum is considered together with evaluated nuclear data on the A = 36, T = 2 quintet. We find agreement with the quadratic form of the IMME given by isospin symmetry, but only coarse qualitative agreement with predictions of the GIMME. Conclusion: A total of three isobaric quintets have their most exotic members measured by Penning trap mass spectrometry. The GIMME predictions in the T = 2 quintet appear to break down for A = 32 and greater.

20 ≤ A ≤ 38↗
DOE OSTI2021

Constraints on the nuclear symmetry energy from asymmetric-matter calculations with chiral 𝑁⁢𝑁 and 3⁢𝑁 interactions

The nuclear symmetry energy is a key quantity in nuclear (astro)physics. It describes the isospin dependence of the nuclear equation of state, which is commonly assumed to be almost quadratic. Here, in this work, we confront this standard quadratic expansion of the equation of state with explicit asymmetric nuclear-matter calculations based on a set of commonly used Hamiltonians including two- and three-nucleon forces derived from chiral effective-field theory. We study, in particular, the importance of nonquadratic contributions to the symmetry energy, including the nonanalytic logarithmic term introduced by Kaiser [Phys. Rev. C 91, 065201 (2015)]. Our results suggest that the nonquadratic contribution to the symmetry energy can be systematically determined from the various Hamiltonians employed, and we obtain 0.74$^{+0.11}_{−0.08}$ MeV (or −1.02$^{+0.11}_{−0.08}$ MeV for the potential term with the effective-mass contribution) at nuclear saturation density, while the logarithmic contribution to the symmetry energy is relatively small and model-dependent. We also employ the meta-model approach to study the impact of the higher-order contributions on the neutron-star crust-core transition density, and find a 5% correction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗