Showing papers on "Coupled cluster published in 2007"

Coupled-cluster theory in quantum chemistry

Abstract: Today, coupled-cluster theory offers the most accurate results among the practical ab initio electronic-structure theories applicable to moderate-sized molecules. Though it was originally proposed for problems in physics, it has seen its greatest development in chemistry, enabling an extensive range of applications to molecular structure, excited states, properties, and all kinds of spectroscopy. In this review, the essential aspects of the theory are explained and illustrated with informative numerical results.

Energy-consistent relativistic pseudopotentials and correlation consistent basis sets for the 4d elements Y-Pd.

Abstract: Scalar-relativistic pseudopotentials and corresponding spin-orbit potentials of the energy-consistent variety have been adjusted for the simulation of the [Ar]3d(10) cores of the 4d transition metal elements Y-Pd. These potentials have been determined in a one-step procedure using numerical two-component calculations so as to reproduce atomic valence spectra from four-component all-electron calculations. The latter have been performed at the multi-configuration Dirac-Hartree-Fock level, using the Dirac-Coulomb Hamiltonian and perturbatively including the Breit interaction. The derived pseudopotentials reproduce the all-electron reference data with an average accuracy of 0.03 eV for configurational averages over nonrelativistic orbital configurations and 0.1 eV for individual relativistic states. Basis sets following a correlation consistent prescription have also been developed to accompany the new pseudopotentials. These range in size from cc-pVDZ-PP to cc-pV5Z-PP and also include sets for 4s4p correlation (cc-pwCVDZ-PP through cc-pwCV5Z-PP), as well as those with extra diffuse functions (aug-cc-pVDZ-PP, etc.). In order to accurately assess the impact of the pseudopotential approximation, all-electron basis sets of triple-zeta quality have also been developed using the Douglas-Kroll-Hess Hamiltonian (cc-pVTZ-DK, cc-pwCVTZ-DK, and aug-cc-pVTZ-DK). Benchmark calculations of atomic ionization potentials and 4d(m-2)5s(2)-->4d(m-1)5s(1) electronic excitation energies are reported at the coupled cluster level of theory with extrapolations to the complete basis set limit.

Understanding of assembly phenomena by aromatic-aromatic interactions: benzene dimer and the substituted systems.

Abstract: Interactions involving aromatic rings are important in molecular/biomolecular assembly and engineering. As a consequence, there have been a number of investigations on dimers involving benzene or other substituted π systems. In this Feature Article, we examine the relevance of the magnitudes of their attractive and repulsive interaction energy components in governing the geometries of several π−π systems. The geometries and the associated binding energies were evaluated at the complete basis set (CBS) limit of coupled cluster theory with singles, doubles, and perturbative triples excitations [CCSD(T)] using a least biased scheme for the given data set. The results for the benzene dimer indicate that the floppy T-shaped structure (center-to-center distance: 4.96 A, with an axial benzene off-centered above the facial benzene) is isoenergetic in zero-point-energy (ZPE) corrected binding energy (D0) to the displaced-stacked structure (vertical interplanar distance: 3.54 A). However, the T-shaped structure i...

PSI3: An open‐source Ab Initio electronic structure package

Abstract: PSI3 is a program system and development platform for ab initio molecular electronic structure computations. The package includes mature programming interfaces for parsing user input, accessing commonly used data such as basis-set information or molecular orbital coefficients, and retrieving and storing binary data (with no software limitations on file sizes or file-system-sizes), especially multi-index quantities such as electron repulsion integrals. This platform is useful for the rapid implementation of both standard quantum chemical methods, as well as the development of new models. Features that have already been implemented include Hartree-Fock, multiconfigurational self-consistent-field, second-order Moller-Plesset perturbation theory, coupled cluster, and configuration interaction wave functions. Distinctive capabilities include the ability to employ Gaussian basis functions with arbitrary angular momentum levels; linear R12 second-order perturbation theory; coupled cluster frequency-dependent response properties, including dipole polarizabilities and optical rotation; and diagonal Born-Oppenheimer corrections with correlated wave functions. This article describes the programming infrastructure and main features of the package. PSI3 is available free of charge through the open-source, GNU General Public License.

Coupling term derivation and general implementation of state-specific multireference coupled cluster theories

Abstract: Simple closed-form expressions are derived for the “same vacuum” renormalization terms that arise in state-specific multireference coupled cluster (MRCC) theories. Explicit equations are provided for these coupling terms through the triple excitation level of MRCC theory, and a general expression is included for arbitrary-order excitations. The first production-level code (PSIMRCC) for state-specific and rigorously size-extensive Mukherjee multireference coupled cluster singles and doubles (MkCCSD) computations has been written. This code is also capable of evaluating analogous Brillouin-Wigner multireference energies (BWCCSD), including a posteriori size-extensivity corrections. Using correlation-consistent basis sets (cc-pVXZ, X=D,T,Q), MkCCSD and BWCCSD were tested and compared on two classic multireference problems: (1) the dissociation potential curve of molecular fluorine (F2) and (2) the structure and vibrational frequencies of ozone. Comparison with experimental data shows that the Mukherjee metho...

The current state of ab initio calculations of optical rotation and electronic circular dichroism spectra.

Abstract: The current ability of ab initio models to compute chiroptical properties such as optical rotatory dispersion and electronic circular dichroism spectra is reviewed. Comparison between coupled cluster linear response theory and experimental data (both gas and liquid phase) yields encouraging results for small to medium-sized chiral molecules including rigid species such as (S)-2-chloropropionitrile and (P)-[4]triangulane, as well as conformationally flexible molecules such as (R)-epichlorohydrin. More problematic comparisons are offered by (S)-methyloxirane, (S)-methylthiirane, and (1S,4S)-norbornenone, for which the comparison between theory and experiment is much poorer. The impact of basis-set incompleteness, electron correlation, zero-point vibration, and temperature are discussed. In addition, future prospects and obstacles for the development of efficient and reliable quantum chemical models of optical activity are discussed, including the problem of gauge invariance, scaling of the coupled cluster approach with system size, and solvation.

Dyson orbitals for ionization from the ground and electronically excited states within equation-of-motion coupled-cluster formalism: Theory, implementation, and examples

Abstract: Implementation of Dyson orbitals for coupled-cluster and equation-of-motion coupled-cluster wave functions with single and double substitutions is described and demonstrated by examples. Both ionizations from the ground and electronically excited states are considered. Dyson orbitals are necessary for calculating electronic factors of angular distributions of photoelectrons, Compton profiles, electron momentum spectra, etc, and can be interpreted as states of the leaving electron. Formally, Dyson orbitals represent the overlap between an initial N-electron wave function and the N1 electron wave function of the corresponding ionized system. For the ground state ionization, Dyson orbitals are often similar to the corresponding Hartree-Fock molecular orbitals MOs; however, for ionization from electronically excited states Dyson orbitals include contributions from several MOs and their shapes are more complex. The theory is applied to calculating the Dyson orbitals for ionization of formaldehyde from the ground and electronically excited states. Partial-wave analysis is employed to compute the probabilities to find the ejected electron in different angular momentum states using the freestanding and Coulomb wave representations of the ionized electron. Rydberg states are shown to yield higher angular momentum electrons, as compared to valence states of the same symmetry. Likewise, faster photoelectrons are most likely to have higher angular momentum. © 2007 American Institute of Physics. DOI: 10.1063/1.2805393

Quintuple-ζ quality coupled-cluster correlation energies with triple-ζ basis sets

Abstract: The explicitly-correlated coupled-cluster method CCSD(T)(R12) is extended to include F12 geminal basis functions that decay exponentially with the interelectronic distance and reproduce the form of the average Coulomb hole more accurately than linear-r12. Equations derived using the Ansatz 2 strong orthogonality projector are presented. The convergence of the correlation energy with orbital basis set for the new CCSD(T)(F12) method is studied and found to be rapid, 98% of the basis set limit correlation energy is typically recovered using triple-ζ orbital basis sets. The performance for reaction enthalpies is assessed via a test set of 15 reactions involving 23 molecules. The title statement is found to hold equally true for total and relative correlation energies.

Vibrational structure theory: new vibrational wave function methods for calculation of anharmonic vibrational energies and vibrational contributions to molecular properties

Abstract: A number of recently developed theoretical methods for the calculation of vibrational energies and wave functions are reviewed. Methods for constructing the appropriate quantum mechanical Hamilton operator are briefly described before reviewing a particular branch of theoretical methods for solving the nuclear Schrodinger equation. The main focus is on wave function methods using the vibrational self-consistent field (VSCF) as starting point, and includes vibrational configuration interaction (VCI), vibrational Moller–Plesset (VMP) theory, and vibrational coupled cluster (VCC) theory. The convergence of the different methods towards the full vibrational configuration interaction (FVCI) result is discussed. Finally, newly developed vibrational response methods for calculation of vibrational contributions to properties, energies, and transition probabilities are discussed.

Extension of the Renormalized Coupled-Cluster Methods Exploiting Left Eigenstates of the Similarity-Transformed Hamiltonian to Open-Shell Systems: A Benchmark Study†

Abstract: The recently formulated completely renormalized coupled-cluster method with singles, doubles, and noniterative triples, exploiting the biorthogonal form of the method of moments of coupled-cluster equations (Piecuch, P.; Wloch, M. J. Chem. Phys. 2005, 123, 224105; Piecuch, P.; Wloch, M.; Gour, J. R.; Kinal, A. Chem. Phys. Lett. 2006, 418, 467), termed CR-CC(2,3), is extended to open-shell systems. Test calculations for bond breaking in the OH radical and the ion and singlet−triplet gaps in the CH2, HHeH, and (HFH)- biradical systems indicate that the CR-CC(2,3) approach employing the restricted open-shell Hartree−Fock (ROHF) reference is significantly more accurate than the widely used CCSD(T) method and other noniterative triples coupled-cluster approximations without making the calculations substantially more expensive. A few molecular examples, including the activation energies of the C2H4 + H → C2H5 forward and reverse reactions and the triplet states of the CH2 and H2Si2O2 biradicals, are used to sho...

Fully automated implementation of the incremental scheme: application to CCSD energies for hydrocarbons and transition metal compounds.

Abstract: A general fully automated implementation of the incremental scheme for molecules and embedded clusters in the framework of the coupled cluster singles and doubles theory is presented. The code can be applied to arbitrary order of the incremental expansion and is parallelized in a master/slave structure. The authors found that the error in the total correlation energy is lower than 1kcal∕mol with respect to the canonical CCSD calculation if the incremental series is truncated in a proper way.

Coupled-cluster theory for three-body Hamiltonians

Abstract: We derive coupled-cluster equations for three-body Hamiltonians. The equations for the one- and two-body cluster amplitudes are presented in a factorized form that leads to an efficient numerical implementation. We employ low-momentum two- and three-nucleon interactions and calculate the binding energy of {sup 4}He. The results show that the main contribution of the three-nucleon interaction stems from its density-dependent zero-, one-, and two-body terms that result from the normal ordering of the Hamiltonian in coupled-cluster theory. The residual three-body terms that remain after normal ordering can be neglected.

The origin of deficiency of the supermolecule second-order Moller-Plesset approach for evaluating interaction energies.

Abstract: Calculations for the complex of thymine and adenine are used to show that the supermolecule second-order Moller-Plesset perturbation theory (MP2) approach for evaluating interaction energies fails in certain cases because of the behavior of one of its components: the uncoupled Hartree-Fock dispersion energy. A simple approach for correcting the MP2 supermolecule interaction energies is proposed. It focuses on correcting a relatively small difference between the MP2 and coupled cluster interaction energies, which is a very appealing feature of the new approach considering a benchmark role played by coupled cluster results.

Basis set convergence of post-CCSD contributions to molecular atomization energies

Abstract: Basis set convergence of correlation effects on molecular atomization energies beyond the coupled cluster with singles and doubles (CCSD) approximation has been studied near the one-particle basis set limit. Quasiperturbative connected triple excitations, (T), converge more rapidly than L−3 (where L is the highest angular momentum represented in the basis set), while higher-order connected triples, T3−(T), converge more slowly—empirically, ∝L−5∕2. Quasiperturbative connected quadruple excitations, (Q), converge smoothly as ∝L−3 starting with the cc-pVTZ basis set, while the cc-pVDZ basis set causes overshooting of the contribution in highly polar systems. Higher-order connected quadruples display only weak, but somewhat erratic, basis set dependence. Connected quintuple excitations converge very rapidly with the basis set, to the point where even an unpolarized double-zeta basis set yields useful numbers. In cases where fully iterative coupled cluster up to connected quintuples (CCSDTQ5) calculations are ...

A Novel Approach to Parallel Coupled Cluster Calculations: Combining Distributed and Shared Memory Techniques for Modern Cluster Based Systems

Abstract: A parallel coupled cluster algorithm that combines distributed and shared memory techniques for the CCSD(T) method (singles + doubles with perturbative triples) is described The implementation of the massively parallel CCSD(T) algorithm uses a hybrid molecular and “direct” atomic integral driven approach Shared memory is used to minimize redundant replicated storage per compute process The algorithm is targeted at modern cluster based architectures that are comprised of multiprocessor nodes connected by a dedicated communication network Parallelism is achieved on two levels: parallelism within a compute node via shared memory parallel techniques and parallelism between nodes using distributed memory techniques The new parallel implementation is designed to allow for the routine evaluation of mid- (500−750 basis function) to large-scale (750−1000 basis function) CCSD(T) energies Sample calculations are performed on five low-lying isomers of water hexamer using the aug-cc-pVTZ basis set

Reliable predictions of the thermochemistry of boron-nitrogen hydrogen storage compounds: BxNxHy, x = 2, 3.

Abstract: Thermochemical data calculated using ab initio molecular orbital theory are reported for 16 BxNxHy compounds with x = 2, 3 and y ≥ 2x. Accurate gas-phase heats of formation were obtained using coupled cluster with single and double excitations and perturbative triples (CCSD(T)) valence electron calculations extrapolated to the complete basis set (CBS) limit with additional corrections including core/valence, scalar relativistic, and spin−orbit corrections to predict the atomization energies and scaled harmonic frequencies to correct for zero point and thermal energies and estimate entropies. Computationally cheaper calculations were also performed using the G3MP2 and G3B3 variants of the Gaussian 03 method, as well as density functional theory (DFT) using the B3LYP functional. The G3MP2 heats of formation are too positive by up to ∼6 kcal/mol as compared with CCSD(T)/CBS values. The more expensive G3B3 method predicts heats of formation that are too negative as compared with the CCSD(T)/CBS values by up t...

Electrostatic energy in the effective fragment potential method: Theory and application to benzene dimer

Abstract: Evaluation of the electrostatic energy within the effective fragment potential (EFP) method is presented. The performance of two variants of the distributed multipole analysis (DMA) together with two different models for estimating the charge penetration energies was studied using six homonuclear dimers. The importance of damping the higher order multipole terms, i.e. charge dipole, was also investigated. Damping corrections recover more than 70% of the charge penetration energy in all dimers, whereas higher order damping introduces only minor improvement. Electrostatic energies calculated by the numerical DMA are less accurate than those calculated by the analytic DMA. Analysis of bonding in the benzene dimer shows that EFP with inclusion of the electrostatic damping term performs very well compared to the high-level coupled cluster singles, doubles, and perturbative triples method. The largest error of 0.4 kcal/mol occurs for the sandwich dimer configuration. This error is about half the size of the corresponding error in second order perturbation theory. Thus, EFP in the current implementation is an accurate and computationally inexpensive method for calculating interaction energies in weakly bonded molecular complexes. © 2006 Wiley Periodicals, Inc. J Comput Chem 28: 276–291, 2007

A Near Linear-Scaling Smooth Local Coupled Cluster Algorithm for Electronic Structure

Abstract: We demonstrate near linear scaling of a new algorithm for computing smooth local coupled-cluster singles-doubles (LCCSD) correlation energies of quantum mechanical systems. The theory behind our approach has been described previously, [J. Subotnik and M. Head-Gordon, J. Chem. Phys. 123, 064108 (2005)], and requires appropriately multiplying standard iterative amplitude equations by a bump function, creating local amplitude equations (which are smooth according to the implicit function theorem). Here, we provide an example that this theory works in practice: we show that our algorithm leads to smooth potential energy surfaces and yields large computational savings. As an example, we apply our LCCSD approach to measure the post-MP2 correction to the energetic gap between two different alanine tetrapeptide conformations.

Quantum-Chemical Investigation of the Structures and Electronic Spectra of the Nucleic Acid Bases at the Coupled Cluster CC2 Level

Abstract: We study the ground-state structures and singlet- and triplet-excited states of the nucleic acid bases by applying the coupled cluster model CC2 in combination with a resolution-of-the-identity approximation for electron interaction integrals. Both basis set effects and the influence of dynamic electron correlation on the molecular structures are elucidated; the latter by comparing CC2 with Hartree-Fock and Moller-Plesset perturbation theory to second order. Furthermore, we investigate basis set and electron correlation effects on the vertical excitation energies and compare our highest-level results with experiment and other theoretical approaches. It is shown that small basis sets are insuffient for obtaining accurate results for excited states of these molecules and that the CC2 approach to dynamic electron correlation is a reliable and efficient tool for electronic structure calculations on medium-sized molecules.

How accurate is the density functional theory combined with symmetry-adapted perturbation theory approach for CH–π and π–π interactions? A comparison to supermolecular calculations for the acetylene–benzene dimer

Abstract: Five different orientations of the acetylene–benzene dimer including the T-shaped global minimum structure are used to assess the accuracy of the density functional theory combined with symmetry adapted perturbation theory (DFT-SAPT) approach in its density-fitting implementation (DF-DFT-SAPT) for the study of CH–π and π–π interactions. The results are compared with the outcome of counterpoise corrected supermolecular calculations employing second-order Moller–Plesset (MP2), spin-component scaled MP2 (SCS-MP2) and single and double excitation coupled cluster theory including perturbative triple excitations (CCSD(T)). For all considered orientations MP2 predicts much deeper potential energy curves with considerably shifted minima compared to CCSD(T) and DFT-SAPT. In spite of being an improvement over the results of MP2, SCS-MP2 tends to underestimate the well depth while DFT-SAPT, employing an asymptotically corrected hybrid exchange–correlation potential in conjunction with the adiabatic local density approximation for the exchange–correlation kernel, is found to be in excellent agreement with CCSD(T). Furthermore, DFT-SAPT provides a detailed understanding of the importance of the electrostatic, induction and dispersion contributions to the total interaction energy and their repulsive exchange corrections.

Transition strengths and first-order properties of excited states from local coupled cluster CC2 response theory with density fitting.

Abstract: A new ab initio method for calculating transition strengths and orbital-unrelaxed first-order properties of singlet ground and excited states of extended molecular systems is presented. It is based on coupled cluster response theory at the level of the CC2 model with local approximations introduced to the doubles-excitation part of the wave function. Density fitting is employed for the calculation of the electron repulsion integrals, so that—with the exception of doubles amplitudes—only three-indexed objects do occur in the formalism. The new method was tested by performing calculations for a set of various molecules and excited states and by comparing the results with corresponding canonical (nonlocal) calculations. It turned out that for calculating transition strengths and properties of excited states the ordinary Boughton-Pulay domains are insufficient in numerous cases. To circumvent this problem a new scheme for extending domains is proposed, which is based on the solution of the coupled perturbed localization and Hartree-Fock equations. When such extended domains are used, a satisfactory agreement between canonical and local results is achieved.

Benchmark full configuration interaction and equation-of-motion coupled-cluster model with single and double substitutions for ionized systems results for prototypical charge transfer systems: noncovalent ionized dimers.

Abstract: Benchmark full configuration interaction and equation-of-motion coupled-cluster model with single and double substitutions for ionized systems (EOM-IP-CCSD) results are presented for prototypical charge transfer species. EOM-IP-CCSD describes these doublet systems based on the closed-shell reference and thus avoids the doublet instability problem. The studied quantities are associated with the quality of the potential energy surface (PES) along the charge transfer coordinate and distribution of the charge between fragments. It is found that EOM-IP-CCSD is capable of describing accurately both the charge-localized and charge-delocalized systems, yielding accurate charge distributions and energies. This is in stark contrast with the methods based on the open-shell reference, which overlocalize the charge and produce a PES cusp when the fragments are indistinguishable.

Energy landscapes of nucleophilic substitution reactions: A comparison of density functional theory and coupled cluster methods

Abstract: We have carried out a detailed evaluation of the performance of all classes of density functional theory (DFT) for describing the potential energy surface (PES) of a wide range of nucleophilic substitution (SN2) reactions involving, amongst others, nucleophilic attack at carbon, nitrogen, silicon, and sulfur. In particular, we investigate the ability of the local density approximation (LDA), generalized gradient approximation (GGA), meta-GGA as well as hybrid DFT to reproduce high-level coupled cluster (CCSD(T)) benchmarks that are close to the basis set limit. The most accurate GGA, meta-GGA, and hybrid functionals yield mean absolute deviations of about 2 kcal/mol relative to the coupled cluster data, for reactant complexation, central barriers, overall barriers as well as reaction energies. For the three nonlocal DFT classes, the best functionals are found to be OPBE (GGA), OLAP3 (meta-GGA), and mPBE0KCIS (hybrid DFT). The popular B3LYP functional is not bad but performs significantly worse than the best GGA functionals. Furthermore, we have compared the geometries from several density functionals with the reference CCSD(T) data. The same GGA functionals that perform best for the energies (OPBE, OLYP), also perform best for the geometries with average absolute deviations in bond lengths of 0.06 A and 0.6 degrees, even better than the best meta-GGA and hybrid functionals. In view of the reduced computational effort of GGAs with respect to meta-GGAs and hybrid functionals, let alone coupled cluster, we recommend the use of accurate GGAs such as OPBE or OLYP for the study of SN2 reactions.

Basis set limit coupled cluster study of h-bonded systems and assessment of more approximate methods.

Abstract: Hydrogen bonds are of utmost importance in both chemistry and biology. As the applicability of density functional theory and ab initio methods extends to ever larger systems and to liquids, an accurate description of such interactions is desirable. However, reference data are often lacking, and ab initio calculations are only possible and done in very small basis sets. Here, we present high level [CCSD(T)] ab initio reference calculations at the basis set limit on a large set of hydrogen-bonded systems and assess the accuracy of second-order perturbation theory (MP2). The possibilities of using basis set extrapolations for geometries and dissociation energies are discussed as well as the results of R12 methods and density functional and local correlation methods.

A Fock space coupled cluster study on the electronic structure of the UO2, UO2+, U4+, and U5+ species

Abstract: The ground and excited states of the UO2 molecule have been studied using a Dirac-Coulomb intermediate Hamiltonian Fock-space coupled cluster approach (DC-IHFSCC). This method is unique in describing dynamic and nondynamic correlation energies at relatively low computational cost. Spin-orbit coupling effects have been fully included by utilizing the four-component Dirac-Coulomb Hamiltonian from the outset. Complementary calculations on the ionized systems UO2+ and UO22+ as well as on the ions U4+ and U5+ were performed to assess the accuracy of this method. The latter calculations improve upon previously published theoretical work. Our calculations confirm the assignment of the ground state of the UO2 molecule as a Φ2u3 state that arises from the 5f17s1 configuration. The first state from the 5f2 configuration is found above 10000cm−1, whereas the first state from the 5f16d1 configuration is found at 5047cm−1.

Quasiclassical Trajectory Calculations of Acetaldehyde Dissociation on a Global Potential Energy Surface Indicate Significant Non-transition State Dynamics

Abstract: A recent experimental study [Houston, P. L.; Kable, S. H. Proc. Natl. Acad. Sci. U.S.A. 2006, 103, 16079] of the photodissociation of acetaldehyde (CH3CHO) has suggested two distinct mechanisms for the production of the molecular products CH4 + CO. One corresponds to the traditional transition state mechanism and the other to a transition state-skirting path similar to the roaming channel previously reported in formaldehyde. To investigate this theoretically, a full-dimensional potential energy surface (PES) has been constructed. The PES was fit with permutationally invariant polynomials to 135 000 points calculated using coupled cluster theory with single and double excitations and a perturbative treatment of triple excitations [CCSD(T)] and correlation consistent basis sets of double- and triple-ζ quality. To test the accuracy of the PES additional CCSD(T) and multireference configuration interaction calculations were carried out. Quasiclassical trajectory calculations were run on the PES starting at th...

Accurate calculation of anharmonic vibrational frequencies of medium sized molecules using local coupled cluster methods.

Abstract: Local coupled cluster methods were applied for the automated generation of accurate multidimensional potential energy surfaces for a set of test molecules ranging from six to nine atoms. Based on these surfaces anharmonic fundamental frequencies were computed using vibrational self-consistent field and configuration interaction methods. The computed vibrational frequencies are compared to those obtained from similar calculations using conventional coupled cluster methods and to experimental values. The results from local and conventional methods are found to be of similar accuracy and in close agreement with experimental values. In addition, an efficient parallelization of the fully automated surface generation code is presented.