Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cluster computer”

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 1,009 records · Page 56

First-principles insights into interaction of CO, NO, and HCN with Ag{sub 8}

We use static as well as time-dependent first-principles computations to study interaction of the CO, NO, and HCN molecules with the Ag{sub 8} nanocluster. The many-body based GW correction is applied for accurate description of the highest occupied (HOMO) and the lowest unoccupied (LUMO) molecular orbital levels. It is argued that the adsorption of these molecules changes the stable structure of Ag{sub 8} from Td to the more chemically active D{sub 2d} symmetry. We discuss that the CO, NO, and HCN molecules prefer to adsorb on the atom of the cluster with significant contribution to both HOMO and LUMO, for the accomplishment of the required charge transfers in the systems. The charge back donation is found to leave an excess energy of about 110 meV on the NO molecular bond, evidencing potential application of silver clusters for NO reduction. It is argued that CO and specially NO exhibit strong physical interaction with the silver cluster and hence significantly modify the electronic and optical properties of the system, while HCN makes very week physical bonds with the cluster. The optical absorption spectra of the Ag{sub 8} cluster before and after molecule adsorption are computed and a nontrivial red shift is observed in the NO and HCN adsorbed clusters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DCA++: A software framework to solve correlated electron problems with modern quantum cluster methods

We present the first open release of the DCA++ project, a high-performance research software framework to solve quantum many-body problems with cutting edge quantum cluster algorithms. DCA++ implements the dynamical cluster approximation (DCA) and its DCA extension with a continuous self-energy. The algorithms capture nonlocal correlations in strongly correlated electron systems, thereby giving insight into high-$T_c$ superconductivity. The code’s scalability allows efficient usage of systems at all scales, from workstations to leadership computers. With regard to the increasing heterogeneity of modern computing machines, DCA++ provides portable performance on conventional and emerging new architectures, such as hybrid CPU–GPU, sustaining multiple petaflops on ORNL’s Titan and CSCS’ Piz Daint supercomputers. Moreover, we show how sustainable and scalable development of the code base has been achieved by adopting standard techniques of the software industry. These include employing a distributed version control system, applying test-driven development and following continuous integration.

97 MATHEMATICS AND COMPUTING↗

Angular-momentum projection in coupled-cluster theory: Structure of 34 Mg

Single- reference coupled-cluster theory is an accurate and affordable computational method for the nuclear many-body problem. For open-shell nuclei, the reference state typically breaks rotational invariance and angular momentum must be restored as a good quantum number. We perform angular-momentum projection after variation and employ the disentangled coupled-cluster formalism and a Hermitian approach. We compare our results with benchmarks for 8 Be and 20 Ne using a two-nucleon interaction from chiral effective field theory and for pf-shell nuclei within the traditional shell model. We compute the rotational band in the exotic nucleus 34 Mg and find agreement with data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DUST IN CLUSTERS: SEPARATING THE CONTRIBUTION OF GALAXIES AND INTRACLUSTER MEDIA

We have analyzed a sample of 327 clusters of galaxies, spanning the range of 0.06–0.70 in redshift. Strong constraints on their mean intracluster emission of dust have been obtained using maps and catalogs from the Herschel MerMES project; within a radius of 5 arcmin centered in each cluster, the 95% C.L. limits obtained are 86.6, 48.2, and 30.9 mJy at the observed frequencies of 250, 350, and 500 μ m. From these restrictions, and assuming physical parameters typical of interstellar media in the Milky Way, we have obtained tight upper limits on the visual extinction of background galaxies due to the intracluster media (ICM): A {sub V} (95% C.L.) ≲ 10{sup −3} mag. Strong constraints are also obtained for the mass of such dust; for instance, using the data at 350 μ m we establish a 95% upper limit of <10{sup 9} M {sub ⊙} within a circle with a radius of 5 arcmin centered in the clusters. This corresponds to a fraction of the total mass of the clusters of 9.5 × 10{sup −6}, and indicates a deficiency in the gas-to-dust ratio in the ICM by about three orders of magnitude in relation to the value found in the Milky Way. Computing the total infrared luminosity of the clusters in three ranges of redshift (0.05–0.24, 0.24–0.42, and 0.42–0.71) and two ranges of mass (<10{sup 14} and >10{sup 14} M {sub ⊙}), respectively, a strong evolution of luminosity in redshift ( L ∼ z {sup 1.5}) for both ranges of masses is found. The results indicate a strong declining in star formation rate with time in the last ∼6 Gyr.

79 ASTRONOMY AND ASTROPHYSICS↗

Reducing the Volume of NASA Earth-Science Data

A computer program reduces data generated by NASA Earth-science missions into representative clusters characterized by centroids and membership information, thereby reducing the large volume of data to a level more amenable to analysis. The program effects an autonomous data-reduction/clustering process to produce a representative distribution and joint relationships of the data, without assuming a specific type of distribution and relationship and without resorting to domain-specific knowledge about the data. The program implements a combination of a data-reduction algorithm known as the entropy-constrained vector quantization (ECVQ) and an optimization algorithm known as the differential evolution (DE). The combination of algorithms generates the Pareto front of clustering solutions that presents the compromise between the quality of the reduced data and the degree of reduction. Similar prior data-reduction computer programs utilize only a clustering algorithm, the parameters of which are tuned manually by users. In the present program, autonomous optimization of the parameters by means of the DE supplants the manual tuning of the parameters. Thus, the program determines the best set of clustering solutions without human intervention.

Lee, Seungwon↗

Coupled-cluster representation of Green function employing modified spectral resolutions of similarity transformed Hamiltonians

In this paper we discuss a new formalism for producing an analytic coupled-cluster (CC) Green's function for an N-electron system by shifting the poles of similarity transformed Hamiltonians represented in N − 1 and N + 1 electron Hilbert spaces. Simple criteria are derived for the states in N − 1 and N + 1 electron spaces that are then corrected in the spectral resolution of the corresponding matrix representations of the similarity transformed Hamiltonian. The accurate description of excited state processes within a Green's function formalism would be of significant importance to a number of scientific communities ranging from physics and chemistry to engineering and the biological sciences. This is because the Green's function methodology provides a direct path for not only calculating properties whose underlying origins come from coupled many-body interactions but also provides a straightforward path for calculating electron transport, response, and correlation functions that allows for a direct link with experiment. As a special case of this general formulation, we discuss the application of this technique for Green's function defined by the CC with singles and doubles representation of the ground-state wave function.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Artificial Neural Network Representation of Dynamic Criticality Excursion Experiment Data

This paper outlines time dependent considerations in applying the artificial neural network (ANN) paradigm to dynamic nuclear criticality excursion experiment data. The ANN amongst a variety of uses, can apply to function fitting of multiple input-multiple output data sets, pattern recognition, cluster analysis, and time series prediction. The ANN is a component of Machine Learning (ML) data analysis. The practical application of ML to criticality safety has been funded through the Y-12 Plant Directed Research and Demonstration (PDRD) program to support criticality safety and emergency planning efforts. The initial effort considered function fitting of criticality excursion data. Recent efforts successfully extended the ANN to criticality safety code validation. The initial ANN work applied to criticality excursions centered around multivariate function fitting of static (tabular) data taken from the French Consequences Radiologiques d'un Accident de Criticite (CRAC) and SILENE reactor campaigns. The former report (LeCorche) summarized important tabular data for the first 44 CRAC experiments. The initial ANN was trained to specific time-independent inputs, then applied to predict peak fission rate and total fissions for a wide variety of fissile solutions process conditions. This work extends that effort to apply the ANN paradigm to dynamic prediction of time-dependent criticality excursion experiment data. A dynamic criticality excursion in fissile solution represents a complex multiphysics phenomenon. Factors include among other things, coupled neutronic, multiphase fluid flow, temperature, pressure, acoustic shock, boiling heat transfer, reactivity feedback, and radiolytic gas considerations. There have been many efforts to model (over time) as coupled ordinary and (space) as partial differential equations (ODE, PDE). The ANN representation noted here does not rely on solving explicit ODE or PDEs. Rather the intent is to demonstrate that as a time series, the criticality excursion can be predicted by training specific ANN architectures to experimental input and output time series exclusively. The rationale for extending this work to the time-dependent regime is important in developing a 'criticality excursion emulator' solely on the basis of learning characteristics from experiment dynamic data. An emulator is different than a computer code simulator driven by the solution of ODEs and PDEs given select 'tuning parameters' that match outputs to experiments. A criticality excursion emulator, derived exclusively from an ANN or ML, could predict most likely excursion characteristics, trained from historic dynamic experiment data, presented with input values within the range of those experiments. The work noted shows that generic time-series ANN architectures can be satisfactorily applied to the prediction of dynamic, extremely chaotic criticality excursion power given an initial time dependent input series. In reality, an excursion emulator could be driven by user-supplied initial (t=0), time independent or constant values to generate time dependent outputs. Thus ongoing work includes developing a custom ANN such that initial dynamic conditions would be predicted. Those results then serve as inputs to a customized ANN architecture currently in development. (author)

42 ENGINEERING↗

A two-stage method for microcalcification cluster segmentation in mammography by deformable models

Purpose: Segmentation of microcalcification (MC) clusters in x-ray mammography is a difficult task for radiologists. Accurate segmentation is prerequisite for quantitative image analysis of MC clusters and subsequent feature extraction and classification in computer-aided diagnosis schemes. Methods: In this study, a two-stage semiautomated segmentation method of MC clusters is investigated. The first stage is targeted to accurate and time efficient segmentation of the majority of the particles of a MC cluster, by means of a level set method. The second stage is targeted to shape refinement of selected individual MCs, by means of an active contour model. Both methods are applied in the framework of a rich scale-space representation, provided by the wavelet transform at integer scales. Segmentation reliability of the proposed method in terms of inter and intraobserver agreements was evaluated in a case sample of 80 MC clusters originating from the digital database for screening mammography, corresponding to 4 morphology types (punctate: 22, fine linear branching: 16, pleomorphic: 18, and amorphous: 24) of MC clusters, assessing radiologists’ segmentations quantitatively by two distance metrics (Hausdorff distance—HDIST{sub cluster}, average of minimum distance—AMINDIST{sub cluster}) and the area overlap measure (AOM{sub cluster}). The effect of the proposed segmentation method on MC cluster characterization accuracy was evaluated in a case sample of 162 pleomorphic MC clusters (72 malignant and 90 benign). Ten MC cluster features, targeted to capture morphologic properties of individual MCs in a cluster (area, major length, perimeter, compactness, and spread), were extracted and a correlation-based feature selection method yielded a feature subset to feed in a support vector machine classifier. Classification performance of the MC cluster features was estimated by means of the area under receiver operating characteristic curve (Az ± Standard Error) utilizing tenfold cross-validation methodology. A previously developed B-spline active rays segmentation method was also considered for comparison purposes. Results: Interobserver and intraobserver segmentation agreements (median and [25%, 75%] quartile range) were substantial with respect to the distance metrics HDIST{sub cluster} (2.3 [1.8, 2.9] and 2.5 [2.1, 3.2] pixels) and AMINDIST{sub cluster} (0.8 [0.6, 1.0] and 1.0 [0.8, 1.2] pixels), while moderate with respect to AOM{sub cluster} (0.64 [0.55, 0.71] and 0.59 [0.52, 0.66]). The proposed segmentation method outperformed (0.80 ± 0.04) statistically significantly (Mann-Whitney U-test, p < 0.05) the B-spline active rays segmentation method (0.69 ± 0.04), suggesting the significance of the proposed semiautomated method. Conclusions: Results indicate a reliable semiautomated segmentation method for MC clusters offered by deformable models, which could be utilized in MC cluster quantitative image analysis.

60 APPLIED LIFE SCIENCES↗

Pair extended coupled cluster doubles

The accurate and efficient description of strongly correlated systems remains an important challenge for computational methods. Doubly occupied configuration interaction (DOCI), in which all electrons are paired and no correlations which break these pairs are permitted, can in many cases provide an accurate account of strong correlations, albeit at combinatorial computational cost. Recently, there has been significant interest in a method we refer to as pair coupled cluster doubles (pCCD), a variant of coupled cluster doubles in which the electrons are paired. This is simply because pCCD provides energies nearly identical to those of DOCI, but at mean-field computational cost (disregarding the cost of the two-electron integral transformation). Here, we introduce the more complete pair extended coupled cluster doubles (pECCD) approach which, like pCCD, has mean-field cost and reproduces DOCI energetically. We show that unlike pCCD, pECCD also reproduces the DOCI wave function with high accuracy. Moreover, pECCD yields sensible albeit inexact results even for attractive interactions where pCCD breaks down.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Tuning Catalytic Reactivity via Wetting Control through Oxygen Vacancies: Ru Clusters on Anatase TiO 2 and CeO 2 Supports

The shape of supported metal particles regulates their catalytic reactivity and is determined by the degree of wetting between the metal particle and the support surface. Flattened particles that wet support surfaces were reported in various catalytic systems, particularly in the subnanometer size regime. Such consequential metal–support wetting phenomena are poorly understood, and methods to study them on powder catalysts under realistic conditions are lacking. Here, we investigate the size-dependent wetting behaviors of Ru particles on two reducible-oxide supports, anatase TiO 2 (TiO 2 -A) and CeO 2 , under reducing catalytic conditions. X-ray absorption spectroscopy (XAS), low-energy ion scattering (LEIS), and density functional theory (DFT) are combined to determine the shape of Ru particles. Ru particles remain three-dimensional without wetting the TiO 2 -A support within the coverage range studied (0.06–0.98 Ru nm –2 ). In contrast, at low coverages (<0.25 Ru nm –2 ), Ru wets the CeO 2 support to form flat, disordered structures. The higher wettability of CeO 2 than TiO 2 -A is attributed to oxygen vacancies in the near-surface region. The shape difference between small Ru particles or clusters on the two supports leads to drastically contrasting catalytic reactivities in polyolefin hydrogenolysis, despite similar diameters. As a result, this work highlights the implications of metal–support wetting, or cluster shape, on catalytic behaviors of small metal clusters, while establishing the foundation for future systematic studies of such a phenomenon in realistic systems, by delivering a multitechnique methodology and revealing governing fundamental principles.

42 ENGINEERING↗

MEDUSA - An overset grid flow solver for network-based parallel computer systems

Continuing improvement in processing speed has made it feasible to solve the Reynolds-Averaged Navier-Stokes equations for simple three-dimensional flows on advanced workstations. Combining multiple workstations into a network-based heterogeneous parallel computer allows the application of programming principles learned on MIMD (Multiple Instruction Multiple Data) distributed memory parallel computers to the solution of larger problems. An overset-grid flow solution code has been developed which uses a cluster of workstations as a network-based parallel computer. Inter-process communication is provided by the Parallel Virtual Machine (PVM) software. Solution speed equivalent to one-third of a Cray-YMP processor has been achieved from a cluster of nine commonly used engineering workstation processors. Load imbalance and communication overhead are the principal impediments to parallel efficiency in this application.

Smith, Merritt H.↗

Status on the Development of a Modeling and Simulation Framework for the Economic Assessment of Nuclear Hybrid Energy Systems

An effort to design and build a modeling and simulation framework to assess the economic viability of Nuclear Hybrid Energy Systems (NHES) was undertaken in fiscal year 2015 (FY15). The purpose of this report is to document the various tasks associated with the development of such a framework and to provide a status on its progress. Several tasks have been accomplished. First, starting from a simulation strategy, a rigorous mathematical formulation has been achieved in which the economic optimization of a Nuclear Hybrid Energy System is presented as a constrained robust (under uncertainty) optimization problem. Some possible algorithms for the solution of the optimization problem are presented. A variation of the Simultaneous Perturbation Stochastic Approximation algorithm has been implemented in RAVEN and preliminary tests have been performed. The development of the software infrastructure to support the simulation of the whole NHES has also moved forward. The coupling between RAVEN and an implementation of the Modelica language (OpenModelica) has been implemented, migrated under several operating systems and tested using an adapted model of a desalination plant. In particular, this exercise was focused on testing the coupling of the different code systems; testing parallel, computationally expensive simulations on the INL cluster; and providing a proof of concept for the possibility of using surrogate models to represent the different NHES subsystems. Another important step was the porting of the RAVEN code under the Windows™ operating system. This accomplishment makes RAVEN compatible with the development environment that is being used for dynamic simulation of NHES components. A very simplified model of a NHES on the electric market has been built in RAVEN to confirm expectations on the analysis capability of RAVEN to provide insight into system economics and to test the capability of RAVEN to identify limit surfaces even for stochastic constraints. This capability will be needed in the future to enforce the stochastic constraints on the electric demand coverage from the NHES. The development team gained experience with many of the tools that are currently envisioned for use in the economic analysis of NHES and completed several important steps. Given the complexity of the project, preference has been given to a structural approach in which several independent efforts have been used to build the cornerstone of the simulation framework. While this is good approach in establishing such a complex framework, it may delay reaching more complete results on the performance of analyzed system configurations. The integration of the previously reported exergy analysis approach was initially proposed as part of this milestone. However, in reality, the exergy-based apportioning of cost will take place only in a second stage of the implementation since it will be used to properly allocate cost among the different NHES subsystems. Therefore, exergy does not appear at the level of the main drivers in the analysis framework; the latter development of the base framework is the focus of this report.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

The Physics and Chemistry of Cluster-based Catalyst Systems

In this project we used first principles calculations to investigate the physical and chemical properties of atomic clusters, small clumps of matter containing up to a few hundred atoms. The calculations are based on density functional theory, the most efficient and accurate method available for computing and thereby predicting the properties of molecules, clusters, and solids. An initial focus of our work involved studying the reactivity of small Pd clusters. We demonstrated that H2 molecules readily dissociate on Pd clusters and that multiple dissociations can occur on clusters larger than a few Pd atoms. For example, up to eight H2 molecules will dissociate on Pd13 at room temperature. A second focus of the project was to study the evolution of the electrical behavior of clusters, in particular, the evolution of metallic behavior in clusters of metal atoms. We developed a technique for using the electric dipole polarizability of the clusters – the response of the cluster to an applied electric field -- as a probe of this behavior. By decomposing the polarizability into so-called charge transfer (or metallic) and local (or insulating) contributions, we obtained the remarkable result that the per atom charge transfer contributions in sodium clusters attain the same value as occurs in bulk sodium metal for clusters containing as few as 10 atoms. Our calculated total polarizabilities matched experimentally measured trends essentially perfectly, validating the reliability of the calculations and supporting our conclusion that tiny clusters already display bulk-like metallic character. Similar results were obtained for other metal atom clusters, including potassium, aluminum, and zinc. We discovered that the evolution toward bulk metallic behavior for each of these cluster species is essentially identical, when judged using these polarizability-based descriptors. This suggests a universality in the approach of metal clusters to bulk-like metallicity. This discovery is a major outcome of our project. We also used our polarizability-based metrics to investigate the evolution of the properties of Si and NaCl clusters. These systems are non-metallic in the bulk limit. For Si, we find the behavior of small clusters to be metal-like and we attribute this to bonding arrangements in the clusters that are very different from the open structure of crystalline silicon. As the clusters get larger, the bonding in the cluster interior becomes more and more like that of the bulk material. We expect a transition from metal-like to insulator-like behavior to accompany this change is structure. Our results point to a transition at sizes beyond 146 atoms, the largest size treated in our study. We hope to study this transition further in the future. A third thrust of our project involved work on a new approach to eliminating unphysical electron self-interaction from density functional theory calculations. Two electrons exert repulsive forces on each other, but a single electron should not interact with itself. Because of approximations made in density functional theory, individual electrons effectively interact with themselves. The effect of this self-interaction is often small, especially when atoms are near their equilibrium positions in molecules. But the effects can be pronounced when bonds are stretched and also in some situations involving transition metals even when the atoms are near equilibrium. Our research involved the further development of a method for removing self-interaction called the Fermi-Löwdin orbital self-interaction correction or FLO-SIC. We conducted the first systematic FLO-SIC calculations for the individual atoms ranging from Li – Kr. We also carried out a critical comparison of the results of FLO-SIC with an older method known as the Perdew-Zunger self-interaction correction (PZ-SIC). FLO-SIC can be thought of as a modified approach to PZ-SIC that is expected to be computationally more efficient. We demonstrated that the two methods are distinct, but give essentially the same results for quantities such as molecular binding energies or electron removal energies. We are continuing to develop the FLO-SIC methodology with a goal of creating a computationally efficient method that is free of self-interaction effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Accelerating Simulation for High-Fidelity PV Inverter System Reliability Assessment with High-Performance Computing

The overall cost of photovoltaic (PV) systems has shown a downward trend during the last decade; however, PV inverter failures account for the highest cost of operation and maintenance. To address this, reliability tools with powerful computation and better accuracy are required for the lifetime prediction and degradation evaluation of PV inverters. This paper proposes an event-driven parallel computing-based simulator. The proposed simulator applies high-performance computing techniques and other accessory optimization techniques-including cluster merging, adaptive model updates, and steady-state identification-to make reliability assessments for PV inverters under given input mission profiles and operating conditions with high efficiency and high fidelity. The main idea of the simulator and its workflow are introduced. Then, a demo PV inverter system simulator is implemented, and the speedup of the total simulations of the switching model reaches 123.03 times.

high-performance computing↗

Parallel Navier-Stokes computations on shared and distributed memory architectures

We study a high order finite difference scheme to solve the time accurate flow field of a jet using the compressible Navier-Stokes equations. As part of our ongoing efforts, we have implemented our numerical model on three parallel computing platforms to study the computational, communication, and scalability characteristics. The platforms chosen for this study are a cluster of workstations connected through fast networks (the LACE experimental testbed at NASA Lewis), a shared memory multiprocessor (the Cray YMP), and a distributed memory multiprocessor (the IBM SPI). Our focus in this study is on the LACE testbed. We present some results for the Cray YMP and the IBM SP1 mainly for comparison purposes. On the LACE testbed, we study: (1) the communication characteristics of Ethernet, FDDI, and the ALLNODE networks and (2) the overheads induced by the PVM message passing library used for parallelizing the application. We demonstrate that clustering of workstations is effective and has the potential to be computationally competitive with supercomputers at a fraction of the cost.

Hayder, M. Ehtesham↗

High-throughput ab initio reaction mechanism exploration in the cloud with automated multi-reference validation

Quantum chemical calculations on atomistic systems have evolved into a standard approach to studying molecular matter. These calculations often involve a significant amount of manual input and expertise, although most of this effort could be automated, which would alleviate the need for expertise in software and hardware accessibility. Here, we present the AutoRXN workflow, an automated workflow for exploratory high-throughput electronic structure calculations of molecular systems, in which (i) density functional theory methods are exploited to deliver minimum and transition-state structures and corresponding energies and properties, (ii) coupled cluster calculations are then launched for optimized structures to provide more accurate energy and property estimates, and (iii) multi-reference diagnostics are evaluated to back check the coupled cluster results and subject them to automated multi-configurational calculations for potential multi-configurational cases. All calculations are carried out in a cloud environment and support massive computational campaigns. Key features of all components of the AutoRXN workflow are autonomy, stability, and minimum operator interference. We highlight the AutoRXN workflow with the example of an autonomous reaction mechanism exploration of the mode of action of a homogeneous catalyst for the asymmetric reduction of ketones.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Assessment of low-scaling approximations to the equation of motion coupled-cluster singles and doubles equations

Methods for fast and reliable computation of electronic excitation energies are in short supply, and little is known about their systematic performance. This work reports a comparison of several low-scaling approximations to the equation of motion coupled cluster singles and doubles (EOM–CCSD) and linear-response coupled cluster singles and doubles (LR–CCSD) equations with other single reference methods for computing the vertical electronic transition energies of 11 small organic molecules. The methods, including second order equation-of-motion many-body perturbation theory (EOM–MBPT2) and its partitioned variant, are compared to several valence and Rydberg singlet states. We find that the EOM–MBPT2 method was rarely more than a tenth of an eV from EOM–CCSD calculated energies, yet demonstrates a performance gain of nearly 30%. The partitioned equation-of-motion approach, P–EOM–MBPT2, which is an order of magnitude faster than EOM–CCSD, outperforms the CIS(D) and CC2 in the description of Rydberg states. CC2, on the other hand, excels at describing valence states where P–EOM–MBPT2 does not. The difference between the CC2 and P–EOM–MBPT2 can ultimately be traced back to how each method approximates EOM–CCSD and LR–CCSD. The results suggest that CC2 and P–EOM–MBPT2 are complementary: CC2 is best suited for the description of valence states while P–EOM–MBPT2 proves to be a superior O(N{sup 5}) method for the description of Rydberg states.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗