Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “harmonic oscillator”

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

Equilibrium carbon isotope fractionation factors of hydrocarbons: Semi-empirical force-field method

Here, we have calculated the reduced partition function ratios for carbon isotopes (β-factor) of 67 hydrocarbons (alkanes, alkenes, alkynes, cycloalkanes, and aromatics), including their 267 single-substituted isotopomers. The calculations were performed using the harmonic oscillator – rigid rotator model and the Urey / Bigeleisen-Mayer approach. Normal frequencies of molecular vibrations of the isotopologues were calculated from the molecular structures data and valence field force constants stored in the Light-handled Elucidation of Vibrations (LEV) database elaborated by Dr. Gribov and his colleagues in Russia. The LEV database is constructed from experimental spectroscopic and structural data by solving inverse problems of molecular vibrations. The LEV is internally consistent and the most comprehensive database available to date. The β-factor were calculated in the temperature range of 200–800 K with a 10 K step. Our calculations predict that the β-factors increase with an increasing number of C atoms within the same groups of hydrocarbons (e.g., alkanes). Our calculations also show a general descending order of 13 C enrichments among the different groups of hydrocarbons: cycloalkanes, aromatics, alkenes (double bonds) and isoalkanes, alkanes, alkynes (triple bonds). Position-specific, intramolecular isotope effects within hydrocarbons are determined by the β-factors of C in the different functional groups in the order: quaternary (C), tertiary (methine, CH), secondary (methylene, CH 2 ), primary (methyl, CH 3 ) and double bond (C=C), saturated bond (C - C), triple bond (HC≡). Our calculations on bulk and position-specific carbon isotope β-factors of the hydrocarbons, which are generally consistent with very limited ab initio calculations in the literature, are internally consistent and the most comprehensive to date for future applications to position-specific isotope geochemistry of hydrocarbons.

58 GEOSCIENCES↗
DOE OSTI2025

SympGNNs: Symplectic Graph Neural Networks for identifying high-dimensional Hamiltonian systems and node classification

Existing neural network models to learn Hamiltonian systems, such as SympNets, although accurate in low-dimensions, struggle to learn the correct dynamics for high-dimensional many-body systems. Herein, we introduce Symplectic Graph Neural Networks (SympGNNs) that can effectively handle system identification in high-dimensional Hamiltonian systems, as well as node classification. SympGNNs combine symplectic maps with permutation equivariance, a property of graph neural networks. Specifically, we propose two variants of SympGNNs: (i) G-SympGNN and (ii) LA-SympGNN, arising from different parameterizations of the kinetic and potential energy. We demonstrate the capabilities of SympGNN on two physical examples: a 40-particle coupled Harmonic oscillator, and a 2000-particle molecular dynamics simulation in a two-dimensional Lennard-Jones potential. Furthermore, we demonstrate the performance of SympGNN in the node classification task, achieving accuracy comparable to the state-of-the-art. Finally, we also empirically show that SympGNN can overcome the oversmoothing and heterophily problems, two key challenges in the field of graph neural networks.

Deep learning↗
DOE OSTI2023

Emergence of 4 H $J^π$ = 1 - resonance in contact theories

We obtain the s- and p-wave low-energy scattering parameters for n 3 H elastic scattering and the position of the 4 H J π = 1 - resonance using the pionless effective field theory at leading order. Results are extracted with three numerical techniques: confining the system in a harmonic oscillator trap, solving the Faddeev-Yakubovsky equations in configuration space, and using an effective two-body cluster approach. The renormalization of the theory for the relevant amplitudes is assessed in a cutoff-regulator range between 1 fm -1 and 10 fm -1 . Most remarkably, we find a cutoff-stable/RG-invariant resonance in the 4 H J π = 1 - system. This p-wave resonance is a universal consequence of a shallow two-body state and the introduction of a three-body s-wave scale set by the triton binding energy. The stabilization of a resonant state in a few-fermion system through pure contact interactions has a significant consequence for the powercounting of the pionless theory. Specifically, it suggests the appearance of similar resonant states also in larger nuclei, like 16-oxygen, in which the theory’s leading order does not predict stable states. Those resonances would provide a starting state to be moved to the correct physical position by the perturbative insertion of sub-leading orders, possibly resolving the discrepancy between data and contact EFT.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Chaotic and integrable magnetic fields in one-dimensional hybrid Vlasov–Maxwell equilibria

The construction of kinetic equilibrium states is important for studying stability and wave propagation in collisionless plasmas. Thus, many studies over the past decades have been focused on calculating Vlasov–Maxwell equilibria using analytical and numerical methods. However, the problem of kinetic equilibrium of hybrid models is less studied, and self-consistent treatments often adopt restrictive assumptions ruling out cases with irregular and chaotic behaviour, although such behaviour is observed in spacecraft observations of space plasmas. In this paper, we develop a one-dimensional (1-D), quasineutral, hybrid Vlasov–Maxwell equilibrium model with kinetic ions and massless fluid electrons and derive associated solutions. The model allows for an electrostatic potential that is expressed in terms of the vector potential components through the quasineutrality condition. The equilibrium states are calculated upon solving an inhomogeneous Beltrami equation that determines the magnetic field, where the inhomogeneous term is the current density of the kinetic ions and the homogeneous term represents the electron current density. We show that the corresponding 1-D system is Hamiltonian, with position playing the role of time, and its trajectories have a regular, periodic behaviour for ion distribution functions that are symmetric in the two conserved particle canonical momenta. For asymmetric distribution functions, the system is nonintegrable, resulting in irregular and chaotic behaviour of the fields. The electron current density can modify the magnetic field phase space structure, inducing orbit trapping and the organization of orbits into large islands of stability. Thus, the electron contribution can be responsible for the emergence of localized electric field structures that induce ion trapping. We also provide a paradigm for the analytical construction of hybrid equilibria using a rotating two-dimensional harmonic oscillator Hamiltonian, enabling the calculation of analytic magnetic fields and the construction of the corresponding distribution functions in terms of Hermite polynomials.

Kaltsas, Dimitrios A. (ORCID:0000000300769015)↗
DOE OSTI2020

Density-Matrix Based Extended Lagrangian Born–Oppenheimer Molecular Dynamics

Extended Lagrangian Born–Oppenheimer molecular dynamics [ Phys. Rev. Lett. 2008, 100, 123004] is presented for Hartree–Fock theory, where the extended electronic degrees of freedom are represented by a density matrix, including fractional occupation numbers at elevated electronic temperatures. In contrast to regular direct Born–Oppenheimer molecular dynamics simulations, no iterative self-consistent field optimization is required prior to the force evaluations. To sample regions of the potential energy landscape where the gap is small or vanishing, which leads to particular convergence problems in regular direct Born–Oppenheimer molecular dynamics simulations, an adaptive integration scheme for the extended electronic degrees of freedom is presented. The integration scheme is based on a tunable, low-rank approximation of a fourth-order kernel, K, that determines the metric tensor, T ≡ K T K, used in the extended harmonic oscillator of the Lagrangian that generates the dynamics of the electronic degrees of freedom. Here, the formulation and algorithms provide a general guide to implement extended Lagrangian Born–Oppenheimer molecular dynamics for quantum chemistry, density functional theory, and semiempirical methods using a density matrix formalism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Entanglement Effect and Angular Momentum Conservation in a Nonseparable Tunneling Treatment

The important, and often dominant, role of tunneling in low temperature kinetics has resulted in numerous theoretical explorations into the methodology for predicting it. Nevertheless, there are still key aspects of the derivations that are lacking, particularly for nonseparable systems in the low temperature regime, and further explorations of the physical factors affecting the tunneling rate are warranted. In this work we obtain a closed-form rate expression for the tunneling rate constant that is a direct analog of the rigid-rotor-harmonic-oscillator expression. This expression introduces a novel "entanglement factor" that modulates the reaction rate. Furthermore, we are able to extend this expression, which is valid for nonseparable systems at low temperatures, to properly account for the conservation of angular momentum. In contrast, previous calculations have considered only vibrational transverse modes and so effectively employ a decoupled rotational partition function for the orientational modes. We also suggest a simple theoretical model to describe the tunneling effects in the vicinity of the crossover temperature (the temperature where tunneling becomes the dominating mechanism). This model allows one to naturally classify, interpret, and predict experimental data. Among other things, it quantitatively explains in simple terms the so-called "quantum bobsled" effect, also known as the negative centrifugal effect, which is related to curvature of the reaction path. Taken together, the expressions obtained here allow one to predict the thermal and E-resolved rate constants over broad ranges of temperatures and energies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Fermi’s Golden Rule Rate Expression for Transitions Due to Nonadiabatic Derivative Couplings in the Adiabatic Basis

Starting from a general molecular Hamiltonian expressed in the basis of adiabatic electronic and nuclear position states, where a compact and complete expression for the nonadiabatic derivative coupling (NDC) Hamiltonian term is obtained, we provide a general analysis of the Fermi’s golden rule (FGR) rate expression for nonadiabatic transitions between adiabatic states. We then consider a quasi-adiabatic approximation that uses crude adiabatic states and NDC couplings, both evaluated at the minimum potential energy configuration of the initial adiabatic state, for the definition of the zeroth and first-order terms of the Hamiltonian. Although the application of this approximation is rather limited, it allows deriving a general FGR rate expression without further approximation while accounting for non-Condon contribution to the FGR rate arising from momentum operators of NDC terms and its coupling with vibronic displacements. For a generic and widely used model where all nuclear degrees of freedom and environmental effects are represented as linearly coupled harmonic oscillators, we derive a closed-form FGR rate expression that requires only Fourier transform. The resulting rate expression includes quadratic contributions of NDC terms and their couplings to Franck–Condon modes, which require evaluation of two additional bath spectral densities in addition to the conventional one that appears in a typical FGR rate theory based on the Condon approximation. Model calculations for the case where nuclear vibrations consist of both a sharp high-frequency mode and an Ohmic bath spectral density illustrate new features and implications of the rate expression. We then apply our theoretical expression to the nonradiative decay from the first excited singlet state of azulene, which illustrates the utility and implications of our theoretical results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Quantum Simulation of Molecular Dynamics Processes─A Benchmark Study Using a Classical Simulator and Present-Day Quantum Hardware

Here, we explore how the fundamental problems in quantum molecular dynamics can be modeled using classical simulators (emulators) of quantum computers and the actual quantum hardware available to us today. The list of problems we tackle includes propagation of a free wave packet, vibration of a harmonic oscillator, and tunneling through a barrier. Each of these problems starts with the initial wave packet setup. Although Qiskit provides a general method for initializing wave functions, in most cases it generates deep quantum circuits. While these circuits perform well on noiseless simulators, they suffer from excessive noise on quantum hardware. To overcome this issue, we designed a shallower quantum circuit for preparing a Gaussian-like initial wave packet, which improves the performance of real hardware. Next, quantum circuits are implemented to apply the kinetic and potential energy operators for the evolution of a wave function over time. The results of our modeling on classical emulators of quantum hardware agree perfectly with the results obtained using the traditional (classical) methods. This serves as a benchmark and demonstrates that the quantum algorithms and Qiskit codes we developed are accurate. However, the results obtained on the actual quantum hardware available today, such as IBM’s superconducting qubits and IonQ’s trapped ions, indicate large discrepancies due to hardware limitations. This work highlights both the potential and challenges of using quantum computers to solve fundamental quantum molecular dynamics problems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Effect of Collective Dynamics and Anharmonicity on Entropy in Heterogenous Catalysis: Building the Case for Advanced Molecular Simulations

In this paper, we present a perspective on the computational determination of entropy and its effects and consequences on heterogeneous catalysis. Special attention is paid to the role of anharmonicity and the deviations from the standard harmonic oscillator approximations which can fail to provide a reliable assessment of entropy. To address these challenges, advanced methodologies are needed that can explicitly account for these thermodynamic drivers through the appropriate statistical sampling of reactive free energy surfaces. We discuss where anharmonicity should be expected, where it has been observed from a theoretical perspective, and the methods currently employed to address it. We concentrate on three types of systems where we have observed major, nonnegligible anharmonic effects: (1) supported nanoparticles, where the migration of metal atoms, complexes, and entire clusters exhibit anharmonic behavior in their dynamic motion; (2) porous solids where confinement effects distort potential energy surfaces and hinder molecular motions resulting in large entropic terms; and (3) solid/liquid interfaces, where interactions between solvent molecules and adsorbed species can result in large solvent organization free energy and unique reactivity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Investigation of Ethane Dehydrogenation and Hydrogenolysis on Pt(111), Pt(211), and Pt(100): Bayesian Quantification and Correction of DFT-Based Enthalpic and Entropic Uncertainties

Computational investigations of heterogeneously catalyzed reactions using density functional theory (DFT) are often inaccurate, largely due to uncertainties in the choice of DFT functional (enthalpic uncertainty) and approximations for modeling adsorbate movement along the catalyst surface (entropic uncertainty). This work illustrates that both uncertainties are significant in the investigation of ethane dehydrogenation (EDH) and hydrogenolysis on Pt catalysts by considering the complete deconstruction of ethane on Pt(111), Pt(211), and Pt(100) using microkinetic modeling (MKM). Hence, this work uses both noncalibrated and Bayesian-calibrated MKMs to quantify and correct inaccuracies in macroscopic properties due to both uncertainties. A Bayesian approach to the correction of entropic errors was introduced using a “Modified Fermi Function (MFF)” to calibrate between the two bounds of entropy represented by the harmonic oscillator (HO) and free translator (FT) approximations. Regardless of enthalpic and entropic uncertainties, all three surfaces are capable of ethane activation; however, Pt(211) was found to be the most active and is largely responsible for methane production. Next, Pt(111) is largely responsible for acetylene production, and Pt(100) has the highest ethylene selectivity but is most susceptible to coking. By comparison of different calibrated models, the FT entropy approximation was found to better describe EDH under typical experimental conditions. Statistical evidence was found to support Pt(111) as the active site for EDH, assuming that one single site is responsible for the chemistry. On the three surfaces, competing second dehydrogenations to CH 2 CH 2 and CH 3 CH were observed as well as isomerization of CH 3 CH back to CH 2 CH 2 and deeper dehydrogenation of CH 3 CH. In conclusion, C–C cleavage was found to largely proceed via the CH 3 C intermediate on Pt(100) and Pt(111), while on Pt(211), it was via both CHC and CH 3 C.

Bayesian model selection↗
DOE OSTI2020

Plasmon-Mediated Coherent Superposition of Discrete Excitons under Strong Exciton–Plasmon Coupling in Few-Layer MoS 2 at Room Temperature

In this paper, we demonstrate room temperature coherent hybridization of the A- and B-excitons in few-layer MoS 2 , mediated by simultaneous strong coupling to surface plasmon polaritons. Few-layer MoS 2 was placed on a tunable plasmonic structure and the system’s dispersion was measured by tuning the plasmon energy across the exciton energies. Strong coupling was observed as double Rabi splitting at the A- and B-excitons of 81 and 93 meV, respectively. A coupled harmonic oscillator model sheds light on the nature of the interaction, revealing a quantum superposition of the A- and B-excitons, mediated by the plasmon interaction. This observation suggests the possibility of room temperature intra- or intervalley quantum information transport and/or spin entanglement. The experiment confirms a previous theoretical prediction of room temperature exciton–exciton hybridization in two-dimensional MoS 2 . Further, through modeling we find that room temperature strong coupling is a general phenomenon among two-dimensional transition metal dichalcogenide exciton–plasmon systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Molecular-Level Insights into the NMR Relaxivity of Gadobutrol Using Quantum and Classical Molecular Simulations

MRI is an indispensable diagnostic tool in modern medicine; however, understanding the molecular-level processes governing NMR relaxation of water in the presence of MRI contrast agents remains a challenge, hindering the molecular-guided development of more effective contrast agents. By using quantum-based polarizable force fields, the first-of-its-kind molecular dynamics (MD) simulations of Gadobutrol are reported where the 1 H NMR longitudinal relaxivity r 1 of the aqueous phase is determined without any adjustable parameters. The MD simulations of r 1 dispersion (i.e., frequency dependence) show good agreement with measurements at frequencies of interest in clinical MRI. Importantly, the simulations reveal key insights into the molecular level processes leading to r 1 dispersion by decomposing the NMR dipole–dipole autocorrelation function G(t) into a discrete set of molecular modes, analogous to the eigenmodes of a quantum harmonic oscillator. The molecular modes reveal important aspects of the underlying mechanisms governing r 1 , such as its multiexponential nature and the importance of the second eigenmodal decay. By simply analyzing the MD trajectories on a parameter-free approach, the Gadobutrol simulations show that the outer-shell water contributes ∼50% of the total relaxivity r 1 compared to the inner-shell water, in contrast to simulations of (nonchelated) gadolinium-aqua where the outer shell contributes only ∼15% of r 1 . The deviation between simulations and measurements of r 1 below clinical MRI frequencies is used to determine the low-frequency electron-spin relaxation time for Gadobutrol, in good agreement with independent studies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Environmental damping and vibrational coupling of confined fluids within isolated carbon nanotubes

Abstract Because of their large surface areas, nanotubes and nanowires demonstrate exquisite mechanical coupling to their surroundings, promising advanced sensors and nanomechanical devices. However, this environmental sensitivity has resulted in several ambiguous observations of vibrational coupling across various experiments. Herein, we demonstrate a temperature-dependent Radial Breathing Mode (RBM) frequency in free-standing, electron-diffraction-assigned Double-Walled Carbon Nanotubes (DWNTs) that shows an unexpected and thermally reversible frequency downshift of 10 to 15%, for systems isolated in vacuum. An analysis based on a harmonic oscillator model assigns the distinctive frequency cusp, produced over 93 scans of 3 distinct DWNTs, along with the hyperbolic trajectory, to a reversible increase in damping from graphitic ribbons on the exterior surface. Strain-dependent coupling from self-tensioned, suspended DWNTs maintains the ratio of spring-to-damping frequencies, producing a stable saturation of RBM in the low-tension limit. In contrast, when the interior of DWNTs is subjected to a water-filling process, the RBM thermal trajectory is altered to that of a Langmuir isobar and elliptical trajectories, allowing measurement of the enthalpy of confined fluid phase change. These mechanisms and quantitative theory provide new insights into the environmental coupling of nanomechanical systems and the implications for devices and nanofluidic conduits.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Geometric control of hyperbolic exciton-polariton condensate dimers

Coupled many-body quantum systems exhibit rich emergent physics with diverse stationary and dynamical behaviours. By engineering platforms with tunable and distinct coupling mechanisms, new insights emerge into the collective behaviour of coupled many body systems. Particles can be exchanged via evanescent or ballistic coupling: the former, based on proximity, yields large spectral splitting, while the latter requires strict phase-matching, analogous to phase-coupled harmonic oscillators and has a smaller impact on the energy landscape. We demonstrate an all-optically tunable quantum fluid dimer based on exciton-polariton condensates in a photonic crystal waveguide with hyperbolic (saddle-like) dispersion. Varying the dimer’s angle relative to the grating tunes the coupling from evanescent to ballistic. We directly observe spectral features and mass flow shaped by the saddle dispersion. This work highlights photonic crystals as powerful platforms to explore condensed matter phenomena lying at the interface between delay-coupled nonlinear oscillators and tight binding physics.

Georgakilas, Ioannis [IBM Research-Zurich, Rüschli↗
DOE OSTI2023

Deep learning for exploring ultra-thin ferroelectrics with highly improved sensitivity of piezoresponse force microscopy

Hafnium oxide-based ferroelectrics have been extensively studied because of their existing ferroelectricity, even in ultra-thin film form. However, studying the weak response from ultra-thin film requires improved measurement sensitivity. In general, resonance-enhanced piezoresponse force microscopy (PFM) has been used to characterize ferroelectricity by fitting a simple harmonic oscillation model with the resonance spectrum. However, an iterative approach, such as traditional least squares (LS) fitting, is sensitive to noise and can result in the misunderstanding of weak responses. In this study, we developed the deep neural network (DNN) hybrid with deep denoising autoencoder (DDA) and principal component analysis (PCA) to extract resonance information. The DDA/PCA-DNN improves the PFM sensitivity down to 0.3 pm, allowing measurement of weak piezoresponse with low excitation voltage in 10-nm-thick Hf 0.5 Zr 0.5 O 2 thin films. Our hybrid approaches could provide more chances to explore the low piezoresponse of the ultra-thin ferroelectrics and could be applied to other microscopic techniques.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Extended Lagrangian Born–Oppenheimer molecular dynamics using a Krylov subspace approximation

It is shown how the electronic equations of motion in extended Lagrangian Born–Oppenheimer molecular dynamics simulations can be integrated using low-rank approximations of the inverse Jacobian kernel. This kernel determines the metric tensor in the harmonic oscillator extension of the Lagrangian that drives the evolution of the electronic degrees of freedom. The proposed kernel approximation is derived from a pseudoinverse of a low-rank estimate of the Jacobian, which is expressed in terms of a generalized set of directional derivatives with directions that are given from a Krylov subspace approximation. The approach allows a tunable and adaptive approximation that can take advantage of efficient preconditioning techniques. The proposed kernel approximation for the integration of the electronic equations of motion makes it possible to apply extended Lagrangian first-principles molecular dynamics simulations to a broader range of problems, including reactive chemical systems with numerically sensitive and unsteady charge solutions. This can be achieved without requiring exact full calculations of the inverse Jacobian kernel in each time step or relying on iterative non-linear self-consistent field optimization of the electronic ground state prior to the force evaluations as in regular direct Born–Oppenheimer molecular dynamics. We note the low-rank approximation of the Jacobian is directly related to Broyden’s class of quasi-Newton algorithms and Jacobian-free Newton–Krylov methods and provides a complementary formulation for the solution of nonlinear systems of equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Ultrasonic oscillatory two-phase flow in microchannels

Experimental and numerical investigations are performed to provide an assessment of the transport behavior of an ultrasonic oscillatory two-phase flow in a microchannel. The work is inspired by the flow observed in an innovative ultrasonic fabric drying device using a piezoelectric bimorph transducer with microchannels, where a water-air two-phase flow is transported by harmonically oscillating microchannels. The flow exhibits highly unsteady behavior as the water and air interact with each other during the vibration cycles, making it significantly different from the well-studied steady flow in microchannels. Here, the computational fluid dynamics (CFD) modeling is realized by combing the turbulence Reynolds-averaged Navier-Stokes (RANS) k – ω model with the phase-field method to resolve the dynamics of the two-phase flow. The numerical results are qualitatively validated by the experiment. Through parametric studies, we specifically examined the effects of vibration conditions (i.e., frequency and amplitude), microchannel taper angle, and wall surface contact angle (i.e., wettability) on the flow rate through the microchannel. The results will advance the potential applications where oscillatory or general unsteady microchannel two-phase flows may be present.

42 ENGINEERING↗
DOE OSTI2022

Charged particle motion in spherically symmetric distributions of magnetic monopoles

The classical equations of motion of a charged particle in a spherically symmetric distribution of magnetic monopoles can be transformed into a system of linear equations, thereby providing a type of integrability. In the case of a single monopole, the solution was given long ago by Poincaré. In the case of a uniform distribution of monopoles, the solution can be expressed in terms of parabolic cylinder functions (essentially the eigenfunctions of an inverted harmonic oscillator). Further, this solution is relevant to recent studies of nonassociative star products, symplectic lifts of twisted Poisson structures, and fluids and plasmas of electric and magnetic charges.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗