Exact property estimation from diffusion Monte Carlo with minimal stochastic reconfiguration /

Our objective is to develop a diffusion Monte Carlo (DMC) algorithm to estimate the exact expectation values, ($o|^|^o), of multiplicative operators, such as polarizabilities and high-order hyperpolarizabilities, for isolated atoms and molecules. The existing forward-walking pure diffusion Monte Carlo (FW-PDMC) algorithm which attempts this has a serious bias. On the other hand, the DMC algorithm with minimal stochastic reconfiguration provides unbiased estimates of the energies, but the expectation values ($o|^|^) are contaminated by ^, an user specified, approximate wave function, when A does not commute with the Hamiltonian. We modified the latter algorithm to obtain the exact expectation values for these operators, while at the same time eliminating the bias. To compare the efficiency of FW-PDMC and the modified DMC algorithms we calculated simple properties of the H atom, such as various functions of coordinates and polarizabilities. Using three non-exact wave functions, one of moderate quality and the others very crude, in each case the results are within statistical error of the exact values.

Calculation of the magnetic field penetration depth for high-Tc cuprate superconductors based on the Interlayer Pair Tunneling model /

In this work, the magnetic field penetration depth for high-Tc cuprate superconductors is calculated using a recent Interlayer Pair Tunneling (ILPT) model proposed by Chakravarty, Sudb0, Anderson, and Strong [1] to explain high temperature superconductivity. This model involves a "hopping" of Cooper pairs between layers of the unit cell which acts to amplify the pairing mechanism within the planes themselves. Recent work has shown that this model can account reasonably well for the isotope effect and the dependence of Tc on nonmagnetic in-plane impurities [2] , as well as the Knight shift curves [3] and the presence of a magnetic peak in the neutron scattering intensity [4]. In the latter case, Yin et al. emphasize that the pair tunneling must be the dominant pairing mechanism in the high-Tc cuprates in order to capture the features found in experiments. The goal of this work is to determine whether or not the ILPT model can account for the experimental observations of the magnetic field penetration depth in YBa2Cu307_a7. Calculations are performed in the weak and strong coupling limits, and the efi"ects of both small and large strengths of interlayer pair tunneling are investigated. Furthermore, as a follow up to the penetration depth calculations, both the neutron scattering intensity and the Knight shift are calculated within the ILPT formalism. The aim is to determine if the ILPT model can yield results consistent with experiments performed for these properties. The results for all three thermodynamic properties considered are not consistent with the notion that the interlayer pair tunneling must be the dominate pairing mechanism in these high-Tc cuprate superconductors. Instead, it is found that reasonable agreement with experiments is obtained for small strengths of pair tunneling, and that large pair tunneling yields results which do not resemble those of the experiments.

A treatment of atomic mean square displacement in higher order perturbation theory

A general derivation of the anharmonic coefficients for a periodic lattice invoking the special case of the central force interaction is presented. All of the contributions to mean square displacement (MSD) to order 14 perturbation theory are enumerated. A direct correspondance is found between the high temperature limit MSD and high temperature limit free energy contributions up to and including 0(14). This correspondance follows from the detailed derivation of some of the contributions to MSD. Numerical results are obtained for all the MSD contributions to 0(14) using the Lennard-Jones potential for the lattice constants and temperatures for which the Monte Carlo results were calculated by Heiser, Shukla and Cowley. The Peierls approximation is also employed in order to simplify the numerical evaluation of the MSD contributions. The numerical results indicate the convergence of the perturbation expansion up to 75% of the melting temperature of the solid (TM) for the exact calculation; however, a better agreement with the Monte Carlo results is not obtained when the total of all 14 contributions is added to the 12 perturbation theory results. Using Peierls approximation the expansion converges up to 45% of TM• The MSD contributions arising in the Green's function method of Shukla and Hubschle are derived and enumerated up to and including 0(18). The total MSD from these selected contributions is in excellent agreement with their results at all temperatures. Theoretical values of the recoilless fraction for krypton are calculated from the MSD contributions for both the Lennard-Jones and Aziz potentials. The agreement with experimental values is quite good.

Molecular dynamics calculation of mean square displacement in alkali metals and rare gas solids and comparison with lattice dynamics

Molec ul ar dynamics calculations of the mean sq ua re displacement have been carried out for the alkali metals Na, K and Cs and for an fcc nearest neighbour Lennard-Jones model applicable to rare gas solids. The computations for the alkalis were done for several temperatures for temperature vol ume a swell as for the the ze r 0 pressure ze ro zero pressure volume corresponding to each temperature. In the fcc case, results were obtained for a wide range of both the temperature and density. Lattice dynamics calculations of the harmonic and the lowe s t order anharmonic (cubic and quartic) contributions to the mean square displacement were performed for the same potential models as in the molecular dynamics calculations. The Brillouin zone sums arising in the harmonic and the quartic terms were computed for very large numbers of points in q-space, and were extrapolated to obtain results ful converged with respect to the number of points in the Brillouin zone.An excellent agreement between the lattice dynamics results was observed molecular dynamics and in the case of all the alkali metals, e~ept for the zero pressure case of CSt where the difference is about 15 % near the melting temperature. It was concluded that for the alkalis, the lowest order perturbation theory works well even at temperat ures close to the melting temperat ure. For the fcc nearest neighbour model it was found that the number of particles (256) used for the molecular dynamics calculations, produces a result which is somewhere between 10 and 20 % smaller than the value converged with respect to the number of particles. However, the general temperature dependence of the mean square displacement is the same in molecular dynamics and lattice dynamics for all temperatures at the highest densities examined, while at higher volumes and high temperatures the results diverge. This indicates the importance of the higher order (eg. ~* ) perturbation theory contributions in these cases.

On the formulation and the calculation of the harmonic contributions to the Debye-Waller factor in metals (sodium)

The algebraic expressions for the anharmonic contributions to the Debye-Waller factor up to 0(A ) and 0 L% ) £ where ^ is the scattering wave-vector] have been derived in a form suitable for cubic metals with small ion cores where the interatomic potential extends to many neighbours. This has been achieved in terms of various wave-vector dependent tensors, following the work of Shukla and Taylor (1974) on the cubic anharmonic Helmholtz free energy. The contribution to the various wave-vector dependent tensors from the coulomb and the electron-ion terms in the interatomic metallic potential has been obtained by the Ewald procedure. All the restricted multiple whole B r i l l o u i n zone (B.Z.) sums are reduced to single whole B.Z. sums by using the plane wave representation of the delta function. These single whole B.Z. sums are further reduced to the •%?? portion of the B.Z. following Shukla and Wilk (1974) and Shukla and Taylor (1974). Numerical calculations have been performed for sodium where the Born-Mayer term in the interatomic potential has been neglected because i t is small £ Vosko (1964)3 • *n o^er to compare our calculated results with the experimental results of Dawton (1937), we have also calculated the r a t io of the intensities at different temperatures for the lowest five reflections (110), (200), (220), (310) and (400) . Our calculated quasi-harmonic results agree reasonably well with the experimental results at temperatures (T) of the order of the Debye temperature ( 0 ). For T » © ^ 9 our calculated anharmonic results are found to be in good agreement with the experimental results.The anomalous terms in the Debye-Waller factor are found not to be negligible for certain reflections even for T ^ ©^ . At temperature T yy Op 9 where the temperature is of the order of the melting temperature (Xm) » "the anomalous terms are found to be important almost for all the f i ve reflections.

The principal in the high school: effectiveness in a climate of sharing

This pilot study developed a climate instrument which was administered in a sample of high schools in one board of education. Several tests were conducted i n order to determine the reliability and internal consistency of the instrument . The ability of the instrument to identify the demographic differences of school and gender was also tested. The relationship between leadership styles and an effective use of authority in creating a productive and rewarding work environment was the f ocus of t his study. Attitudes to leadership and perceived school morale were investigated in a demographic study, a climate survey, as well as a body of related literature. In light of the empirical research, an attempt was made to determine the extent to which the authority figure's behaviour and adopted leadership style contributed to a positive school climate : one in which t eachers were motivated to achieve to t he best of their abilities by way of their commitment and service. The tone of authority assumed by t he leader not only shapes the mood of the school environment but ultimately determines the efficiency and morale of t he teaching staff.

An Interactive Theorem Prover for First-Order Dynamic Logic

Dynamic logic is an extension of modal logic originally intended for reasoning about computer programs. The method of proving correctness of properties of a computer program using the well-known Hoare Logic can be implemented by utilizing the robustness of dynamic logic. For a very broad range of languages and applications in program veri cation, a theorem prover named KIV (Karlsruhe Interactive Veri er) Theorem Prover has already been developed. But a high degree of automation and its complexity make it di cult to use it for educational purposes. My research work is motivated towards the design and implementation of a similar interactive theorem prover with educational use as its main design criteria. As the key purpose of this system is to serve as an educational tool, it is a self-explanatory system that explains every step of creating a derivation, i.e., proving a theorem. This deductive system is implemented in the platform-independent programming language Java. In addition, a very popular combination of a lexical analyzer generator, JFlex, and the parser generator BYacc/J for parsing formulas and programs has been used.

High-Spin and Emissive Molecular Materials: Synthesis and Characterization of New Polynuclear Ni(II) Complexes from the Use of Aromatic Schiff Bases as Bridging Ligands

The employment of the bridging/chelating Schiff bases, N-salicylidene-4-methyl-o-aminophenol (samphH2) and N-naphthalidene-2-amino-5-chlorobenzoic acid (nacbH2), in nickel cluster chemistry has afforded eight polynuclear Ni(II) complexes with new structural motifs, interesting magnetic and optical properties, and unexpected organic ligand transformations. In the present thesis, Chapter 1 deals with all the fundamental aspects of polynuclear metal complexes, molecular magnetism and optics, while research results are reported in Chapters 2 and 3. In the first project (Chapter 2), I investigated the coordination chemistry of the organic chelating/bridging ligand, N-salicylidene-4-methyl-o-aminophenol (samphH2). The general NiII/tBuCO2-/samphH2 reaction system afforded two new tetranuclear NiII clusters, namely [Ni4(samph)4(EtOH)4] (1) and [Ni4(samph)4(DMF)2] (2), with different structural motifs. Complex 1 possessed a cubane core while in complex 2 the four NiII ions were located at the four vertices of a defective dicubane. The nature of the organic solvent was found to be of pivotal importance, leading to compounds with the same nuclearity, but different structural topologies and magnetic properties. The second project, the results of which are summarized in Chapter 3, included the systematic study of a new optically-active Schiff base ligand, N-naphthalidene-2-amino-5-chlorobenzoic acid (nacbH2), in NiII cluster chemistry. Various reactions between NiX2 (X- = inorganic anions) and nacbH2 were performed under basic conditions to yield six new polynuclear NiII complexes, namely (NHEt3)[Ni12(nacb)12(H2O)4](ClO4) (3), (NHEt3)2[Ni5(nacb)4(L)(LH)2(MeOH)] (4), [Ni5(OH)2(nacb)4(DMF)4] (5), [Ni5(OMe)Cl(nacb)4(MeOH)3(MeCN)] (6), (NHEt3)2[Ni6(OH)2(nacb)6(H2O)4] (7), and [Ni6(nacb)6(H2O)3(MeOH)6] (8). The nature of the solvent, the inorganic anion, X-, and the organic base were all found to be of critical importance, leading to products with different structural topologies and nuclearities (i.e., {Ni5}, {Ni6} and {Ni12}). Magnetic studies on all synthesized complexes revealed an overall ferromagnetic behavior for complexes 4 and 8, with the remaining complexes being dominated by antiferromagnetic exchange interactions. In order to assess the optical efficiency of the organic ligand when bound to the metal centers, photoluminescence studies were performed on all synthesized compounds. Complexes 4 and 5 show strong emission in the visible region of the electromagnetic spectrum. Finally, the ligand nacbH2 allowed for some unexpected organic transformations to occur; for instance, the pentanuclear compound 5 comprises both nacb2- groups and a new organic chelate, namely the anion of 5-chloro-2-[(3-hydroxy-4-oxo-1,4-dihydronaphthalen-1-yl)amino]benzoic acid. In the last section of this thesis, an attempt to compare the NiII cluster chemistry of the N-naphthalidene-2-amino-5-chlorobenzoic acid ligand with that of the structurally similar but less bulky, N-salicylidene-2-amino-5-chlorobenzoic acid (sacbH2), was made.