Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “laser cluster”

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 109 records · Page 6

Two Cycling Centers in One Molecule: Communication by Through-Bond Interactions and Entanglement of the Unpaired Electrons

Many applications in quantum information science (QIS) rely on the ability to laser- cool molecules. The scope of applications can be expanded if laser-coolable molecules possess two or more cycling centers, i.e., moieties capable of scattering photons via multiple absorption-emission events. Here in this paper we employ equation-of-motion coupled-cluster method for double electron attachment (EOM-DEA-CCSD) to study the electronic structure of hypermetallic molecules with two alkaline earth metals connected by an acetylene linker. The electronic structure of the molecules is similar to that of two separated alkali metals; however, the interaction between the two electrons is weak and largely dominated by through-bond interactions. The communication between the two cycling centers is quantified by the extent of the entanglement of the two unpaired electrons associated with each center. This contribution highlights the rich electronic structure of hypermetallic molecules that may advance various applications in QIS and beyond.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Conservation laws in coupled cluster dynamics at finite temperature

We extend the finite-temperature Keldysh non-equilibrium coupled cluster theory (Keldysh-CC) [A. F. White and G. K.-L. Chan, J. Chem. Theory Comput. 15, 6137–6253 (2019)] to include a time-dependent orbital basis. When chosen to minimize the action, such a basis restores local and global conservation laws (Ehrenfest’s theorem) for all one-particle properties while remaining energy conserving for time independent Hamiltonians. We present the time-dependent Keldysh orbital-optimized coupled cluster doubles method in analogy with the formalism for zero-temperature dynamics, extended to finite temperatures through the time-dependent action on the Keldysh contour. To demonstrate the conservation property and understand the numerical performance of the method, we apply it to several problems of nonequilibrium finite-temperature dynamics: a 1D Hubbard model with a time-dependent Peierls phase, laser driving of molecular H2, driven dynamics in warm-dense silicon, and transport in the single impurity Anderson model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular dynamics simulations of radiation cascade evolution near cellular dislocation structures in additively manufactured stainless steels

The dislocation dense cellular structures observed in additively manufactured (AM) stainless steels can allow these materials to have enhanced strength and ductility compared to conventionally manufactured materials. Effective design of radiation resistant materials often necessitates creating microstructures containing of a high density of dislocations that act as sinks for radiation induced defects. This work uses molecular dynamics simulations to study the impact of pre-existing defects on radiation damage in stainless steel 316L fabricated by the laser powder bed fusion process. The evolution of synthetically generated dislocations in a crystal system in response to multiple radiation collision cascades suggest that dislocation dense regions reduce the probability of surviving point defects forming clusters. This finding seems to result from the dislocation cores’ ability to absorb point defects, thereby limiting the formation of clusters during the recombination phase. However, after successive cascades, the ability of the dislocation entanglement to prevent defect formation is diminished and defect cluster formation trends to an equilibrium. Accordingly, it is hypothesized that in an experimental setting the dislocation cells will initially act as neutral sinks for point defects, potentially delaying the onset of radiation damage effects; but, these structures will be degraded due to radiation-enhanced diffusion and lose their effectiveness at doses above 1.5–2 displacements per atom. Nevertheless, the ability to produce complex and feature-specific microstructures suggests promise for the design of radiation tolerant materials using AM methods.

36 MATERIALS SCIENCE↗

Deciphering Photochemical Reaction Pathways of Aqueous Tetrachloroauric Acid by X-ray Transient Absorption Spectroscopy

Photolysis reaction pathways of [Au(III)Cl 4 ] - in aqueous solution have been investigated by time-resolved X-ray absorption spectroscopy. Ultraviolet excitation directly breaks the Au-Cl bond in [Au(III)Cl 4 ] - to form [Au(II)Cl 3 ] - that becomes highly reactive within 79 ps. Disproportionation of [Au(II)Cl 3 ] - generates [Au(I)Cl 2 ] - , which is stable for ≤ 10 mu s. In contrast, intense near-infrared lasers photolyze water to generate hydrated electrons, which then reduce [Au(III)Cl 4 ] - to [Au(II)Cl 3 ] - at 5 ns. Hydrated electrons further induce a chain reaction from [Au(II)Cl 3 ] - to [Au(0)Cl] - by successively removing one Cl - . The zero-valency Au anions quickly polymerize and condense to form Au nanoparticles, which become the dominating product after 400 s. Here our results reveal that the condensation of zero-valency Au starts with dimerization of gold clusters coordinated with chloride ions rather than direct condensation of pristine Au atoms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Cluster Dynamics Modeling Needs for the Advanced Materials and Manufacturing Technologies Program

This milestone report aims to identify and assess the cluster dynamics (CD) modeling requirements within the Department of Energy's Office of Nuclear Energy (DOE-NE) Advanced Materials and Manufacturing Technologies (AMMT) program and to communicate these needs to the DOE-NE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. The goal is to ensure NEAMS is well-informed about the CD modeling requirements to support AMMT's mission of accelerating the development, qualification, demonstration, and deployment of advanced structural materials and manufacturing for nuclear energy applications. CD modeling is an essential tool for predicting the degradation of structural materials under irradiation, which is a key component of AMMT's accelerated qualification process. The AMMT program focuses on both additively manufactured and wrought structural alloys, such as laser powder-bed fusion 316H austenitic stainless steel, alloy 709, Haynes 244, and alloy 617. These materials require a generalized CD modeling framework to facilitate rapid model development and computational simulation. A flexible, generalized CD software, similar to the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element framework, would enable modeling of various cluster types, including defect clusters, defect-solute clusters, and multicomponent clusters, incorporating thermodynamics and kinetics parameters. Radiation effects, microstructural feature evolution, and multi-dimensional modeling are critical considerations for the CD model. The usability of the CD code should allow for easy modification and coupling with MOOSE-based simulations. Additionally, the software should adhere to Nuclear Quality Assurance-1 standards, include a testing suite for verification and validation, and be version-controlled within a national laboratory-managed Git repository. Benchmark problems are needed to assess code predictions and performance.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Performance Analysis of PIConGPU: Particle-in-Cell on GPUs using NVIDIA’s NSight Systems and NSight Compute

PIConGPU, Particle In Cell on GPUs, is an open source simulations framework for plasma and laser-plasma physics used to develop advanced particle accelerators for radiation therapy of cancer, high energy physics and photon science. While PIConGPU has been optimized for at least 5 years to run well on NVIDIA GPU-based clusters, there has been limited exploration by the development team of potential scalability bottlenecks using recently updated and new tools including NVIDIA’s NVProf tool and the brand-new NVIDIA NSight Suite (Systems and Compute) tools. PIConGPU is a highly optimized application that runs production jobs at scale on a system Oak Ridge Leadership Facility’s (OLCF) Summit supercomputer (using the full machine at 4600 nodes; at 98% of GPU utilization on all ~28000 NVIDIA Volta GPUs). PIConGPU has been selected as one of the the eight applications for OLCF’s coveted Center for Accelerated Application Readiness (CAAR) program aimed at the facility’s Frontier supercomputer (OLCF’s first exascale system to launch in 2021), to partner with our vendors (primary vendors: AMD and Cray/HPE) ensuring that Frontier will be able to perform large-scale science when it opens to users in 2022. To this effect, performance engineers on the PIConGPU team wanted to dive deep into the application to understand at the finest granularity, which portions of the code could be further optimized to exploit the hardware on Summit at it’s maximum potential and also to elucidate which key kernels should be tracked and optimized for the CAAR effort to port this code to Frontier. Any bottlenecks that are observed via performance profiling on Summit are likely to also impact scalability on the Frontier-dev system and the Frontier Early Access (EA) system. Additionally, the engineers wanted to take a closer look at the newest NVIDIA profiling tools which allows us to identify the most useful features on these tools and will provide an opportunity to compare it to new AMD and Cray’s performance analysis tool releases and provide feedback to our vendor partners on what features are most important and mission critical for CAAR efforts. The primary goal of this report is to focus on the evaluation of PIConGPU’s most time-intensive kernels using NVProf and NSight Suite. Three kernels, Current Deposition (also known as Compute Current), Particle Push (Move and Mark), and Shift Particles are known to be some of the most time-consuming kernels in PIConGPU. The Current xi Deposition kernel and Particle Push kernel both set up the particle attributes for running any physics simulation with PIConGPU, so it is crucial to improve the performance of these two kernels. In this report, we measure single GPU metrics for the three kernels, offer high level takeaways from the conducted analysis, and compare the profiling data from NSight Compute to that of NVProf. This analysis was performed using a grid size of 240 x 272 x 224, and 10 time steps with the Mid-November Figure of Merit (FOM) run setup. The Traveling Wave Electron Acceleration (TWEAC) science case used in this run is a representative science case for PIConGPU. This execution can also be used for baseline analysis on AMD MI50/ MI60 systems. As of the time of writing, the PIConGPU application has limited use for features of NSight Systems, so this report will mainly focus on insights garnered from NSight Compute. For this analysis, we run the “full” metric set available in NSight Compute version 2020.1.2 and use NSight Systems version 2020.3.1 to generate the application timeline.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Uncovering acoustic signatures of pore formation in laser powder bed fusion

Abstract We present a machine learning workflow to discover signatures in acoustic measurements that can be utilized to create a low-dimensional model to accurately predict the location of keyhole pores formed during additive manufacturing processes. Acoustic measurements were sampled at 100 kHz during single-layer laser powder bed fusion (LPBF) experiments, and spatio-temporal registration of pore locations was obtained from post-build radiography. Power spectral density (PSD) estimates of the acoustic data were then decomposed using non-negative matrix factorization with custom $$\varvec{k}$$ k -means clustering (NMF $$\varvec{k}$$ k ) to learn the underlying spectral patterns associated with pore formation. NMF $$\varvec{k}$$ k returned a library of basis signals and matching coefficients to blindly construct a feature space based on the PSD estimates in an optimized fashion. Moreover, the NMF $$\varvec{k}$$ k decomposition led to the development of computationally inexpensive machine learning models which are capable of quickly and accurately identifying pore formation with classification accuracy of supervised and unsupervised label learning greater than 95% and 90%, respectively. The intrinsic data compression of NMF k , the relatively light computational cost of the machine learning workflow, and the high classification accuracy makes the proposed workflow an attractive candidate for edge computing toward in-situ keyhole pore prediction in LPBF.

36 MATERIALS SCIENCE↗

High-Resolution XFEL Structure of the Soluble Methane Monooxygenase Hydroxylase Complex with its Regulatory Component at Ambient Temperature in Two Oxidation States

Soluble methane monooxygenase (sMMO) is a multicomponent metalloenzyme that catalyzes the conversion of methane to methanol at ambient temperature using a nonheme, oxygen-bridged dinuclear iron cluster in the active site. Structural changes in the hydroxylase component (sMMOH) containing the diiron cluster caused by complex formation with a regulatory component (MMOB) and by iron reduction are important for the regulation of O 2 activation and substrate hydroxylation. Structural studies of metalloenzymes using traditional synchrotron-based X-ray crystallography are often complicated by partial X-ray-induced photoreduction of the metal center, thereby obviating determination of the structure of the enzyme in pure oxidation states. Here, microcrystals of the sMMOH:MMOB complex from Methylosinus trichosporium OB3b were serially exposed to X-ray free electron laser (XFEL) pulses, where the ≤35 fs duration of exposure of an individual crystal yields diffraction data before photoreduction-induced structural changes can manifest. Merging diffraction patterns obtained from thousands of crystals generates radiation damage-free, 1.95 Å resolution crystal structures for the fully oxidized and fully reduced states of the sMMOH:MMOB complex for the first time. The results provide new insight into the manner by which the diiron cluster and the active site environment are reorganized by the regulatory protein component in order to enhance the steps of oxygen activation and methane oxidation. Furthermore, this study also emphasizes the value of XFEL and serial femtosecond crystallography (SFX) methods for investigating the structures of metalloenzymes with radiation sensitive metal active sites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Discriminative Dimensionality Reduction using Deep Neural Networks for Clustering of LIGO Data

In this paper, leveraging the capabilities of neural networks for modeling the non-linearities that exist in the data, we propose several models that can project data into a low dimensional, discriminative, and smooth manifold. The proposed models can transfer knowledge from the domain of known classes to a new domain where the classes are unknown. A clustering algorithm is further applied in the new domain to find potentially new classes from the pool of unlabeled data. The research problem and data for this paper originated from the Gravity Spy project which is a side project of Advanced Laser Interferometer Gravitational-wave Observatory (LIGO). The LIGO project aims at detecting cosmic gravitational waves using huge detectors. However non-cosmic, non-Gaussian disturbances known as "glitches", show up in gravitational-wave data of LIGO. This is undesirable as it creates problems for the gravitational wave detection process. Gravity Spy aids in glitch identification with the purpose of understanding their origin. Since new types of glitches appear over time, one of the objective of Gravity Spy is to create new glitch classes. Towards this task, we offer a methodology in this paper to accomplish this.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Affinity propagation clustering of full-field, high-spatial-dimensional measurements for robust output-only modal identification: A proof-of-concept study

Determination of the model order is a challenging problem in system identification, especially in output-only or operational modal identification where some modes are weakly excited. While existing methods such as the stabilization diagram method (spectral information) are effective, they do not scale to high-dimensional data, which is usually needed for high-fidelity characterization of structural dynamics and has been made available in the emerging full-field measurement techniques using optical methods such as photogrammetry and laser vibrometers. In this proof-of-concept study we present a new non-parametric, data-driven approach for robust output-only identification of high-spatial-dimensional modal parameters of basic structures by efficiently processing and interactively exploiting the full-field measurement (i.e., very dense spatial measurement points). Specifically, we first over-estimate the system model once, producing a pool of candidate modes associated with their modal frequencies and full-field, high-spatial-dimensional mode shapes. This is accomplished by a data-driven method termed affinity propagation clustering (APC), where the active clusters, which are the active modes in our formulations, emerge from the "message-passing" procedure and does not require a pre-determination of the cluster number (mode or model order). Next, rather than using the spectral information to distinguish the physical and spurious modes in the stabilization diagram method, we exploit and visualize the spatial, full-field mode shape associated with each candidate mode to do so. We conduct extensive experiments on basic structural models with comparisons to a few existing methods. The results indicate that the new method is computationally efficient for identifying high-spatial-dimensional modal parameters, and robust to identify weak modes by exploiting the full-field measurement. We also discuss its applicability and limitations for structures with complex geometry (shapes).

42 ENGINEERING↗

Peculiarities of the formation of photonic nanojets by a matrix of dielectric microtoroids

We report the results of theoretical modelling of photonic nanojets (PNJs) formed by laser radiation scattering on a single-layer ordered assembly of dielectric microtoroids placed on the surface of a transparent matrix (silicone film). Using the method of computational electrodynamics (FDTD), the main PNJ parameters (length, width, peak intensity) are analysed under the conditions of mutual influence of the light fields of neighbouring microparticles. It is shown that the main factor affecting the PNJ characteristics under study is the spatial configuration of the radiation-scattering particle, namely, the internal diameter of its cross section. It is found that for certain configurations of toroid placement in a cluster, a PNJ ensemble is implemented with parameters significantly better than those for a single toroid. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

SEED LIPID DROPLET PROTEIN1, SEED LIPID DROPLET PROTEIN2, and LIPID DROPLET PLASMA MEMBRANE ADAPTOR mediate lipid droplet–plasma membrane tethering

Membrane contact sites (MCSs) are interorganellar connections that allow for the direct exchange of molecules, such as lipids or Ca 2+ between organelles, but can also serve to tether organelles at specific locations within cells. Here, we identified and characterized three proteins of Arabidopsis thaliana that form a lipid droplet (LD)–plasma membrane (PM) tethering complex in plant cells, namely LD-localized SEED LD PROTEIN (SLDP) 1 and SLDP2 and PM-localized LD-PLASMA MEMBRANE ADAPTOR (LIPA). Using proteomics and different protein–protein interaction assays, we show that both SLDPs associate with LIPA. Disruption of either SLDP1 and SLDP2 expression, or that of LIPA, leads to an aberrant clustering of LDs in Arabidopsis seedlings. Ectopic co-expression of one of the SLDPs with LIPA is sufficient to reconstitute LD–PM tethering in Nicotiana tabacum pollen tubes, a cell type characterized by dynamically moving LDs in the cytosolic streaming. Furthermore, confocal laser scanning microscopy revealed both SLDP2.1 and LIPA to be enriched at LD–PM contact sites in seedlings. These and other results suggest that SLDP and LIPA interact to form a tethering complex that anchors a subset of LDs to the PM during post-germinative seedling growth in Arabidopsis.

59 BASIC BIOLOGICAL SCIENCES↗

Regional surrogates for predictive control of digital twins

Digital twins of complex systems must involve a model that is fast, generalizable, and usable for real-time control. For example, high-fidelity nonlinear multiphysics simulations can capture laser-material interactions, but are too slow for optimization or model predictive control (MPC). Reduced-order models, used to accelerate such computation, frequently fail to generalize to unseen inputs or control states. We show theoretically that this failure is intrinsic, i.e., that a learned model is non-unique outside the sampled subspace when its low-rank structure arises from limited excitation and clustered eigenvalues, rather than from a user-imposed truncation alone. Motivated by this result, we propose a control-ready regional surrogate-construction framework for both autonomous and nonautonomous dynamics; it employs Koopman lifting to represent nonlinearities, while preserving spatial locality. We illustrate our approach by constructing a control-ready surrogate for the digital twin of a thermal component of additive-manufacturing process. Our surrogate, localized in space through a von Neumann stencil, is learned from noisy high-fidelity simulations that emulate thermal-camera images collected during the manufacturing. It is linear in thermo-physically augmented states so that MPC reduces to a convex quadratic program. The surrogate requires no online correction, generalizes to unseen scan paths and power profiles of the laser, and is more than three orders of magnitude faster than a finite-difference solver. Furthermore, when the MPC sequence computed on the digital twin is applied to this solver, closed-loop temperature regulation is recovered, showing that the surrogate preserves control-relevant input-output behavior.

Data-driven model↗

Structural insights into the mechanism of rhodopsin phosphodiesterase

Crystal structures of TMD of Rhodopsin phosphodiesterase (Rh-PDE). 6160 small-wedge (10°/crystal) datasets collected from loop-harvested microcrystals using EIGER X 9M detector at a wavelength of 1 Å on BL32XU, SPring-8. Beam size was around 15×10 µm2 and oscillation step was 0.1°. The crystals belonged to space group P212121 with unit cell parameters a~66, b~74, c~117 Å. 864 datasets were indexed and integrated using XDS with KAMO pipeline. Finally 271 integrated results were merged at 2.6 Å resolution after CC-based clustering and outlier rejection by KAMO in the published result.

BL32XU↗

Metallofullerene photoswitches driven by photoinduced fullerene-to-metal electron transfer

We report on the discovery and detailed exploration of the unconventional photo-switching mechanism in metallofullerenes, in which the energy of the photon absorbed by the carbon cage π-system is transformed to mechanical motion of the endohedral cluster accompanied by accumulation of spin density on the metal atoms. Comprehensive photophysical and electron paramagnetic resonance (EPR) studies augmented by theoretical modelling are performed to address the phenomenon of the light-induced photo-switching and triplet state spin dynamics in a series of Y x Sc 3-x N@C 80 (x = 0–3) nitride clusterfullerenes. Variable temperature and time-resolved photoluminescence studies revealed a strong dependence of their photophysical properties on the number of Sc atoms in the cluster. All molecules in the series exhibit temperature-dependent luminescence assigned to the near-infrared thermally-activated delayed fluorescence (TADF) and phosphorescence. The emission wavelengths and Stokes shift increase systematically with the number of Sc atoms in the endohedral cluster, whereas the triplet state lifetime and S 1 –T 1 gap decrease in this row. For Sc 3 N@C 80 , we also applied photoelectron spectroscopy to obtain the triplet state energy as well as the electron affinity. Spin distribution and dynamics in the triplet states are then studied by light-induced pulsed EPR and ENDOR spectroscopies. The spin–lattice relaxation times and triplet state lifetimes are determined from the temporal evolution of the electron spin echo after the laser pulse. Well resolved ENDOR spectra of triplets with a rich structure caused by the hyperfine and quadrupolar interactions with 14 N, 45 Sc, and 89 Y nuclear spins are obtained. The systematic increase of the metal contribution to the triplet spin density from Y 3 N to Sc 3 N found in the ENDOR study points to a substantial fullerene-to-metal charge transfer in the excited state. These experimental results are rationalized with the help of ground-state and time-dependent DFT calculations, which revealed a substantial variation of the endohedral cluster position in the photoexcited states driven by the predisposition of Sc atoms to maximize their spin population.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Interrupted DNA and Slow Silver Cluster Luminescence

A DNA–silver cluster conjugate is a hierarchical chromophore with a partly reduced silver core embedded within the DNA nucleobases that are covalently linked by the phosphodiester backbone. Specific sites within a polymeric DNA can be targeted to spectrally tune the silver cluster. Here, the repeated (C 2 A) 6 strand is interrupted with a thymine, and the resulting (C 2 A) 2 -T-(C 2 A) 4 forms only Ag 10 6+ , a chromophore with both prompt (~1 ns) green and sustained (~10 2 μs) red luminescence. Thymine is an inert placeholder that can be removed, and the two fragments (C 2 A) 2 and (C 2 A) 4 also produce the same Ag 10 6+ adduct. In relation to (C 2 A) 2 T(C 2 A) 4 , the (C 2 A) 2 + (C 2 A) 4 pair is distinguished because the red Ag 10 6+ luminescence is ~6× lower, relaxes ~30% faster, and is quenched ~2× faster with O 2 . These differences suggest that a specific break in the phosphodiester backbone can regulate how a contiguous vs broken scaffold wraps and better protects its cluster adduct.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Lagrangian statistics of a shock-driven turbulent dynamo in decaying turbulence

Small-scale fluctuating magnetic fields of order nG are observed in supernova shocks and galaxy clusters, where its amplification is likely caused by the Biermann battery mechanism. However, these fields cannot be amplified further without the turbulent dynamo, which generates magnetic energy through the stretch-twist-fold (STF) mechanism. Thus, we present here novel 3D magnetohydrodynamic (MHD) simulations of a laser-driven shock propagating into a stratified, multiphase medium, to investigate the post-shock turbulent magnetic field amplification via the turbulent dynamo. The configuration used here is currently being tested in the shock tunnel at the National Ignition Facility (NIF). In order to probe the statistical properties of the post-shock turbulent region, we use 384 × 512 × 384 tracers to track its evolution through the Lagrangian framework, thus providing a high-fidelity analysis of the shocked medium. Our simulations indicate that the growth of the magnetic field, which accompanies the near-Saffman kinetic energy decay (E kin ∝ t -1.15 ) without turbulence driving, exhibits slightly different characteristics as compared to periodic box simulations. Seemingly no distinct phases exist in its evolution because the shock passage and time to observe the magnetic field amplification during the turbulence decay are very short (~0.3 of a turbulent turnover time). Yet, the growth rate is still consistent with those expected for compressive (curl-free) turbulence driving in subsonic, compressible turbulence. Phenomenological understanding of the dynamics of the magnetic and velocity fields are also elucidated via Lagrangian frequency spectra, which are consistent with the expected inertial range scalings in the Eulerian–Lagrangian bridge.

79 ASTRONOMY AND ASTROPHYSICS↗

Predicting the oxidation states of Mn ions in the oxygen-evolving complex of photosystem II using supervised and unsupervised machine learning

Abstract Serial Femtosecond Crystallography at the X-ray Free Electron Laser (XFEL) sources enabled the imaging of the catalytic intermediates of the oxygen evolution reaction of Photosystem II (PSII). However, due to the incoherent transition of the S-states, the resolved structures are a convolution from different catalytic states. Here, we train Decision Tree Classifier and K-means clustering models on Mn compounds obtained from the Cambridge Crystallographic Database to predict the S-state of the X-ray, XFEL, and CryoEM structures by predicting the Mn’s oxidation states in the oxygen-evolving complex. The model agrees mostly with the XFEL structures in the dark S 1 state. However, significant discrepancies are observed for the excited XFEL states (S 2 , S 3, and S 0 ) and the dark states of the X-ray and CryoEM structures. Furthermore, there is a mismatch between the predicted S-states within the two monomers of the same dimer, mainly in the excited states. We validated our model against other metalloenzymes, the valence bond model and the Mn spin densities calculated using density functional theory for two of the mismatched predictions of PSII. The model suggests designing a more optimized sample delivery and illumiation systems are crucial to precisely resolve the geometry of the advanced S-states to overcome the noncoherent S-state transition. In addition, significant radiation damage is observed in X-ray and CryoEM structures, particularly at the dangler Mn center (Mn4). Our model represents a valuable tool for investigating the electronic structure of the catalytic metal cluster of PSII to understand the water splitting mechanism.

Plant Sciences↗