Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Point Optimization”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Gamma-ray burst detection prospects for next generation ground-based VHE facilities

ABSTRACT Gamma-ray Bursts (GRB) were discovered by satellite-based detectors as powerful sources of transient γ-ray emission. The Fermi satellite detected an increasing number of these events with its dedicated Gamma-ray Burst Monitor (GBM), some of which were associated with high energy photons ($E \gt 10$ GeV), by the Large Area Telescope (LAT). More recently, follow-up observations by Cherenkov telescopes detected very high energy emission ($E \gt 100$ GeV) from GRBs, opening up a new observational window with implications on the interpretation of their central engines and on the propagation of very energetic photons across the Universe. Here, we use the data published in the 2nd Fermi-LAT Gamma Ray Burst Catalogue to characterize the duration, luminosity, redshift, and light curve of the high energy GRB emission. We extrapolate these properties to the very high energy domain, comparing the results with available observations and with the potential of future instruments. We use observed and simulated GRB populations to estimate the chances of detection with wide-field ground-based γ-ray instruments. Our analysis aims to evaluate the opportunities of the Southern Wide-field-of-view Gamma-ray Observatory (SWGO), to be installed in the Southern Hemisphere, to complement CTA. We show that a low-energy observing threshold ($E_{low} \lt 200$ GeV), with good point source sensitivity ($F_{lim} \approx 10^{-11} \,\mathrm{erg\, cm^{-2}\, s^{-1}}$ in $1$ yr), are optimal requirements to work as a GRB trigger facility and to probe the burst spectral properties down to time-scales as short as $10$ s, accessing a time domain that will not be available to Imaging Atmospheric Cherenkov Telescopes instruments.

La Mura, G. (ORCID:000000018553499X)↗
DOE OSTI2026

Simulation of Divertor Performance in ST40 Under Dynamic Double-Null Plasmas

A power fraction model was implemented for the simultaneous prediction of 3-D surface temperature evolution at all four divertor targets in near-double-null (DN) tokamak configurations, which is especially important for compact high-field devices that may not have the ability to dissipate large amounts of power on the high-field side. Evaluating the power-sharing between the four divertor strike points in a disconnected DN configuration is important for understanding the overall power balance, as well as for optimizing the power exhaust performance and prolonging the survivability of the plasma facing components (PFCs). This power-sharing is typically evaluated in terms of the separation between the primary and secondary separatrices at the outboard midplane, $\textit {dR}_{\text {sep}}$. The Heat flux Engineering Analysis Toolkit (HEAT) is coupled with Brunner’s power fraction model to simulate the deposited heat flux and resultant temperature change on 3-D divertor targets in a dynamic DN (DDN) pulse operation in ST40, a high-field spherical tokamak. The simulation results showed that with DDN operation, the operation time has significantly increased compared with single-null geometry configurations.

ST40↗
DOE OSTI2023

Final Cooling with Thick Wedges for a Muon Collider

A muon collider is a particle accelerator that collides muons rather than protons or electrons. An important component of such a collider is cooling, in which the emittance of the muon beam is reduced to make it suitable for use in a collider. In the last stage of this process, called final 4D cooling, emittance in the transverse axes is reduced while emittance in the longitudinal axis is allowed to grow. In previous studies, Final 4-D cooling used absorbers within very high field solenoids to cool low-momentum muons. Simulations of the systems did not reach the desired cooling design goals. In this study, we develop and optimize a different conceptual design for the final 4D cooling channel, which is based on single thick wedges. We used G4Beamline to simulate the channel and Python to generate and analyze particle distributions. We characterized the effects of various parameters on the performance of this system, including the initial emittances and Twiss parameters, the length and angle of the wedges, the length of the drift channel, and the length, frequency, gradient, and phase of the RF cavities. We optimized the design parameters of the cooling channel and produced two conceptual designs (corresponding to two possible starting points for the input beam) which achieve transverse cooling in both x and y by a factor of ≈ 3.5. These channels achieve a lower transverse and longitudinal emittance than the best design previously published.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Evaluating Portable Parallelization Strategies for Heterogeneous Architectures in High Energy Physics

High-energy physics (HEP) experiments have developed millions of lines of code over decades that are optimized to run on traditional x86 CPU systems. However, we are seeing a rapidly increasing fraction of floating point computing power in leadership-class computing facilities and traditional data centers coming from new accelerator architectures, such as GPUs. HEP experiments are now faced with the untenable prospect of rewriting millions of lines of x86 CPU code, for the increasingly dominant architectures found in these computational accelerators. This task is made more challenging by the architecture-specific languages and APIs promoted by manufacturers such as NVIDIA, Intel and AMD. Producing multiple, architecture-specific implementations is not a viable scenario, given the available person power and code maintenance issues. The Portable Parallelization Strategies team of the HEP Center for Computational Excellence is investigating the use of Kokkos, SYCL, OpenMP, std::execution::parallel and alpaka as potential portability solutions that promise to execute on multiple architectures from the same source code, using representative use cases from major HEP experiments, including the DUNE experiment of the Long Baseline Neutrino Facility, and the ATLAS and CMS experiments of the Large Hadron Collider. This cross-cutting evaluation of portability solutions using real applications will help inform and guide the HEP community when choosing their software and hardware suites for the next generation of experimental frameworks. We present the outcomes of our studies, including performance metrics, porting challenges, API evaluations, and build system integration.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

The effect of filament peaking on electron beam weld quality

Electron beam welding uses high intensity beams that are generated by an electron emitting cathode under high vacuum conditions. The cathode, typically a ribbon or wire filament of refractory metal, must be heated to sufficiently high temperatures for electron emission to occur and must be carefully controlled to produce the desired beam current while at the same time optimizing filament lifetime. Traditional methods for doing this in commonly used triode gun assemblies consist of first finding the filament knee point in the filament current versus beam current relationship that produces the desired beam current for a given beam voltage. The filament current is then increased by some amount, typically 5–10%, per conventional wisdom to produce a peaked beam. This investigation studies and quantifies the beam and weld quality produced by underkneed and overpeaked filaments using electron beam diagnostics to measure the power density distribution of these beams. Underkneed filaments are shown to have power densities that drop off very quickly below the traditional knee, resulting in poor weld penetration. Kneed filaments do not reach the full beam intensity and also produce shallow welds. Increasing the filament current 10% above the traditional knee was shown to produce a circular Gaussian-like beam with peaked intensity optimized for welding and for filament lifetime. Further increases in the filament current that produce overpeaked beams do not significantly change the beam or weld quality and would only contribute to a reduced filament lifetime.

Electron Beam Welding↗
DOE OSTI2021

Train Like a (Var)Pro: Efficient Training of Neural Networks with Variable Projection

Deep neural networks (DNNs) have achieved state-of-the-art performance across a variety of traditional machine learning tasks, e.g., speech recognition, image classification, and segmentation. The ability of DNNs to efficiently approximate high-dimensional functions has also motivated their use in scientific applications, e.g., to solve partial differential equations and to generate surrogate models. In this paper, we consider the supervised training of DNNs, which arises in many of the above applications. We focus on the central problem of optimizing the weights of the given DNN such that it accurately approximates the relation between observed input and target data. Devising effective solvers for this optimization problem is notoriously challenging due to the large number of weights, nonconvexity, data sparsity, and nontrivial choice of hyperparameters. To solve the optimization problem more efficiently, we propose the use of variable projection (VarPro), a method originally designed for separable nonlinear least-squares problems. Our main contribution is the Gauss--Newton VarPro method (GNvpro) that extends the reach of the VarPro idea to nonquadratic objective functions, most notably cross-entropy loss functions arising in classification. These extensions make GNvpro applicable to all training problems that involve a DNN whose last layer is an affine mapping, which is common in many state-of-the-art architectures. In our four numerical experiments from surrogate modeling, segmentation, and classification, GNvpro solves the optimization problem more efficiently than commonly used stochastic gradient descent (SGD) schemes. Finally, GNvpro finds solutions that generalize well, and in all but one example better than well-tuned SGD methods, to unseen data points.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Adaptive sampling for accelerating neutron diffraction-based strain mapping *

Abstract Neutron diffraction is a useful technique for mapping residual strains in dense metal objects. The technique works by placing an object in the path of a neutron beam, measuring the diffracted signals and inferring the local lattice strain values from the measurement. In order to map the strains across the entire object, the object is stepped one position at a time in the path of the neutron beam, typically in raster order, and at each position a strain value is estimated. Typical dwell times at neutron diffraction instruments result in an overall measurement that can take several hours to map an object that is several tens of centimeters in each dimension at a resolution of a few millimeters, during which the end users do not have an estimate of the global strain features and are at risk of incomplete information in case of instruments outages. In this paper, we propose an object adaptive sampling strategy to measure the significant points first. We start with a small initial uniform set of measurement points across the object to be mapped, compute the strain in those positions and use a machine learning technique to predict the next position to measure in the object. Specifically, we use a Bayesian optimization based on a Gaussian process regression method to infer the underlying strain field from a sparse set of measurements and predict the next most informative positions to measure based on estimates of the mean and variance in the strain fields estimated from the previously measured points. We demonstrate our real-time measure-infer-predict workflow on additively manufactured steel parts—demonstrating that we can get an accurate strain estimate even with 30%–40% of the typical number of measurements—leading the path to faster strain mapping with useful real-time feedback. We emphasize that the proposed method is general and can be used for fast mapping of other material properties such as phase fractions from time-consuming point-wise neutron measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Minimizing Optimal Transport for Functions with Fixed-Size Nodal Sets

Consider the class of zero-mean functions with fixed L ∞ and L 1 norms and exactly N ϵ N nodal points. Which functions f minimize W p (f + ,f – ), the Wasserstein distance between the measures whose densities are the positive and negative parts? We provide a complete solution to this minimization problem on the line and the circle, which provides sharp constants for previously proven “uncertainty principle”-type inequalities, i.e., lower bounds on N • W p (f + ,f – ). We further show that, while such inequalities hold in many metric measure spaces, they are no longer sharp when the non-branching assumption is violated; indeed, for metric star-graphs, the optimal lower bound on W p (f + ,f – ) is not inversely proportional to the size of the nodal set, N. Here, based on similar reductions, we make connections between the analogous problem of minimizing W p (f + ,f – ) for f defined on Ω C R d with an equivalent optimal domain partition problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Surpassing the Energy Resolution Limit with Ferromagnetic Torque Sensors

We discuss the fundamental noise limitations of a ferromagnetic torque sensor based on a levitated magnet in the tipping regime. We evaluate the optimal magnetic field resolution taking into account the thermomechanical noise and the mechanical detection noise at the standard quantum limit. We find that the energy resolution limit, pointed out in recent literature as a relevant benchmark for most classes of magnetometers, can be surpassed by many orders of magnitude. Moreover, similarly to the case of a ferromagnetic gyroscope, it is also possible to surpass the standard quantum limit for magnetometry with independent spins, arising from spin-projection noise. Our finding indicates that magnetomechanical systems optimized for magnetometry can achieve a magnetic field resolution per unit volume several orders of magnitude better than any conventional magnetometer. Here, we discuss possible implications, focusing on fundamental physics problems such as the search for exotic interactions beyond the standard model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

𝜂 and 𝜂′ Production in 𝐽/𝜓 Radiative Decays from Quantum Chromodynamics

We present a first principles calculation within quantum chromodynamics (QCD) of the radiative decays of the 𝐽/𝜓 into the light pseudoscalar mesons 𝜂 and 𝜂′. Within a lattice computation we obtain the transition form factors as a function of photon virtuality from the timelike region, accessible experimentally via the “Dalitz” decay 𝐽/𝜓 →𝑒 + ⁢𝑒 − ⁢𝜂 (′) , through to the real photon point corresponding to 𝐽/𝜓 →𝛾⁢𝜂 (′) . This is the first calculation in lattice QCD with two (heavier than physical) degenerate flavors of light quark and a heavier strange quark, in which the 𝜂′ appears as the first excited state with pseudoscalar isoscalar quantum numbers. We access it reliably by using variationally optimized operators , the use of which also improves the purity of the 𝐽/𝜓 and 𝜂 signals, reducing systematic uncertainties. High quality results at a large number of kinematic points are obtained in a typically noisy disconnected process by using a novel correlator averaging procedure. Our results show the expected enhanced production of the 𝜂′ over the 𝜂 in this process, and suggest that the demonstrated lattice technology is suitable for future calculations considering processes in which light meson resonances are produced.

form factors↗
DOE OSTI2020

LiDAR-Aided Interior Orientation Parameters Refinement Strategy for Consumer-Grade Cameras Onboard UAV Remote Sensing Systems

Unmanned aerial vehicles (UAVs) are quickly emerging as a popular platform for 3D reconstruction/modeling in various applications such as precision agriculture, coastal monitoring, and emergency management. For such applications, LiDAR and frame cameras are the two most commonly used sensors for 3D mapping of the object space. For example, point clouds for the area of interest can be directly derived from LiDAR sensors onboard UAVs equipped with integrated global navigation satellite systems and inertial navigation systems (GNSS/INS). Imagery-based mapping, on the other hand, is considered to be a cost-effective and practical option and is often conducted by generating point clouds and orthophotos using structure from motion (SfM) techniques. Mapping with photogrammetric approaches requires accurate camera interior orientation parameters (IOPs), especially when direct georeferencing is utilized. Most state-of-the-art approaches for determining/refining camera IOPs depend on ground control points (GCPs). However, establishing GCPs is expensive and labor-intensive, and more importantly, the distribution and number of GCPs are usually less than optimal to provide adequate control for determining and/or refining camera IOPs. Moreover, consumer-grade cameras with unstable IOPs have been widely used for mapping applications. Therefore, in such scenarios, where frequent camera calibration or IOP refinement is required, GCP-based approaches are impractical. To eliminate the need for GCPs, this study uses LiDAR data as a reference surface to perform in situ refinement of camera IOPs. The proposed refinement strategy is conducted in three main steps. An image-based sparse point cloud is first generated via a GNSS/INS-assisted SfM strategy. Then, LiDAR points corresponding to the resultant image-based sparse point cloud are identified through an iterative plane fitting approach and are referred to as LiDAR control points (LCPs). Finally, IOPs of the utilized camera are refined through a GNSS/INS-assisted bundle adjustment procedure using LCPs. Seven datasets over two study sites with a variety of geomorphic features are used to evaluate the performance of the developed strategy. The results illustrate the ability of the proposed approach to achieve an object space absolute accuracy of 3–5 cm (i.e., 5–10 times the ground sampling distance) at a 41 m flying height.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Optimization of Preparation Method, Nucleating Agent, and Stabilizers for Synthesizing Calcium Chloride Hexahydrate (CaCl2.6H2O) Phase Change Material

This study investigates improvements in low-cost latent heat storage material calcium chloride hexahydrate (CaCl2.6H2O). Its melting point is between 25 and 28 °C, with relatively high enthalpy (170–190 J/g); however, this phase change material (PCM) shows supercooling and phase separation. In CaCl2.6H2O incongruent melting causes lower hydrates of CaCl2 to form, which affects the overall energy storage capacity and long-term durability. In this work, PCM performance enhancement was achieved by adding SrCl2.6H2O as a nucleating agent and NaCl/KCl as a stabilizer to prevent supercooling and phase separation, respectively. We investigated the PCM preparation method and optimized the proportions of SrCl2.6H2O and NaCl/KCl. Thermal testing for 25 cycles combined with DSC and T-history testing was performed to observe changes in enthalpy, phase transitions and supercooling over the extended period of usage. X-ray diffraction was used to verify crystalline structure in the compounds. It was found that the addition of 2 wt.% of SrCl2.6H2O reduced supercooling from 12 °C to 0 °C compared to unmodified CaCl2.6H2O. The addition of 5 wt.% NaCl or KCl proved to effectively suppress separation and the melting enthalpy achieved was 169 J/g–178 J/g with congruent melting over 25 cycles, with no supercooling and almost no reduction in the latent heat.

42 ENGINEERING↗
DOE OSTI2025

Demonstration of an AI-driven workflow for dynamic x-ray spectroscopy

X-ray absorption near edge structure (XANES) spectroscopy is a powerful technique for characterizing the chemical state and symmetry of individual elements within materials, but requires collecting data at many energy points which can be time-consuming. While adaptive sampling methods exist for efficiently collecting spectroscopic data, they often lack domain-specific knowledge about the structure of XANES spectra. Here we demonstrate a knowledge-injected Bayesian optimization approach for adaptive XANES data collection that incorporates understanding of spectral features like absorption edges and pre-edge peaks. We show this method accurately reconstructs the absorption edge of XANES spectra using only 15–20% of the measurement points typically needed for conventional sampling, while maintaining the ability to determine the x-ray energy of the sharp peak after the absorption edge with errors less than 0.03 eV, the absorption edge with errors less than 0.1 eV; and overall root-mean-square errors less than 0.005 compared to traditionally sampled spectra. Our experiments on battery materials and catalysts demonstrate the method’s effectiveness for both static and dynamic XANES measurements, improving data collection efficiency and enabling better time resolution for tracking chemical changes. This approach advances the degree of automation in XANES experiments, reducing the common errors of under- or over-sampling points near the absorption edge and enabling dynamic experiments that require high temporal resolution or limited measurement time.

Bayesian optimization↗
DOE OSTI2025

ASCR Workshop Position Paper: Challenges and Opportunities in High Energy Physics

High energy particle physics and cosmology concern themselves with estimating fundamental parameters of nature, such as the masses and interactions of fundamental particles like the Higgs boson and the rate of expansion of the universe. In doing so, they analyze exabyte-scale datasets, some of the largest in all of science, and face many challenges in subsequent data analysis. These challenges are shared between the two disciplines, but we focus on particle physics to highlight one specific domain. In particle physics, the standard method for estimating parameters involves performing Monte Carlo (MC) integration as a function of both parameters of interest and nuisance parameters using an expensive simulator, counting the number of observed collision events (i.i.d. samples) from an experiment in the corresponding integration domains, and forming a Poisson likelihood function. This likelihood function is then used in a Frequentist manner to construct a maximum likelihood point estimate (MLE) and confidence set for the parameters. To sufficiently populate the high-dimensional integration domains, simulators consume billions of CPU-hours annually and produce hundreds of petabytes of intermediate output data. Several techniques have been developed to: optimize definitions of the integration domains so as to be maximally sensitive to a particular subset of parameters, efficiently estimate the integrals, and build robust surrogate models by interpolating between integral evaluations at different parameter points. One can view this whole endeavor as classical Simulation-Based Inference (SBI).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Optimal Hairpin Winding Configuration for EV Traction Motors to Enhance Active CMV Cancellation Capability of NPL.X Inverter

Active common-mode voltage (CMV) cancellation in neutral-point-less (NPL.X) inverter is an effective approach for reducing common-mode electromagnetic interference and alleviating voltage stress associated with high-speed switching in EV traction drives. However, the existing hairpin winding configuration often exhibits impedance imbalance due to slot leakage inductance and manufacturing constraints, limiting the effectiveness of active CMV cancellation. This paper investigates symmetric hairpin winding configurations for a dual-three-phase EV traction machine to improve impedance symmetry. Considering the manufacturing constraints of six-layer hairpin windings, four feasible symmetric winding configurations are analyzed. Two-dimensional finite-element analysis is performed to extract the impedance parameters of each configuration. The extracted impedance parameters are used to quantify the impedance mismatch between complementary phase pairs and evaluate its influence on active CMV cancellation. The results show that the fully interleaved winding configuration achieves perfectly balanced phase impedances with 0 % mismatch, whereas the mismatch in the other configurations ranges from approximately 2.4 % to 0.8 %. Consequently, the resulting total CMV overshoot is significantly suppressed to near zero voltage under high-voltage and high switching frequency operating conditions, demonstrating the importance of winding symmetry for maximizing the active CMV cancellation capability of the NPL.X inverter.

Lee, Kangbeen [Purdue Univ., West Lafayette, IN (U↗
DOE OSTI2025

The role of moisture in MgCl 2 salt: A multiscale approach to TES performance

This study presents a comprehensive multiscale analysis to evaluate the influence of moisture on the thermal performance of thermal energy storage systems using magnesium chloride (MgCl 2 ) as the phase-changing material. The system uses graphite foam with 90% relative density to enhance thermal conductivity. The analysis includes thermal conductivity calculations and specific heat capacity for different hydrate phases of MgCl 2 : Anhydrous, Mono, Di, Tetra, and Hexa. These properties were evaluated for the first time using the phonon density of states from Density Functional Theory simulations. Results showed that thermal conductivity decreased, while specific heat capacity increased by a factor of two as the phase changed from anhydrous to hexahydrate. Meso-scale models were created to account for the anisotropy of graphite foam and property variations of the MgCl 2 hydrate phases. Asymptotic Expansion homogenization simulations determined the anisotropic thermal conductivity for all phases. This unique methodology improved simulation accuracy, which matched experimental data for anhydrous MgCl 2 . A parametric study examined various operating conditions and their effect on TES performance. It revealed that higher charging temperatures did not enhance exergy efficiency, but increased discharge mass flow rates improved it due to better heat transfer. Thermal performance evaluated by the exergy efficiency remained consistent across hydrate systems under the tested conditions. Furthermore, the study suggests that this uniformity is linked to missing key data, particularly latent heat and melting point for different hydrates. Overall, it highlights the importance of multiscale effects and accurate material properties in designing and optimizing TES systems.

Molten salt degradation↗
DOE OSTI2020

PyFLOSIC: Python-based Fermi–Löwdin orbital self-interaction correction

We present pyflosic, an open-source, general-purpose python implementation of the Fermi–Löwdin orbital self-interaction correction (FLO-SIC), which is based on the python simulation of chemistry framework (pyscf) electronic structure and quantum chemistry code. Thanks to pyscf, pyflosic can be used with any kind of Gaussian-type basis set, various kinds of radial and angular quadrature grids, and all exchange-correlation functionals within the local density approximation, generalized-gradient approximation (GGA), and meta-GGA provided in the libxc and xcfun libraries. A central aspect of FLO-SIC is the Fermi-orbital descriptors, which are used to estimate the self-interaction correction. Importantly, they can be initialized automatically within pyflosic; they can also be optimized within pyflosic with an interface to the atomic simulation environment, a python library that provides a variety of powerful gradient-based algorithms for geometry optimization. Although pyflosic has already facilitated applications of FLO-SIC to chemical studies, it offers an excellent starting point for further developments in FLO-SIC approaches, thanks to its use of a high-level programming language and pronounced modularity.

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

A nearsighted force-training approach to systematically generate training data for the machine learning of large atomic structures

A challenge of atomistic machine-learning (ML) methods is ensuring that the training data are suitable for the system being simulated, which is particularly challenging for systems with large numbers of atoms. Most atomistic ML approaches rely on the nearsightedness principle (“all chemistry is local”), using information about the position of an atom’s neighbors to predict a per-atom energy. Here, in this work, we develop a framework that exploits the nearsighted nature of ML models to systematically produce an appropriate training set for large structures. We use a per-atom uncertainty estimate to identify the most uncertain atoms and extract chunks centered around these atoms. It is crucial that these small chunks are both large enough to satisfy the ML’s nearsighted principle (that is, filling the cutoff radius) and are large enough to be converged with respect to the electronic structure calculation. We present data indicating when the electronic structure calculations are converged with respect to the structure size, which fundamentally limits the accuracy of any nearsighted ML calculator. These new atomic chunks are calculated in electronic structures, and crucially, only a single force—that of the central atom—is added to the growing training set, preventing the noisy and irrelevant information from the piece’s boundary from interfering with ML training. The resulting ML potentials are robust, despite requiring single-point calculations on only small reference structures and never seeing large training structures. We demonstrated our approach via structure optimization of a 260-atom structure and extended the approach to clusters with up to 1415 atoms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗