2026-09-28 | | Total: 118
Skyrmions experience a sideways gyroscopic force when they move. Bound molecules can cancel this force; we call them gyroscopically neutral. We show that the quantum states available on each site of a magnet select the composition needed for cancellation. In two models of $SU(3)$ magnetism, this composition changes from a pair to a trimer. Both models assign exactly the same energy to every configuration of their common classical order parameter, isolating the effect of the local quantum states. Neutral molecules can also unwind when additional local states are allowed. Mean field calculations nevertheless find both compositions bound and locally stable, with the neutral trimer surviving sampled interaction changes. Comparing their motion reveals the changed cancellation rule, and prepared neutral molecules translate over short distances with little deformation.
Isoelectronic compounds are generally expected to exhibit related electronic, structural, and magnetic behavior. Chromium monopnictides CrX (X= Sb, As, P), however, display a striking evolution from high-temperature altermagnetism in NiAs-type CrSb, through double-helical antiferromagnetism in MnP-type CrAs, to the absence of resolved long-range magnetic order in MnP-type CrP. Here, combining first-principles calculations in the non-collinear Curie-paramagnetic state, structural analysis, magnetic phase-diagram calculations, and exchange-parameter extraction, we establish a unified microscopic framework for this evolution. All three compounds share a formal high-spin $d^5$ configuration, nominally Cr$^{1+}$, coupled to itinerant states of predominantly pnictogen character. From Sb to P, chemical pressure enhances ligand itinerancy and drives the structural evolution from the high-symmetry NiAs lattice to increasingly distorted MnP-type networks. The accompanying reconstruction of competing Cr-Cr exchanges favors A-type antiferromagnetism in CrSb and double-helical order in CrAs, while placing CrP near a frustrated regime with closely competing magnetic tendencies. Our results establish chemical-pressure-driven exchange reconstruction as a general mechanism through which a common spin-fermion electronic structure can generate contrasting magnetic phases in strongly correlated materials.
Geometrically frustrated magnetism conventionally arises from pre-existing magnetic moments on lattices whose geometry prevents their interactions from being simultaneously satisfied, giving rise to highly degenerate states and rich collective behavior. Creating such frustration in an intrinsically non-magnetic material presents a fundamentally different challenge, requiring both the magnetism and the competing interactions to emerge from correlated electrons. Here we show that this can be realized in graphene through lithographically programmable strain engineering. Patterning strain and the associated pseudo-magnetic field (PMF) into superlattices creates a correlated electronic system with flat bands and strong interactions. Transport measurements reveal interaction-driven insulating behavior, anisotropic magnetic hysteresis, and slow relaxation dynamics reminiscent of spin freezing, qualitatively captured by Monte Carlo simulations of competing magnetic moments on the PMF-defined ruby superlattice. Cryogenic magnetic force microscopy further reveals magnetic textures associated with the PMF landscape. These observations demonstrate frustrated magnetism emerging from correlated electrons in otherwise non-magnetic graphene. Our results establish lithographic strain engineering as a general and scalable route to flat bands and correlated states in van der Waals materials, providing a versatile platform for programmable quantum matter.
Strain engineering of magnetic properties offers a promising route toward flexible spintronic applications. Here, we report an unconventional magnetoelastic response in Al rich Co Fe Al thin films on flexible substrates using strain dependent magneto optical Kerr effect microscopy and magnetometry measurements. While the films exhibit a conventional positive magnetostriction coefficient, consistent with standard in plane easy axis rotation under stress, their saturation magnetization increases under compressive strain and decreases under tensile strain. First principles calculations reveal that this unconventional response originates from a strain induced competition between exchange splitting and crystal field effects in an inverse Heusler like local environment created by Al enrichment. This leads to a highly sensitive, sublattice dependent magnetic state, consistent with a strain induced reconfiguration of Co and Fe moments. Our results demonstrate that local compositional tuning can fundamentally alter magnetoelastic behavior, establishing strain controlled sublattice compensation as a route toward programmable magnetic functionality in flexible spintronic systems.
The structural and optical properties of the doubly yttrium-doped boron cluster B$_{18}$Y$_2$ have been systematically investigated using density functional theory calculations. The lowest-energy structure was identified through extensive basin-hopping searches accelerated by a pre-trained MACE machine-learning potential and subsequently refined and validated at the DFT level. The resulting global-minimum structure adopts a double-ring geometry, consisting of two fused B$_9$ rings stabilized by yttrium atoms positioned above and below the boron framework. Vibrational frequency calculations confirm the dynamical stability of the optimized structure, while the calculated infrared and UV--Vis spectra provide characteristic signatures of the Y--B interactions and the electronic structure of the boron framework. These results demonstrate how yttrium doping can stabilize unusual double-ring boron architectures and highlight the effectiveness of machine-learning-assisted basin-hopping searches for exploring the complex potential-energy landscapes of doped boron clusters.
Frequently, nonextensive systems are described through $q$-exponential functions, where the $q$-parameter determines how far a system is from equilibrium. In this work, we study systems constrained to $q-1\ll1$ to describe systems near thermal equilibrium, and derive the probability density function and temperature fluctuations. To this end, we derive the series expansion of the $q$-exponential function, yielding a perturbative series that is a superposition of gamma distributions. On the other hand, the temperature fluctuations are expressed as a superposition of Dirac delta functions and their derivatives, meaning that the bulk of the system remains at a constant temperature, while point-temperature variations emerge, marking the system's departure from thermal equilibrium. The statistics of the temperature fluctuations exhibit a transition from a symmetric to a nonsymmetric distribution as the nonextensive parameter increases. These results provide evidence of the association between the emergence of nonextensivity and the symmetry breaking of temperature fluctuations.
Sintering is widely used to manufacture ceramics, but coupled densification and grain growth, material-dependent kinetics, and sparse measurements complicate predictive modeling and process design. We present Sinter-PiNDiff, a retrainable physics-integrated neural differentiable framework for predicting density and grain-size evolution. Two neural networks learn densification and grain-growth coefficients within coupled rate equations, while a smooth saturation factor attenuates densification near theoretical density. The same governing structure, network architecture, and training procedure were fitted independently to published data for MgO, Al-doped ZnO, and CaO-doped ThO2. Tests at held-out temperatures and compositions yielded the lowest mean error in all twelve material-metric comparisons against multilayer perceptron and residual network baselines. For MgO, Al-doped ZnO, and CaO-doped ThO2, respectively, density normalized root-mean-square errors were 14.6%, 10.8%, and 14.4%, and grain-size errors using the same metric were 8.6%, 12.1%, and 19.3%. Removing evolving density from both neural-network inputs increased density and grain-size trajectory errors in all three systems and ten of twelve aggregate errors, supporting density-dependent kinetic feedback. Deep ensembles estimated model disagreement, but empirical coverage showed that the uncertainty bands were not calibrated and did not capture all model-data discrepancies. These results establish Sinter-PiNDiff as a retrainable framework for sparse-data prediction and uncertainty-informed selection of sintering conditions.
Interaction between Rydberg excitons in core-shell Cu$_2$O nanostructures and surface plasmons is investigated. An enhancement of a quadrupole transitions to the Rydberg exciton states in copper oxide plasmonic nanaostructures is examined. A possibility of a modification of plasmonic enhancment (i.e., an average amplification factor) of electric field distributions in the vicinity of nanostructures, which depends on its shape and geometrical size, is also discussed.
The influence of crystallographic orientation on the early stages of the oxidation of pure Zr single crystals was studied by X ray photoelectron spectroscopy (XPS). Two samples, cut from the same a-Zr single crystal, were oxidized simultaneously at room temperature and O2 pressures of 1 X 10-8, 1 X 10-7 and 1 X 10-6 torr: one with its surface normal parallel to the c-axis (Z1, basal (0001) plane) and the other with its normal 12deg from that of a prismatic plane (Z2, near-prismatic orientation). At all three pressures the oxidation kinetics of both samples followed a three stage, logarithmic type behaviour, but the near prismatic sample (Z2) incorporated oxygen faster and to a greater extent than the basal sample (Z1) in every case. Deconvolution of the Zr 3d core level spectra resolved, in addition to metallic Zr and ZrO2, two sub-stoichiometric Zr O compounds, denoted (ZrO)a and (ZrO)b, with binding-energy shifts of 1.3 and 2.3 eV with respect to metallic Zr, consistent with sub oxides previously reported for polycrystalline Zr. Angle resolved XPS showed that ZrO2 is the outermost compound in both orientations, while the two sub-oxides are distributed nearly homogeneously through the film. Oxide thicknesses, calculated from the attenuation of the metallic Zr 3d signal, ranged from 10 to 13 A for Z1 and from 14 to 19 A for Z2 depending on the oxidation pressure, in close agreement with earlier measurements on cold worked polycrystalline Zr, whose surface texture is dominated by prismatic-oriented grains.
Ferroelectric switching in heterostructures couples composition, layer topology, temperature, and interfacial boundary conditions. ReaxFF molecular dynamics isolates these variables in ZnO/ Zn1- xMgxO/ZnO and Zn1-xMgxO/ZnO/ Zn1-xMgxO stacks. Within pristine, initially single-domain models, coupling to switchable Zn1-xMgxO reduces the applied field required to reverse ZnO by up to fivefold. Temperature generally lowers the coercive field, whereas Mg concentration produces a nonmonotonic response. At equal ZnO and Zn1-xMgxO (ZMO) proportions, structures with ZnO at the center switch at lower fields than those with ZMO at the center at all four temperatures examined, demonstrating a topology-dependent response. Layer-resolved trajectories reveal topology-dependent switching sequences with direction-dependent redistribution of normal stress near the heterointerfaces, consistent with a stress-assisted cooperative pathway. Limiting MgO-containing structures exhibit sequential multilevel switching or low-polarity trapping, depending on thickness and temperature. Fixed-charge atomistic simulations complement previous continuum descriptions by resolving structural, energetic, and local stress evolution under a common applied field. TEM and STEM-EDS observations provide experimental structural context for the modeled architectures. Together, the results establish layer topology and interfacial mechanical confinement as design variables for wurtzite ferroelectric heterostructures.
Non-Hermitian skin effect, the accumulation of bulk eigenstates at edges under open boundary conditions, has been observed across classical, quantum, and synthetic platforms. Whether it can be generalized to the genuine time dimension remains open. We propose a spatiotemporal NHSE in the Floquet phase space via Floquet engineering of a ring-trapped ultracold atomic gas. Placing spatial and temporal degrees of freedom on an equal footing, we show that a single parent skin effect projects onto either domain. The spatial projection yields an NHSE with arbitrarily tunable skin locations, eliminating the need for physical open boundaries. The temporal projection gives time-domain nonreciprocity and periodic temporal funneling, consistent with causality. This establishes a unified spatiotemporal paradigm for NHSE and opens a route to non-Hermitian topological physics in time.
Finite-range hopping constrains both the number of Dirac points and the masses that gap them. We determine the maximal absolute Chern numbers of an isotropic finite-range extension of the Qi--Wu--Zhang model on a square lattice through thirteenth-neighbor hopping. Algebraic zero counting and mass-weighted sum rules yield global bounds, and we construct explicit Hamiltonians that saturate them. We find $|C|_{\max}=36,40,51$ for eleventh-, twelfth-, and thirteenth-neighbor hopping. Equal numbers of Dirac points can produce different maxima because finite-range mass harmonics distinguish opposite vortex charges differently. Twice the squared hopping radius bounds every maximum through twelfth-neighbor hopping and is saturated at the eleventh and twelfth, but the thirteenth-neighbor maximum exceeds it. Thus finite-range hopping constrains the maximal Chern number in two complementary ways: through the number of Dirac vortices it can generate and through the mass-sign patterns it can assign to them. This geometric viewpoint provides a design principle for high-Chern bands and a route to analogous bounds in moiré and other lattice systems.
Linear dynamics of the small amplitude perturbations in the antiferromagnetic multiferroics is analytically considered. Equilibrium state in the considered regime is assumed to be cycloidal structure both for the antiferromagnetic vector and the nonzero magnetization vector (presence of weak ferromagnetism in the system is included). The electric polarization of spin origin related to the "noncollinear" spins is considered and its influence on the spin dynamics is included. The dynamical electric susceptibility tensor is obtained to demonstrate response of the multiferroics on the alternating electric field. Existence of two polarization projections is found, response of each of them on two projections of the electric field perpendicular to the cycloid direction is presented. Different roles of the in plane and out of plane electric perturbations is demonstrated analytically. Essential role of the low frequency gapless spin wave in the dielectric response is described. Contributions of the spiral structure and the weak ferromagnetism in the dielectric response are specified.
HfO$_2$-based ferroelectric tunnel junctions (FTJs) are promising candidates for scalable non-volatile memory, but simultaneously achieving a high tunneling electro-resistance ratio (TER) and a low resistance-area (RA) product remains challenging. To address this challenge, this work introduces non-uniform quantum well (QW) and barrier thickness engineering in HfO$_2$-based multi-QW FTJs using a self-consistent Preisach-based ferroelectric (FE) model integrated with the coherent and inelastic non-equilibrium Green's function (NEGF) formalism. The non-uniform well and barrier configuration produces a wide range of FTJ design landscapes due to closely spaced, broad resonant states in the low resistance state (LRS) and a larger separation in the high resistance state (HRS), resulting in strong polarization-dependent resonant transmission with TER reaching the order of $\mathbf{1\times10^{8}\%}$ and an LRS RA product as low as $\mathbf{1~Ω\cdot\mathrm{cm}^{2}}$ at read bias and in the presence of scattering. By incorporating self-consistent elastic scattering into the NEGF framework, we also show that elastic scattering can positively influence the TER and RA performance of non-uniform FTJs by progressively increasing the overlap of the closely spaced resonances. Overall, the results establish non-uniform QW and barrier thickness as an effective and robust design parameter for controlling resonant-state alignment and achieving a favorable combination of high TER and low RA in HfO$_2$-based multi-QW FTJs.
We use density functional theory (DFT) to demonstrate rearrangement of charge disproportionation as a new mechanism for screening polar discontinuities at heterointerfaces. We focus on SmNiO$_3$, a III-III perovskite characterized by an insulating ground state that originates from disproportionation of the electrons at the Ni sites. When interfaced with ferroelectric, II-IV perovskite BaTiO$_3$, this system can display a polar discontinuity whose magnitude is determined by the nature of the interface (NiO$_2$-BaO or SmO-TiO$_2$) and the orientation of the BaTiO$_3$ spontaneous polarization. In the orientation maximizing the interface charge, SmNiO$_3$ compensates by transferring 1 $e$/f.u. between the two interfaces and modifying the disproportionation sequence from Ni$^{4+}$/Ni$^{2+}$ to Ni$^{2+}$/Ni$^{2+}$ at the SmO-TiO$_2$ interface and to Ni$^{4+}$/Ni$^{4+}$ at the NiO$_2$-BaO interface. While this mechanism is potentially universal in all perovskites exhibiting charge disproportionation, our formal analysis of BiNiO$_3$/BaTiO$_3$ superlattices suggests that larger polar discontinuities may need more complex disproportionation patterns for screening.
Grain boundaries influence hydrogen transport in polycrystalline metals by acting as both trapping sites and interconnected diffusion pathways. Conventional interpretations of thermal desorption spectroscopy (TDS) and permeation experiments often treat traps as isolated defects, neglecting grain boundary connectivity and potentially misinterpreting experimental data. Here, we develop a quantitative phase-field model for hydrogen diffusion, grain boundary trapping, and grain boundary-assisted transport. The formulation preserves the physical grain boundary volume independently of the numerical interface thickness and accounts for anisotropic diffusion parallel and perpendicular to grain boundaries. Benchmark simulations verify quantitative behavior under interface upscaling. The model shows that grain boundary diffusion can significantly shift TDS peaks, making Kissinger-type analyses unreliable when grain boundary transport is significant. Effective diffusion coefficients are also shown to depend strongly on grain size, trapping free energy, temperature, and grain boundary mobility. These results highlight the coupled role of trapping thermodynamics, bulk diffusion, and grain-boundary transport.
The decay of a density wave in an interacting adlayer obeys an exact memory equation of Mori's form. Being a Fourier component of a conserved density, its correlation function is a characteristic function, and the kernel is built from a second one, that of a single jump. The local rate obeys a sum rule whose momentum dependence is one minus the characteristic function of a single jump, divided by the static structure factor. Detailed balance makes the decay a superposition of relaxation modes of non-negative weight whose moments are static equilibrium averages, so that the truncations of the associated continued fraction form a convergent hierarchy of closed forms, each of them a bound on the decay itself. The single rate a fit returns is shown to be the mean of the instantaneous rate over a time window, hence never above the local rate. At long wavelength, the two coincide, the relaxation spectrum separating there, one collective mode taking almost the whole weight of the density wave and the single exponential usually fitted to spin-echo data being the decay itself. What an experiment reaches at every wavevector is also the area under the decay, a correlation time. Together with the static structure factor and the single-jump geometry, it forms one combination which the theory identifies with the effective hop rate divided by the product of the local rate and the correlation time. That product never falls below unity and tends to unity at long wavelength whenever the rate at which an adparticle leaves a site is independent of the site it leaves for. In the hydrodynamic limit, the sum rule becomes a relation of Darken type in which the rate multiplied by the thermodynamic factor is the effective hop rate and not the mobility of a labelled adparticle.
In a recent paper Nahum argued on the basis of a renormalization-group analysis that, contrary to standard lore, ground states of 1D spin chains with short-range interactions can spontaneously break U(1) ``easy-plane'' spin rotation symmetry. True long-range order of $(S^x,S^y)$ arises at the phase transition between two quasi-long-range-ordered phases. Here we present detailed numerical results for an anisotropic spin-1 chain model. We find a magnetization exponent $β\approx0.29$, consistent with the $ε$-expansion prediction $β\approx0.28$, a divergence of the Luttinger parameter on approaching the transition from both sides with an exponent that we relate to Nahum's transition, and finite-size spectra consistent with a dynamical critical exponent $z\approx2$. Our results support the critical behavior predicted by Nahum's field theory.
P-type transparent conducting materials (TCMs) remain scarce, limiting progress in transparent electronics and tandem photovoltaics. Here, we perform a high-throughput computational search across CuI-based ternary X--Cu--I compounds using the minima hopping method combined with density functional theory. Filtering for thermodynamic stability, optical transparency, and light hole effective mass, we identify 58 candidate p-type TCMs, of which 49 are previously unreported. Twenty are p-type degenerate semiconductors with intrinsic hole carriers, including Cu$_2$FI$_2$ (calculated gap 3.75 eV, hole effective mass $0.262\,m_0$). We further use this dataset to revisit the chemical modulation of the valence band (CMVB) design rules, originally formulated for Cu-based delafossite oxides. Four generalizations emerge: linear Cu coordination is not required, as tetrahedral coordination (CN = 4) is compatible with p-type TCM character; the mechanism extends to halide-based systems with I $5p$ anion states; the Cu $d$ state energy center $ε_d$ is a quantitative descriptor for valence band dispersion; and Cu undercoordination in monovalent-anion hosts drives intrinsic p-type carrier generation through stoichiometric charge balance. These results establish a generalized computational framework for the design of p-type halide TCMs.
Polymers are widely used as functional insulating materials, but predicting their dielectric behavior is challenging because multiple mechanisms, acting across different spatial and temporal scales, contribute to their response. In polyethylene (PE), one of the most common polymers, electronic, atomic, and mesoscale dynamics must all be considered. In this work, semicrystalline models of polyethylene were developed and investigated using a multiscale approach. Quantum simulations were employed to describe electronic and vibrational contributions to the dielectric response, while classical molecular dynamics simulations captured the behavior of polymer chains at room temperature and frequencies down to the 10-100 MHz range. The study focuses on the impact of radio-oxidation defects on the dielectric properties of polyethylene. Quantum calculations reveal how electronic and vibrational contributions depend on the local atomic environment surrounding these defects. Comparison with molecular dynamics results highlights the influence of temperature on the static dielectric constant. Structural analyses further assess the effect of each defect type on PE crystallinity. Analysis of dipolar correlation functions shows that different defects interact with one another and with the polymer matrix in distinct ways, affecting permittivity and dielectric loss peaks. For all defect types, defect-defect interactions provide the largest contribution to dielectric response. While most oxidized groups exhibit positive coupling with the PE matrix, alcohol defects display a negative cross-correlation, partially offsetting their impact on the dielectric properties. Our results also show that ketone groups produce the largest dielectric loss between the defects studied.
We show that a linear damped elastic system can display a Mpemba-like reversal of vibrational relaxation without any nonlinear constitutive response or amplitude-dependent damping. The mechanism is purely modal. For a one-dimensional Kelvin--Voigt elastic medium, the decay rate of the $n$-th vibrational mode scales as $γ_n\propto n^2$. We formulate relaxation in terms of a positive modal envelope energy, thereby excluding crossings caused merely by oscillation phase. A one-parameter family of initial states is then constructed such that increasing the initial elastic excitation simultaneously reduces the projection onto the slow fundamental mode. In the invariant two-mode subspace, the relaxation-order crossing is obtained analytically and, for this family, occurs at a universal dimensionless time independent of the selected pair of initial states. In the limiting case where the slowest mode is absent, the dominant relaxation rate changes discontinuously, yielding a direct vibrational analogue of the strong Mpemba effect. We further show that the same relaxation-order reversal occurs in the ordinary mechanical vibrational energy, which decreases monotonically despite retaining the oscillatory dynamics of the underlying modes.
Recent advances in experiments using individually trapped atoms have enabled precise control over circular Rydberg electrons with exceptionally long lifetimes. We show that these giant and stable electron orbits can form toroidal lattice potentials for ultracold atoms with a period set by the electron's de Broglie wavelength. Unlike conventional static optical lattices, this electron lattice is formed via the electron-atom interaction, which mixes Rydberg circular states with opposite azimuthal phase winding into a standing electronic matter wave. The lattice phase is thereby intrinsically coupled to the atom position. For a pair of atoms, this results in ballistic tunneling motion in the two-atom spatial correlations along the ring, while the single particle dynamics is essentially free rotation. We simulate the dynamics for an experimentally realistic setting that exploits optical tweezers for individual atom control. Our results open a route toward incorporating long range atom-atom interactions mediated by a single electron and, ultimately, realizing small Bose and Fermi gases confined in these microscopic electronic atom traps.
Computing non-dilute diffusion coefficients from Onsager transport coefficients of a multicomponent alloy requires the thermodynamic factor, which measures the curvature of the free energy with respect to composition. As diffusion in metals is mediated by vacancies, the curvature that matters is not that of the vacancy-free alloy, which is straightforward to compute, but that of an alloy carrying a dilute concentration of vacancies. The resulting matrix is nearly singular, and composition fluctuations from which it is obtained converge slowly in Monte Carlo simulations. Here we develop two routes that avoid sampling vacancy composition fluctuations. The first constructs the thermodynamic factor from two quantities, the free-energy curvature of the vacancy-free alloy and the vacancy concentration. The second truncates the semi-grand canonical partition function at a single vacancy and recovers the full thermodynamic factor of the vacancy-containing alloy. The same separation shows that the largest eigenvalue of the diffusion matrix is the vacancy tracer diffusion coefficient for any number of components and any degree of non-ideality, while the remaining eigenvalues scale with the vacancy concentration. Across binary through quinary alloys of the Hf--Mo--Nb--Ti--Zr system, both routes reproduce eigenvalues and eigenvectors of diffusion matrices obtained from fully converged Monte Carlo simulations. The single-vacancy expansion reaches this accuracy with roughly one tenth of the sampling effort. A scaling relation extends both routes to local vacancy chemical potentials away from equilibrium, without additional simulations. These results allow thermodynamic descriptions from atomistic models or CALPHAD assessments to be used directly in mesoscale simulations of mass transport.
A ball bouncing on a thicker cushion should bounce less. We show that this is not always so . Measuring the coefficient of restitution $e$ of a steel ball dropped on a stack of $N$ sheets of standard A4 printer paper, we find that $e$ first drops, then rises to a pronounced maximum, and only then decreases towards a minimum. The maximum occurs where the acoustic round trip through the stack matches the contact time of the impact. For a contact set by the stack itself, one recovers the classical bar-impact condition of a stack-to-ball mass ratio of order one. Bar impact, however, predicts only a weak recovery, insensitive to the surrounding gas. For paper stacks, the recovery is large and no maximum is resolved when the stack is evacuated: the air trapped between the sheets, not the paper alone, makes the stack a trampoline.
We report an electron paramagnetic resonance (EPR) study of Dy${^3+}$ in CaWO$_4$ single crystals spanning dopant concentrations from 80~ppb to 247~ppm. The $g$ factors of the ground Kramers doublet at the $S_4$ point-group site, $g_\parallel = 7.22 \pm 0.01$ and $g_\perp = 5.45 \pm 0.01$, and the hyperfine constants of \ce{^{161}Dy} and $^{163}$Dy are determined with an order-of-magnitude improvement in precision over the only previous report. A crystal-field analysis constrained jointly by published optical levels and by the measured $g$ tensor identifies the ground doublet as the lower branch of an anticrossing between the near-degenerate $\|{\pm}11/2>$ and $\|{\mp}13/2>$ states, and explains why optical data alone left the $g$ values unconstrained. \add{The same parameter set, tested on the excited multiplets of the $^6H$ term, reproduces the ${}^6H_{13/2}$ crystal-field levels better than the optical parameters themselves.} The angular dependence of the linewidth in the $ab$ plane reveals broadening by random internal electric fields through the linear electric field effect. Calibrating these fields with the residual Er$^{3+}$ present in the same crystals yields the first electric-field coupling parameters of Dy$^3+$, $(B_{31}^2+B_{36}^2)^{1/2} = (34\pm5)\times10^{-6}$~cm/V, three times the Er$^{3+}$ value and the largest reported for a rare-earth ion in this host.The residual linewidth and the concentration-dependent lineshape asymmetry follow the quadratic field scaling expected for a second-order Stark shift, which escapes the inversion-site cancellation constraining the linear effect. Both couplings trace back to the $\sim 20$~cm$^{-1}$ gap to the first opposite-symmetry doublet, a consequence of the dense crystal-field structure of the $J=15/2$ manifold.