Implementation and application of basis set superposition error-correction schemes to the theoretical modeling of weak intermolecular interactions

This thesis deals with the so-called Basis Set Superposition Error (BSSE) from both a methodological and a practical point of view. The purpose of the present thesis is twofold: (a) to contribute step ahead in the correct characterization of weakly bound complexes and, (b) to shed light the understanding of the actual implications of the basis set extension effects in the ab intio calculations and contribute to the BSSE debate. The existing BSSE-correction procedures are deeply analyzed, compared, validated and, if necessary, improved. A new interpretation of the counterpoise (CP) method is used in order to define counterpoise-corrected descriptions of the molecular complexes. This novel point of view allows for a study of the BSSE-effects not only in the interaction energy but also on the potential energy surface and, in general, in any property derived from the molecular energy and its derivatives A program has been developed for the calculation of CP-corrected geometry optimizations and vibrational frequencies, also using several counterpoise schemes for the case of molecular clusters. The method has also been implemented in Gaussian98 revA10 package. The Chemical Hamiltonian Approach (CHA) methodology has been also implemented at the RHF and UHF levels of theory for an arbitrary number interacting systems using an algorithm based on block-diagonal matrices. Along with the methodological development, the effects of the BSSE on the properties of molecular complexes have been discussed in detail. The CP and CHA methodologies are used for the determination of BSSE-corrected molecular complexes properties related to the Potential Energy Surfaces and molecular wavefunction, respectively. First, the behaviour of both BSSE-correction schemes are systematically compared at different levels of theory and basis sets for a number of hydrogen-bonded complexes. The Complete Basis Set (CBS) limit of both uncorrected and CP-corrected molecular properties like stabilization energies and intermolecular distances has also been determined, showing the capital importance of the BSSE correction. Several controversial topics of the BSSE correction are addressed as well. The application of the counterpoise method is applied to internal rotational barriers. The importance of the nuclear relaxation term is also pointed out. The viability of the CP method for dealing with charged complexes and the BSSE effects on the double-well PES blue-shifted hydrogen bonds is also studied in detail. In the case of the molecular clusters the effect of high-order BSSE effects introduced with the hierarchical counterpoise scheme is also determined. The effect of the BSSE on the electron density-related properties is also addressed. The first-order electron density obtained with the CHA/F and CHA/DFT methodologies was used to assess, both graphically and numerically, the redistribution of the charge density upon BSSE-correction. Several tools like the Atoms in Molecules topologycal analysis, density difference maps, Quantum Molecular Similarity, and Chemical Energy Component Analysis were used to deeply analyze, for the first time, the BSSE effects on the electron density of several hydrogen bonded complexes of increasing size. The indirect effect of the BSSE on intermolecular perturbation theory results is also pointed out It is shown that for a BSSE-free SAPT study of hydrogen fluoride clusters, the use of a counterpoise-corrected PES is essential in order to determine the proper molecular geometry to perform the SAPT analysis.

Gaussian basis sets to the theoretical study of the electronic structure of perovskite (LaMnO3)

The Generator Coordinate Hartree-Fock (GCHF) method is applied to generate extended (20s14p), (30s19p13d), and (31s23p18d) Gaussian basis sets for the 0, Mn, and La atoms, respectively. The role of the weight functions (WFs) in the assessment of the numerical integration range of the GCHF equations is shown. These basis sets are then contracted to [5s3p] and [11s6p6d] for 0 and Mn atoms, respectively, and [17s11p7d] for La atom by a standard procedure. For quality evaluation of contracted basis sets in molecular calculations, we have accomplished calculations of total and orbital energies in the Hartree-Fock-Roothaan (HFR) method for (MnO1+)-Mn-5 and (LaO1+)-La-1 fragments. The results obtained with the contracted basis sets are compared with values obtained with the extended basis sets. The addition of one d polarization function in the contracted basis set for 0 atom and its utilization with the contracted basis sets for Mn and La atoms leads to the calculations of dipole moment and total atomic charges of perovskite (LaMnO3). The calculations were performed at the HFR level with the crystal [LaMnO3](2) fragment in space group C-2v the values of dipole moment, total energy, and total atomic charges showed that it is reasonable to believe that LaMnO3 presents behaviour of piezoelectric material. (C) 2003 Elsevier B.V. All rights reserved.

Design of Gaussian basis sets to the theoretical interpretation of IR-spectrum of hexaaquachromium (III) ion, tetraoxochromium (IV) ion, and tetraoxochromium (VI) ion

The Generator Coordinate Hartree-Fock (GCHF) method is employed to design 16s, 16s10p, 24s17p13d, 25s17p13d, and 26s17p Gaussian basis sets for the H ((2)S), O ((3)P), O(2-) ((1)S), Cr(3+) ((4)F), Cr(4+) ((3)F), and Cr(6+) ((1)S) atomic species. These basis sets are then contracted to (4s) for H ((2)S), (6s4p) for O ((3)P), and O(2-) ((1)S), (986p3d) for Cr(3+) ((4)F), (10s8p3d) for Cr(4+) ((3)F), and (13s7p) for Cr(6+) (1S) by a standard procedure. For evaluation of the quality of those basis sets in molecular calculations, we have accomplished studies of total and orbital (HOMO and HOMO-1) energies at the HF-Roothaan level for the molecular species of our interest. The results obtained with the contracted basis sets are compared to the values obtained with our extended basis sets and to the standard 6-311G basis set from literature. Finally, the contracted basis sets are enriched with polarization function and then utilized in the theoretical interpretation of IR-spectrum of hexaaquachromium (III) ion, [Cr(H(2)O)(6)](3+), tetraoxochromium (IV) ion, [CrO(4)](4-), and tetraoxochromium (VI) ion, [CrO(4)](2-). The respective theoretical harmonic frequencies and IR-intensities were computed at the density functional theory (DFT) level. In the DFT calculations we employed the Becke's 1988 functional using the LYP correlation functional. The comparison between the results obtained and the corresponding experimental values indicates a very good description of the IR-spectra of the molecular ions studied, and that the GCHF method is still a legitimate alternative for selection of Gaussian basis sets. (C) 2003 Elsevier B.V. All rights reserved.

Basis sets applied to the theoretical study of the vibrational structure of hexaaquaaluminum(III) ion

Gaussian basis sets were developed with the Generator Coordinate Hartree-Fock (GCHF) method for the atoms from H (14s), O (23s16p), and Al (29sl9p) in the ground state. These basis sets were then contracted to 3s (12,1,1), 5s3p (18,2,1,1,1/14,1,1), and 7s5p (20,3,2,1,1,1,1/14,2,1,1,1) for H, O and Al atoms, respectively, by a standard procedure. The quality of contracted basis sets in molecular calculations was evaluated through studies of the total and orbital (epsilon(HOMO) and epsilon(HOMO-1)) energies at the HF level for the hexaaquaaluminum(III) ion, [AI(H(2)O)(6)](3+). For the O atom, the 5s3p was supplemented with d polarization function and it was used in combination with 3s, and 7s5p for H and Al atoms was used to the theoretical interpretation of the Infrared (IR) spectrum of hexaaquaaluminum(III) ion. The calculations of the IR-spectrum were also performed at the HF level and it showed that the basis sets obtained with the aid of GCHF method lead to the selection of useful contracted Gaussian basis sets for the theoretical study of vibrational property of ionic specie of our interest. (C) 2004 Elsevier B.V. All rights reserved.

Designing Gaussian basis sets to the theoretical study of the piezoelectric effect of perovskite (BaTiO3)

Gaussian basis sets (24s14p, 30s19p14d, and 33s21p14d for O (P-3), Ti (S-5), and Ba (S-1) atoms, respectively), are designed with the strategy of the Generator Coordinate Hartree-Fock method. The basis sets are then contracted to [6s4p], [10s5p4d], and [16s9p5d] to O, Ti, and Ba atoms, respectively, and used in calculations of total and orbital energies of (TiO+2)-Ti-1 and (BaO)-Ba-1 fragments for quality evaluation in molecular studies. For O atom, the [6s4p] basis set is enriched with d polarization function and used along with the [10s5p4d] and [16s9p5d] basis sets for the theoretical study of the piezoelectric effect of perovskite (BaTiO3). The results of this work evidence that the piezoelectric properties in BaTiO3 can be caused by electrostatic interactions. (C) 2003 Elsevier B.V. All rights reserved.

Contracted GTF basis sets applied to the theoretical interpretation of the Raman spectrum of hexaaquachromium(III) ion

Contracted GTF basis sets designed with aid of the Generator Coordinate Hartree-Fock (GCHF) method for H(2S), O2-(1S), and Cr3+(4F) atomic species are applied to perform theoretical interpretation of the Raman spectrum of hexaaquachromium(III) ion. The 16s, 16s 10p, and 24s17p13d GTF basis sets were contracted to [4s] for H atom, [6s4p], and [9s6p3d] for O2- and Cr3+, respectively, by Dunning's scheme. For Cr3+, the [9s6p3d] basis set was enriched with f polarization function and used in combination com [4s] and [6s4p] in the study of our interest. The results obtained in this report show that the contracted GTF basis sets used are a useful alternative for the theoretical interpretation of Raman spectrum of hexaaquachromium(III) ion and that GCHF method is an effective alternative to selection of GTF basis sets for theoretical study of vibrational properties of poliatomic species. © 2003 Elsevier Science B.V. All rights reserved.

A constraints-led perspective to understanding skill acquisition and game play: A basis for integration of motor learning theory and physical education praxis?

Background: In order to design appropriate environments for performance and learning of movement skills, physical educators need a sound theoretical model of the learner and of processes of learning. In physical education, this type of modelling informs the organization of learning environments and effective and efficient use of practice time. An emerging theoretical framework in motor learning, relevant to physical education, advocates a constraints-led perspective for acquisition of movement skills and game play knowledge. This framework shows how physical educators could use task, performer and environmental constraints to channel acquisition of movement skills and decision making behaviours in learners. From this viewpoint, learners generate specific movement solutions to satisfy the unique combination of constraints imposed on them, a process which can be harnessed during physical education lessons. Purpose: In this paper the aim is to provide an overview of the motor learning approach emanating from the constraints-led perspective, and examine how it can substantiate a platform for a new pedagogical framework in physical education: nonlinear pedagogy. We aim to demonstrate that it is only through theoretically valid and objective empirical work of an applied nature that a conceptually sound nonlinear pedagogy model can continue to evolve and support research in physical education. We present some important implications for designing practices in games lessons, showing how a constraints-led perspective on motor learning could assist physical educators in understanding how to structure learning experiences for learners at different stages, with specific focus on understanding the design of games teaching programmes in physical education, using exemplars from Rugby Union and Cricket. Findings: Research evidence from recent studies examining movement models demonstrates that physical education teachers need a strong understanding of sport performance so that task constraints can be manipulated so that information-movement couplings are maintained in a learning environment that is representative of real performance situations. Physical educators should also understand that movement variability may not necessarily be detrimental to learning and could be an important phenomenon prior to the acquisition of a stable and functional movement pattern. We highlight how the nonlinear pedagogical approach is student-centred and empowers individuals to become active learners via a more hands-off approach to learning. Summary: A constraints-based perspective has the potential to provide physical educators with a framework for understanding how performer, task and environmental constraints shape each individual‟s physical education. Understanding the underlying neurobiological processes present in a constraints-led perspective to skill acquisition and game play can raise awareness of physical educators that teaching is a dynamic 'art' interwoven with the 'science' of motor learning theories.

Theoretical comparison of innovative window daylighting devices for a sub-tropical climate using radiance

Daylighting in tropical and sub-tropical climates presents a unique challenge that is generally not well understood by designers. In a sub-tropical region such as Brisbane, Australia the majority of the year comprises of sunny clear skies with few overcast days and as a consequence windows can easily become sources of overheating and glare. The main strategy in dealing with this issue is extensive shading on windows. However, this in turn prevents daylight penetration into buildings often causing an interior to appear gloomy and dark even though there is more than sufficient daylight available. As a result electric lighting is the main source of light, even during the day. Innovative daylight devices which redirect light from windows offer a potential solution to this issue. These devices can potentially improve daylighting in buildings by increasing the illumination within the environment decreasing the high contrast between the window and work regions and deflecting potentially glare causing sunlight away from the observer. However, the performance of such innovative daylighting devices are generally quantified under overcast skies (i.e. daylight factors) or skies without sun, which are typical of European climates and are misleading when considering these devices for tropical or sub-tropical climates. This study sought to compare four innovative window daylighting devices in RADIANCE; light shelves, laser cut panels, micro-light guides and light redirecting blinds. These devices were simulated in RADIANCE under sub-tropical skies (for Brisbane) within the test case of a typical CBD office space. For each device the quantity of light redirected and its distribution within the space was used as the basis for comparison. In addition, glare analysis on each device was conducted using Weinold and Christoffersons evalglare. The analysis was conducted for selected hours for a day in each season. The majority of buildings that humans will occupy in their lifetime are already constructed, and extensive remodelling of most of these buildings is unlikely. Therefore the most effective way to improve daylighting in the near future will be through the alteration existing window spaces. Thus it will be important to understand the performance of daylighting systems with respect to the climate it is to be used in. This type of analysis is important to determine the applicability of a daylighting strategy so that designers can achieve energy efficiency as well the health benefits of natural daylight.

Genesis and role of standards : theoretical foundations and socio-economical model for the construction and use of standards

A key issue in the economic development and performance of organizations is the existence of standards. Their definition and control are sources of power and it is important to understand their concept, as it gives standards their direction and their legitimacy, and to explore how they are represented and applied. The difficulties posed by classical micro-economics in establishing a theory of standardization that is compatible with its fundamental axiomatic are acknowledged. We propose to reconsider the problem by taking the opposite perspective in questioning its theoretical base and by reformulating assumptions about the independent and autonomous decisions taken by actors. The Theory of Conventions will offer us a theoretical framework and tools enabling us to understand the systemic dimension and dynamic structure of standards. These will be seen as a special case of conventions. This work aims to provide a sound basis and promote a better consciousness in the development of global project management standards. It aims also to emphasize that social construction is not a matter of copyright but a matter of open minds, collective cognitive process and freedom for the common wealth.

Genesis and role of standards : theoretical foundations and socio-economical model for the construction and use of standards

A key issue for the economic development and for performance of organizations is the existence of standards. As their definitions and control are source of power, it seems to be important to understand the concept and to wonder about the representations authorized by the concept which give their direction and their legitimacy. The difficulties of classical microeconomics of establishing a theory of standardisation compatible with its fundamental axiomatic are underlined. We propose to reconsider the problem by carrying out the opposite way: to question the theoretical base, by reformulating assumptions on the autonomy of the choice of the actors. The theory of conventions will offer us both a theoretical framework and tools, enabling us to understand the systemic dimension and dynamic structure of standards seen as special case of conventions. This work aims thus to provide a sound basis and promote a better consciousness in the development of global project management standards, aiming also to underline that social construction is not a matter of copyright but a matter of open minds, collective cognitive process and freedom for the common wealth.

Theoretical and experimental characterisation of energy in an electrostatic discharge

Electrostatic discharges have been identified as the most likely cause in a number of incidents of fire and explosion with unexplained ignitions. The lack of data and suitable models for this ignition mechanism creates a void in the analysis to quantify the importance of static electricity as a credible ignition mechanism. Quantifiable hazard analysis of the risk of ignition by static discharge cannot, therefore, be entirely carried out with our current understanding of this phenomenon. The study of electrostatics has been ongoing for a long time. However, it was not until the wide spread use of electronics that research was developed for the protection of electronics from electrostatic discharges. Current experimental models for electrostatic discharge developed for intrinsic safety with electronics are inadequate for ignition analysis and typically are not supported by theoretical analysis. A preliminary simulation and experiment with low voltage was designed to investigate the characteristics of energy dissipation and provided a basis for a high voltage investigation. It was seen that for a low voltage the discharge energy represents about 10% of the initial capacitive energy available and that the energy dissipation was within 10 ns of the initial discharge. The potential difference is greatest at the initial break down when the largest amount of the energy is dissipated. The discharge pathway is then established and minimal energy is dissipated as energy dissipation becomes greatly influenced by other components and stray resistance in the discharge circuit. From the initial low voltage simulation work, the importance of the energy dissipation and the characteristic of the discharge were determined. After the preliminary low voltage work was completed, a high voltage discharge experiment was designed and fabricated. Voltage and current measurement were recorded on the discharge circuit allowing the discharge characteristic to be recorded and energy dissipation in the discharge circuit calculated. Discharge energy calculations show consistency with the low voltage work relating to discharge energy with about 30-40% of the total initial capacitive energy being discharged in the resulting high voltage arc. After the system was characterised and operation validated, high voltage ignition energy measurements were conducted on a solution of n-Pentane evaporating in a 250 cm3 chamber. A series of ignition experiments were conducted to determine the minimum ignition energy of n-Pentane. The data from the ignition work was analysed with standard statistical regression methods for tests that return binary (yes/no) data and found to be in agreement with recent publications. The research demonstrates that energy dissipation is heavily dependent on the circuit configuration and most especially by the discharge circuit's capacitance and resistance. The analysis established a discharge profile for the discharges studied and validates the application of this methodology for further research into different materials and atmospheres; by systematically looking at discharge profiles of test materials with various parameters (e.g., capacitance, inductance, and resistance). Systematic experiments looking at the discharge characteristics of the spark will also help understand the way energy is dissipated in an electrostatic discharge enabling a better understanding of the ignition characteristics of materials in terms of energy and the dissipation of that energy in an electrostatic discharge.

A basis for instrusion detection in distributed systems using kernel-level data tainting

This project was a step forward in developing intrusion detection systems in distributed environments such as web services. It investigates a new approach of detection based on so-called "taint-marking" techniques and introduces a theoretical framework along with its implementation in the Linux kernel.

A practical and theoretical definition of very small field size for radiotherapy output factor measurements

Purpose This work introduces the concept of very small field size. Output factor (OPF) measurements at these field sizes require extremely careful experimental methodology including the measurement of dosimetric field size at the same time as each OPF measurement. Two quantifiable scientific definitions of the threshold of very small field size are presented. Methods A practical definition was established by quantifying the effect that a 1 mm error in field size or detector position had on OPFs, and setting acceptable uncertainties on OPF at 1%. Alternatively, for a theoretical definition of very small field size, the OPFs were separated into additional factors to investigate the specific effects of lateral electronic disequilibrium, photon scatter in the phantom and source occlusion. The dominant effect was established and formed the basis of a theoretical definition of very small fields. Each factor was obtained using Monte Carlo simulations of a Varian iX linear accelerator for various square field sizes of side length from 4 mm to 100 mm, using a nominal photon energy of 6 MV. Results According to the practical definition established in this project, field sizes < 15 mm were considered to be very small for 6 MV beams for maximal field size uncertainties of 1 mm. If the acceptable uncertainty in the OPF was increased from 1.0 % to 2.0 %, or field size uncertainties are 0.5 mm, field sizes < 12 mm were considered to be very small. Lateral electronic disequilibrium in the phantom was the dominant cause of change in OPF at very small field sizes. Thus the theoretical definition of very small field size coincided to the field size at which lateral electronic disequilibrium clearly caused a greater change in OPF than any other effects. This was found to occur at field sizes < 12 mm. Source occlusion also caused a large change in OPF for field sizes < 8 mm. Based on the results of this study, field sizes < 12 mm were considered to be theoretically very small for 6 MV beams. Conclusions Extremely careful experimental methodology including the measurement of dosimetric field size at the same time as output factor measurement for each field size setting and also very precise detector alignment is required at field sizes at least < 12 mm and more conservatively < 15 mm for 6 MV beams. These recommendations should be applied in addition to all the usual considerations for small field dosimetry, including careful detector selection.

Conversion of linear to rhombic C-4 in the gas phase: A joint experimental and theoretical study

Charge reversal (CR) and neutralization reionization (NR) experiments carried out on a 4-sector mass spectrometer demonstrate that isotopically labeled, linear C-4 anion rearranges upon collisional oxidation. The cations and neutrals formed in these experiments exhibit differing degrees of isotopic scrambling in their fragmentation patterns, indicative of (at least) partial isomerization of both states. Theoretical studies, employing the CCSD(T)/aug-cc-pVDZ//B3LYP/6-31G(d) level of theory, favor conversion to the rhombic C-4 isomer on both cationic and neutral potential-energy surfaces with the rhombic structures predicted to be slightly more stable than the linear forms in each case. The combination of experiment with theory indicates that the elusive rhombic C-4 is formed as a cation and as a neutral following charge stripping of linear C-4(-)