title: Microwave detected, microwave-optical double resonance of NH3, NH2D, NHD2, and ND3. I. Sturcture and force field of the A state

abstract: Microwave detected, microwave-optical double resonance was used to record the A state electronic spectrum of NH3, NH2D, and NHD2 with both vibrational and rotational resolution. To investigate ND3 with the same resolution as we had with our hydrogen containing isotopomers, a strip-line cell was constructed allowing the simultaneous passage of radio-frequency and ultraviolet radiation. Rotational constants were obtained as a function of nu2 excitation and an A state equilibrium bond length was estimated at 1.055(8) angstroms. In addition, the harmonic force field for the A state has been experimentally determined. fhh, falpha,alpha-falpha,alpha', and frr were found to be 1.06(4)aJ/angstrom2, 0.25(2)aJ, and 4.9aJ/angstrom^2 respectively. This calculated harmonic force field predicts that the asymmetry observed in the NH3 2(4) band is due to a strong anharmonic interaction with the 4(3) level and the broad feature observed in the dispersed fluorescence spectrum previously assigned to the 1(1) band is more likely attributable to the 4(2) level.


title: Specific heat in the thermodynamics of clusters

abstract: The thermodynamic properties such as the specific heat are uniquely determined by the second moments of the energy distribution for a given ensemble averaging. However for small particle numbers the results depend on the ensemble chosen. We calculated the higher moments of the distributions of some observables for both the canonical and the microcanonical ensemble of the same van der Waals clusters. The differences of the resulting thermodynamic observables for the two ensembles are calculated in terms of the higher moments. We demonstrate how for increasing particle number these terms decrease to vanish for bulk material. For the calculation of the specific heat within the microcanonical ensemble we give a new method based on an analysis of histograms.


title: Two-Color Coherent Photodissociation of Nitrogen Oxide in Intense Laser Fields

abstract: A simple one-dimensional semi-classical model with a Morse potential is used to investigate the possibility of two-color infrared multi-photon dissociation of vibrationally excited nitrogen oxide. The amplitude ratio effects and adiabatic effects are investigated. Some initial states are found to have thresholds smaller than expected from single-mode considerations and multiple thresholds exist for initial states up to 32. PACS: 42.50.Hz


title: Integral equation theory for the electrode-electrolyte interface with the central force water model. Results for an aqueous solution of sodium chloride

abstract: The structure of an aqueous solution of sodium chloride at a planar electrode is investigated by integral equation techniques. With the central force water model the aqueous electrolyte is modelled as a mixture of sodium and chloride ions and partially charged hydrogen and oxygen atoms interacting via effective spherically symmetric pair potentials. The bulk correlation functions are obtained from the Ornstein-Zernike equation. With the Wertheim-Lovett-Mou-Buff equation we have calculated the density profiles at the uncharged and charged electrode. Steric interactions between the differently sized ions and the ice-like water structure near the electrode dominates the ionic distribution. This model electrolyte also responds differently to opposite charges on the electrode.


title: Structure of Si12 Cluster

abstract: Tight-binding molecular dynamic simulations have revealed that Si$_{12}$ is an icosahedron with all atoms on the surface of an approximately 5 \AA~ diameter sphere. This is the most spherical cage structure for silicon clusters in the 2-13 atom size range.


title: Structure of Silicon Clusters

abstract: We determined the structures of silicon clusters in the 11-14 atom size range using the tight-binding molecular dynamics method. These calculations reveal that \Si{11} is an icosahedron with one missing cap, \Si{12} is a complete icosahedron, \Si{13} is a surface capped icosahedron, and \Si{14} is a 4-4-4 layer structure with two caps. The characteristic feature of these clusters is that they are all surface.


title: Spectral Shifts of Semiconductor Clusters

abstract: The shifts of the electronic absorption spectra of GaAs and GaP semiconductor clusters are calculated using accurate pseudopotentials. In the absence of experimental data at present, these calculations provide estimates of expected spectral shifts in these clusters. In addition, these calculations show that Coulomb interaction between the electron and hole dominates over the confinement energy in small clusters, with the result that the electronic absorption spectra of small clusters exhibit redshift instead of blueshift as the cluster size is decreased.


title: "The numerical accuracy of truncated Ewald sums for periodic systems with long-range Coulomb interactions"

abstract: Ewald summation is widely used to calculate electrostatic interactions in computer simulations of condensed-matter systems. We present an analysis of the errors arising from truncating the infinite real- and Fourier-space lattice sums in the Ewald formulation. We derive an optimal choice for the Fourier-space cutoff given a screening parameter $\eta$. We find that the number of vectors in Fourier space required to achieve a given accuracy scales with $\eta^3$. The proposed method can be used to determine computationally efficient parameters for Ewald sums, to assess the quality of Ewald-sum implementations, and to compare different implementations.


title: Critical Behavior of the Widom-Rowlinson Lattice Model

abstract: We report extensive Monte Carlo simulations of the Widom-Rowlinson lattice model in two and three dimensions. Our results yield precise values for the critical activities and densities, and clearly place the critical behavior in the Ising universality class.


title: Sequence Dependence of Self-Interacting Random Chains

abstract: We study the thermodynamic behavior of the random chain model proposed by Iori, Marinari and Parisi, and how this depends on the actual sequence of interactions along the chain. The properties of randomly chosen sequences are compared to those of designed ones, obtained through a simulated annealing procedure in sequence space. We show that the transition to the folded phase takes place at a smaller strength of the quenched disorder for designed sequences. As a result, folding can be relatively fast for these sequences.


title: A Variational Approach for Minimizing Lennard-Jones Energies

abstract: A variational method for computing conformational properties of molecules with Lennard-Jones potentials for the monomer-monomer interactions is presented. The approach is tailored to deal with angular degrees of freedom, {\it rotors}, and consists in the iterative solution of a set of deterministic equations with annealing in temperature. The singular short-distance behaviour of the Lennard-Jones potential is adiabatically switched on in order to obtain stable convergence. As testbeds for the approach two distinct ensembles of molecules are used, characterized by a roughly dense-packed ore a more elongated ground state. For the latter, problems are generated from natural frequencies of occurrence of amino acids and phenomenologically determined potential parameters; they seem to represent less disorder than was previously assumed in synthetic protein studies. For the dense-packed problems in particular, the variational algorithm clearly outperforms a gradient descent method in terms of minimal energies. Although it cannot compete with a careful simulating annealing algorithm, the variational approach requires only a tiny fraction of the computer time. Issues and results when applying the method to polyelectrolytes at a finite temperature are also briefly discussed.


title: Application of Gradient-Corrected Density Functional Theory to the Structures and Thermochemistries of $\mbox{ScF}_3$, $\mbox{TiF}_4$, $\mbox{VF}_5$, and $\mbox{CrF}_6$

abstract: Density functional theory(DFT) and Hartree-Fock(HF) calculations are reported for the family of transition metal fluorides $\mbox{ScF}_3$, $\mbox{TiF}_4$, $\mbox{VF}_5$, and $\mbox{CrF}_6$. Both HF and the local-density-aproximation (LDA) yield excellent agreement with experimental bond lengths, while the B-LYP gradient-corrected density functional gives bond lengths $0.04-0.05$ \AA\ too long. An investigation of various combinations of exchange and correlation functionals shows that, for this series, the origin of this behavior lies in the Becke exchange functional. Much improved bond distances are found using the hybrid HF/DFT functional advocated by Becke. This approximation also leads to much improved thermochemistries. The LDA overestimates average bond energies in this series by $30-40$ kcal/mol, whereas the B-LYP functional overbinds by only $\sim8-12$ kcal/mol, and the hybrid HF/DFT method overbinds by only $\sim 2$ kcal/mol. The hybrid method predicts the octahedral isomer of $\mbox{CrF}_6$ to be more stable than the trigonal prismatic form by $14$ kcal/mol. Comparison of theoretical vibrational frequencies with experiment supports the assignment of an octahedral geometry.


title: A Monte Carlo study of temperature-programmed desorption spectra with attractive lateral interactions

abstract: We present results of a Monte Carlo study of temperature-programmed desorption in a model system with attractive lateral interactions. It is shown that even for weak interactions there are large shifts of the peak maximum temperatures with initial coverage. The system has a transition temperature below which the desorption has a negative order. An analytical expression for this temperature is derived. The relation between the model and real systems is discussed.


title: Nonlinearity and trapping in excitation transfer: Dimers and Trimers.

abstract: We study the interplay between nonlinearity in exciton transport and trapping due to a sink site for the dimer and the trimer with chain configuration by a numerical integration of the discrete nonlinear Schroedinger equation. Our results for the dimer show, that the formation of a self trapped state due to the nonlinear coupling increases the life time of the exciton substantially. Self trapping can be enhanced by the sink for short times, but for long times it disappears. In the trimer consisting of a subdimer extended by a sink site exists a transition between states localized on the two sites of the subdimer before for larger nonlinear coupling self trapping on one site of the subdimer is observed. For large trapping rates the fear of death effect leads to an increasing life time of the excitation on both, the dimer and the trimer. The sink site is then effectively decoupled. We explain this effect using an asymptotic theory for strong trapping and demonstrate it by direct numerical computation.


title: Excitonic - vibronic coupled dimers: A dynamic approach

abstract: The dynamical properties of exciton transfer coupled to polarization vibrations in a two site system are investigated in detail. A fixed point analysis of the full system of Bloch - oscillator equations representing the coupled excitonic - vibronic flow is performed. For overcritical polarization a bifurcation converting the stable bonding ground state to a hyperbolic unstable state which is basic to the dynamical properties of the model is obtained. The phase space of the system is generally of a mixed type: Above bifurcation chaos develops starting from the region of the hyperbolic state and spreading with increasing energy over the Bloch sphere leaving only islands of regular dynamics. The behaviour of the polarization oscillator accordingly changes from regular to chaotic.


title: Temperature Dependence of the Optical Response of Small Sodium Clusters

abstract: We present an analysis of the temperature dependence of the optical response of small sodium clusters in a temperature range bracketing the melting phase transition. When the temperature increases, the mean excitation energy undergoes a red shift and the plasmon is significantly broadened, in agreement with recent experimental data. We show that the single--particle levels acquire a prominent width and the HOMO--LUMO gap as well as the width of the occupied band are reduced due to large thermal cluster size and shape fluctuations. This results in a sharp increase of the static polarizability with temperature.


title: Metal-Insulator Transition in C60-Polymers

abstract: Variations in the band structures of C60-polymers are studied, when pi-conjugation conditions are changed. We look at band structures in order to discuss a metal-insulator transition, using a semi-empirical model with the Su-Schrieffer-Heeger type electron-phonon interactions. We find that electronic structures change among direct-gap insulators and the metal, depending on the degree of pi-conjugations. High pressure experiments could observe such pressure-induced metal-insulator transitions.

categories: chem-ph cond-mat mtrl-th physics.chem-ph


title: Application of A Distributed Nucleus Approximation In Grid Based Minimization of the Kohn-Sham Energy Functional

abstract: In the distributed nucleus approximation we represent the singular nucleus as smeared over a smallportion of a Cartesian grid. Delocalizing the nucleus allows us to solve the Poisson equation for theoverall electrostatic potential using a linear scaling multigrid algorithm.This work is done in the context of minimizing the Kohn-Sham energy functionaldirectly in real space with a multiscale approach. The efficacy of the approximation is illustrated bylocating the ground state density of simple one electron atoms and moleculesand more complicated multiorbital systems.


title: Density functional theory using an optimized exchange-correlation potential

abstract: We have performed self-consistent calculations for first and second row atoms using a variant of density-functional theory, the optimized effective potential method, with an approximation due to Krieger, Li and Iafrate and a correlation-energy functional developed by Colle and Salvetti. The mean absolute deviation of first-row atomic ground-state energies from the exact non-relativistic values is 4.7 mH in our scheme, as compared to 4.5 mH in a recent configuration-interaction calculation. The proposed scheme is significantly more accurate than the conventional Kohn-Sham method while the numerical effort involved is about the same as for an ordinary Hartree-Fock calculation.


title: On finite-size effects in computer simulations using the Ewald potential

abstract: We discuss the origin and relevance for computer simulations of a strong finite-size effect that appears when using the Ewald summation formula. It can be understood as arising from a volume-dependent shift of the potential in a finite, periodic box relative to the infinite volume limit. This shift is due to the fact that the ``zero of energy'' for a periodic system cannot be defined by letting the interacting particles be separated by an infinite distance; the correct definition corresponds to setting its $\bbox k=\bbox 0$ Fourier mode to zero. The implications of this effect for computer simulations are discussed.


title: Global geometry optimization of clusters using a growth strategy optimized by a genetic algorithm

abstract: A new strategy for global geometry optimization of clusters is presented. Important features are a restriction of search space to favorable nearest-neighbor distance ranges, a suitable cluster growth representation with diminished correlations, and easy transferability of the results to larger clusters. The strengths and possible limitations of the method are demonstrated for Si10 using an empirical potential.


title: Studies of an Off-Lattice Model for Protein Folding: Sequence Dependence and Improved Sampling at Finite Temperature

abstract: We study the thermodynamic behavior of a simple off-lattice model for protein folding. The model is two-dimensional and has two different ``amino acids''. Using numerical simulations of all chains containing eight or ten monomers, we examine the sequence dependence at a fixed temperature. It is shown that only a few of the chains exist in unique folded state at this temperature, and the energy level spectra of chains with different types of behavior are compared. Furthermore, we use this model as a testbed for two improved Monte Carlo algorithms. Both algorithms are based on letting some parameter of the model become a dynamical variable; one of the algorithms uses a fluctuating temperature and the other a fluctuating monomer sequence. We find that by these algorithms one gains large factors in efficiency in comparison with conventional methods.


title: Sequence Dependence in an Off-Lattice Model for Protein Folding

abstract: See chem-ph/9505003.


title: On the free energy of ionic hydration

abstract: The hydration free energies of ions exhibit an approximately quadratic dependence on the ionic charge, as predicted by the Born model. We analyze this behavior using second-order perturbation theory. This provides effective methods to calculating free energies from equilibrium computer simulations. The average and the fluctuation of the electrostatic potential at charge sites appear as the first coefficients in a Taylor expansion of the free energy of charging. Combining the data from different charge states allows calculation of free-energy profiles as a function of the ionic charge. The first two Taylor coefficients of the free-energy profiles can be computed accurately from equi- librium simulations; but they are affected by a strong system-size dependence. We apply corrections for these finite-size effects by using Ewald lattice sum- mation and adding the self-interactions consistently. Results are presented for a model ion with methane-like Lennard-Jones parameters in SPC water. We find two very closely quadratic regimes with different parameters for positive and negative ions. We also studied the hydration free energy of potassium, calcium, fluoride, chloride, and bromide ions. We find negative ions to be solvated more strongly compared to positive ions of equal size, in agreement with experimen- tal data. We ascribe this preference of negative ions to their strong interac- tions with water hydrogens, which can penetrate the ionic van der Waals shell without direct energetic penalty in the models used. We also find a positive electrostatic potential at the center of uncharged Lennard-Jones particles in water, which favors negative ions. Regarding finite-system-size effects, we show that even using only 16 water molecules it is possible to calculate accu- rately the hydration free energy of sodium if self-interactions are considered.


title: Thermodynamics of helix-coil transitions studied by multicanonical algorithms

abstract: Thermodynamics of helix-coil transitions in amino-acid homo-oligomers are studied by the recently proposed multicanonical algorithms. Homo-oligomers of length 10 are considered for three characteristic amino acids, alanine (helix former), valine (helix indifferent), and glycine (helix breaker). For alanine other lengths (15 and 20) are also considered in order to examine the length dependence. From one multicanonical production run with completely random initial conformations, we have obtained the lowest-energy conformations and various thermodynamic quantities (average helicity, Zimm-Bragg $s$ and $\sigma$ parameters, free energy differences between helix and coil states, etc.) as functions of temperature. The results confirm the fact that alanine is helix-forming, valine is helix-indifferent, and glycine is helix-breaking.

categories: chem-ph cond-mat hep-lat physics.chem-ph


title: Quantum Confinement Effects in Semiconductor Clusters II

abstract: The band gaps and spectral shifts of CdS, CdSe, CdTe, AlP, GaP, GaAs, and InP semiconductor clusters are calculated from band structure calculations using accurate local and non-local empirical pseudopotentials. The effect of spin-orbit coupling on the band structures is included in the calculations when they are important. The complete set of pseudopotential parameters and full computational details are reported for all these semiconductors. The calculated spectral shifts of zinc-blende and wurtzite CdS, wurtzite CdSe, and zinc-blende InP clusters are in good agreement with experiments over a range of cluster sizes. The effect of crystal structure on the band gaps is small in large clusters but becomes important in small clusters. In the absence of experimental data, our calculations provide reasonable estimates of expected spectral shifts for the other clusters. These results demonstrate that the empirical pseudopotential method yields unique insights into the quantum confinement effects and is a powerful tool for calculating the spectral shifts of semiconductor clusters.


title: Chemistry of Nanoscale Semiconductor Clusters

abstract: The ground state structures of small silicon clusters are determined through exhaustive tight-binding molecular dynamics simulation studies. These simulations revealed that \Si{11} is an icosahedron with one missing cap, \Si{12} is a complete icosahedron, \Si{13} is a surface capped icosahedron, \Si{14} is a 4-4-4 layer structure with two caps, \Si{15} is a 1-5-3-5-1 layer structure, and \Si{16} is a partially closed cage consisting of five-membered rings. The characteristic feature of these clusters is that they are all surface. Smalley and co-workers discovered that chemisorption reactivities of silicon clusters vary over three orders of magnitude as a function of cluster size. In particular, they found that \Si{33}, \Si{39}, and \Si{45} clusters are least reactive towards various reagents compared to their immediate neighbors in size. We provide insights into this observed reactivity pattern through our stuffed fullerene model. This structural model consists of bulk-like core of five atoms surrounded by fullerene-like surface. Reconstruction of the ideal fullerene geometry gives rise to four-fold coordinated crown atoms and $\pi$-bonded dimer pairs. This model yields unique structures for \Si{33}, \Si{39}, and \Si{45} clusters without any dangling bonds and thus explains their lowest reactivity towards chemisorption of closed shell reagents. We also explain why a) these clusters are substantially unreactive compared to bulk surfaces and b) dissociative chemisorption occurs on bulk surfaces while molecular chemisorption occurs on cluster surfaces. Finally, experiments on Si$_x$X$_y$ (X = B, Al, Ga, P, As, AlP, GaAs) are suggested as a means of verifying the proposed model.


title: Growth Pattern of Silicon Clusters

abstract: Tight-binding molecular dynamics simulated annealing technique is employed to search for the ground state geometries of silicon clusters containing 11-17 atoms. These studies revealed that layer formation is the dominant growth pattern in all these clusters. Fullerene-like precursor structures consisting of fused pentagon rings are also observed. The atoms in all these clusters exhibit pronounced preference for residing on the surface.


title: Titrating Polyelectrolytes - Variational Calculations and Monte Carlo Simulations

abstract: Variational methods are used to calculate structural and thermodynamical properties of a titrating polyelectrolyte in a discrete representation. The Coulomb interactions are emulated by harmonic repulsive forces, the force constants being used as variational parameters to minimize the free energy. For the titrating charges, a mean field approach is used. The accuracy is tested against Monte Carlo data for up to 1000 monomers. For an unscreened chain, excellent agreement is obtained for the end-to-end distance and the apparent dissociation constant. With screening, the thermodynamical properties are invariably well described, although the structural agreement deteriorates. A very simple rigid-rod approximation is also considered, giving surprisingly good results for certain properties.


title: The Markoff-Automaton - a New Algorithm for Simulating the Time--Evolution of Large Stochastic Dynamic Systems

abstract: We describe a new algorithm for simulating complex Markoff-processes. We have used a reaction-cell method in order to simulate arbitrary reactions. It can be used for any kind of RDS on arbitrary topologies, including fractal dimensions or configurations not being related to any spatial geometry. The events within a single cell are managed by an event handler which has been implemented independently of the system studied. The method is exact on the Markoff level including the correct treatment of finite numbers of molecules. To demonstrate its properties, we apply it on a very simple reaction-diffusion-systems (RDS). The chemical equations A+A -> inert and A+B -> inert in 1 to 4 dimensions serve as models for systems whose dynamics on an intermediate time scale are governed by fluctuations. We compare our results to the analytic approach by the scaling ansatz. The simulations confirm the exponents of the A+B system within statistical errors, including the logarithmic corrections in the dimension d=2. The method is capable to simulate the crossover from the reaction to diffusion limited regime, which is defined by a crossover time depending on the system size.

categories: chem-ph cond-mat hep-lat mtrl-th physics.chem-ph


title: Spectrum Generating Algebra for X$_{3}$ Molecules

abstract: A new spectrum generating algebra for a unified description of rotations and vibrations in polyatomic molecules is introduced. An application to nonlinear X$_3$ molecules shows that this model (i) incorporates exactly the relevant point group, (ii) provides a complete classification of oblate top states, and (iii) treats properly both degenerate and nondegenerate vibrations.


title: Kinetic and thermodynamic analysis of proteinlike heteropolymers: Monte Carlo histogram technique

abstract: Using Monte Carlo dynamics and the Monte Carlo Histogram Method, the simple three-dimensional 27 monomer lattice copolymer is examined in depth. The thermodynamic properties of various sequences are examined contrasting the behavior of good and poor folding sequences. The good (fast folding) sequences have sharp well-defined thermodynamic transitions while the slow folding sequences have broad ones. We find two independent transitions: a collapse transition to compact states and a folding transition from compact states to the native state. The collapse transition is second order-like, while folding is first order. The system is also studied as a function of the energy parameters. In particular, as the average energetic drive toward compactness is reduced, the two transitions approach each other. At zero average drive, collapse and folding occur almost simultaneously; i.e., the chain collapses directly into the native state. At a specific value of this energy drive the folding temperature falls below the glass point, indicating that the chain is now trapped in local minimum. By varying one parameter in this simple model, we obtain a diverse array of behaviors which may be useful in understanding the different folding properties of various proteins.


title: The hydration free energy of water

abstract: We study the chemical potential of water as a function of charge based on perturbation theory. By calculating the electrostatic-energy fluctuations of two states (fully charged and uncharged) we are able to determine accurate values for the dependence of the chemical potential on charge. We find identical results for the chemical-potential difference of fully charged and uncharged water from overlapping-histogram and acceptance-ratio methods and by smoothly connecting the curves of direct exponential averages. Our results agree with those of Rick and Berne (J. Am. Chem. Soc., 1994, 116, 3949) with respect to both the chemical-potential difference and its dependence on the charge coupling parameter. We observe significant deviations from simple Gaussian-fluctuation statistics. The dependence on the coupling parameter is not quadratic, as would be inferred from linear continuum models of electrostatics.


title: Bifurcation in Rotational Spectra of Nonlinear AB$_2$ Molecules

abstract: A classical microscopic theory of rovibrational motion at high angular momenta in symmetrical non-linear molecules AB$_2$ is derived within the framework of small oscillations near the stationary states of a rotating molecule. The full-dimensional analysis including stretching vibrations has confirmed the existence of the bifurcation predicted previously by means of the rigid-bender model. The formation of fourfold energy clusters has already been experimentally verified for H$_2$Se and it has been demonstrated in fully-dimensional quantum mechanical calculations using the MORBID computer program. We show in the present work that apart from the level clustering, the bifurcation produces physically important effects including molecular symmetry-breaking and a transition from the normal mode to the local mode limit for the stretching vibrations due to rovibrational interaction. The application of the present theory with realistic molecular potentials to the H$_2$Te, H$_2$Se and H$_2$S hydrides results in predictions of the bifurcation points very close to those calculated previously. However for the lighter H$_2$O molecule we find that the bifurcation occurs at higher values of the total angular momentum than obtained in previous estimations. The present work shows it to be very unlikely that the bifurcation in H$_2$O will lead to clustering of energy levels. This result is in agreement with recent variational calculations.


title: A generalized local approximation to the exchange potential

abstract: A new method to obtain a local parameterization for the exchange term in the many--body electronic problem is presented. The approach amounts to the introduction of a coordinate dependent electron effective mass. Numerical results for metallic clusters in the jellium model are compared with other standard methods.


title: Electrostatic potential of a homogeneously charged square and cube in two and three dimensions

abstract: A closed form of the electrostatic potential of a homogeneously charged cube is derived by integration. The exact result is compared with multipole expansions for the exterior and interior of the cube. The electrostatic potential of a homogeneously charged square in two-dimensional electrostatics is also determined.


title: Coherent Control of Isotope Separation in HD+ Photodissociation by Strong Fields

abstract: The photodissociation of the HD+ molecular ion in intense short- pulsed linearly polarized laser fields is studied using a time- dependent wave-packet approach where molecular rotation is fully included. We show that applying a coherent superposition of the fundamental radiation with its second harmonic can lead to asymmetries in the fragment angular distributions, with significant differences between the hydrogen and deuterium distributions in the long wavelength domain where the permanent dipole is most efficient. This effect is used to induce an appreciable isotope separation.


title: Bond Orientational Order, Molecular Motion and Free Energy of High Density DNA Mesophases

abstract: By equilibrating condensed DNA arrays against reservoirs of known osmotic stress and examining them with several structural probes, it has been possible to achieve a detailed thermodynamic and structural characterization of the change between two distinct regions on the liquid crystalline phase digram: a higher-density hexagonally packed region with long-range bond orientational order in the plane perpendicular to the average molecular direction; and a lower-density cholesteric region with fluid-like positional order. X-rays scattering on highly ordered DNA arrays at high density and with the helical axis oriented parallel to the incoming beam showed a six-fold azimuthal modulation of the first order diffraction peak that reflects the macroscopic bond-orientational order. Transition to the less-dense cholesteric phase through osmotically controlled swelling shows the loss of this bond orientational order that had been expected from the change in optical birefringence patterns and that is consistent with a rapid onset of molecular positional disorder. This change in motion was previously inferred from intermolecular force measurements and is now confirmed by $\rm ^{31}P$ NMR. Controlled reversible swelling and compaction under osmotic stress, spanning a range of densities between $\sim 120$ mg/ml to $\sim 600$ mg/ml, allows measurement of the free energy changes throughout each phase and at the phase transition, essential information for theories of liquid-crystalline states.


title: Iterative Solution of the Ornstein-Zernike Equation with Various Closures Using Vector Extrapolation

abstract: The solution of the Ornstein-Zernike equation with various closure approximations is studied. This problem is rewritten as an integral equation that can be solved iteratively on a grid. The convergence of the fixed point iterations is relatively slow. We consider transformations of the sequence of solution vectors using non-linear sequence transformations, so-called vector extrapolation processes. An example is the vector J transformation. The transformed vector sequences turn out to converge considerably faster than the original sequences.


title: A Heuristic Molecular Model of Hydrophobic Interactions

abstract: Hydrophobic interactions provide driving forces for protein folding, membrane formation, and oil-water separation. Motivated by information theory, the poorly understood nonpolar solute interactions in water are investigated. A simple heuristic model of hydrophobic effects in terms of density fluctuations is developed. This model accounts quantitatively for the central hydrophobic phenomena of cavity formation and association of inert gas solutes; it therefore clarifies the underlying physics of hydrophobic effects and permits important applications to conformational equilibria of nonpolar solutes and hydrophobic residues in biopolymers.


title: Mixed Direct-Iterative Methods for Boundary Integral Formulations of Dielectric Solvation Models

abstract: This paper describes a mixed direct-iterative method for boundary integral formulations of dielectric solvation models. We give an example for which a direct solution at thermal accuracy is nontrivial and for which Gauss-Seidel iteration diverges in rare but reproducible cases. This difficulty is analyzed by obtaining the eigenvalues and the spectral radius of the iteration matrix. This establishes that the nonconvergence is due to inaccuracies of the asymptotic approximations for the matrix elements for accidentally close boundary element pairs on different spheres. This difficulty is cured by checking for boundary element pairs closer than the typical spatial extent of the boundary elements and for those pairs performing an `in-line' Monte Carlo integration to evaluate the required matrix elements. This difficulty are not expected and have not been observed when only a direct solution is sought. Finally, we give an example application of these methods to deprotonation of monosilicic acid in water.


title: Boundary Integral Methods for the Poisson Equation of Continuum Dielectric Solvation Models

abstract: This paper tests a dielectric model for variation of hydration free energy with geometry of complex solutes in water. It works out some basic aspects of the theory of boundary integral methods for these problems. One aspect of the algorithmic discussion lays the basis for multigrid methods of solution, methods that are likely to be necessary for similarly accurate numerical solution of these models for much larger solutes. Other aspects of the algorithmic work show how macroscopic surfaces such as solution interfaces and membranes may be incorporated and also show how these methods can be transferred directly to periodic boundary conditions. This dielectric model is found to give interesting and helpful results for the variation in solvation free energy with solute geometry. However, it typically significantly over-stabilizes classic attractive ion-pairing configurations. On the basis of the examples and algorithmic considerations, we make some observations about extension of this continuum model incrementally to reintroduce molecular detail of the solvation structure.


title: A Blocking Technique for Emulating Very Large Polyelectrolytes

abstract: A new Monte Carlo method for computing thermodynamical properties of very large polyelectrolytes is presented. It is based on a renormalization group relating the original polymer to a smaller system, where in addition to the naively rescaled forces, a corrective nearest-neighbor interaction originating from the short distance Coulomb cutoff is introduced. The method is derived for low T but is in the unscreened case valid for all T. Large polymers with N monomers are emulated by Monte Carlo calculations on smaller systems, K=N/Q. The computational gain of the method is Q^3. It is explored with emphasis on room temperature. Results for N=10000 are presented.


title: The Atomic and Electronic Structure of Liquid N- Methylformamide as Determined from Diffraction Experiments

abstract: The structure of liquid N-methylformamide (NMF) has been investigated using synchrotron radiation at 77 and 95 keV. The use of high energy photons has several advantages, in this case especially the large accessible momentum transfer range, the low absorption and the direct comparability with neutron diffraction. The range of momentum transfer covered is 0.6 \AA$^{-1} <$ Q $<$24.0 \AA$^{-1}$. Neutron diffraction data on the same sample in the same momentum transfer range have been published previously. In that study two differently isotope - substituted species were investigated. In order to compare neutron and photon diffraction data properly Reverse Monte Carlo (RMC-) simulations have been performed. Some modifications had to be added to the standard RMC- code introducing different constraints for inter- and intramolecular distances as these distances partly overlap in liquid NMF. RMC- simulations having only the neutron data as input were carried out in order to test the quality of the X-ray data. The photon structure factor calculated from the RMC- configurations is found to agree well with the present experimental data, while it deviates considerably from earlier X-ray work using low energy photons (17 keV). Finally we discuss whether the different interaction mechanisms of neutrons and photons can be used to directly access the electronic structure in the liquid. Evidence is presented that the elastic self scattering part of liquid NMF is changed with respect to the independent atom approximation. This modification can be accounted for by a simple charged atoms model.


title: Influence of surfactants on the structure of titanium oxide gels : experiments and simulations

abstract: We report here on experimental and numerical studies of the influence of surfactants on mineral gel synthesis. The modification of the gel structure when the ratios water-precursor and water-surfactant vary is brought to the fore by fractal dimension measures. A property of {\em polydispersity of the initial hydrolysis} is proposed to explain these results, and is successfuly tested through numerical experiments of three dimensional chemically limited aggregation.


title: Monte-Carlo studies of bosonic van der Waals clusters

abstract: In a previous paper (http://www.phys.uri.edu/people/nightingale/publications.html, chem-ph/9406003) we developed a form of variational trial wave function and applied it to van der Waals clusters: five or less atoms of Ar and Ne modeled by the Lennard-Jones potential. In addition, we tested the trial functions for a hypothetical, light atom resembling Ne but with only half its mass. We did not study atoms such as He with larger de Boer parameters, i.e., systems in which the zero point energy plays a more important role relative to the potential energy. This is the main purpose of the present paper. In fact, we study clusters to the very limit where the zero-point energy destroys the ground state as a bound state. A simple picture of this un-binding transition predicts the power law with which the energy vanishes as the de Boer parameter approaches its critical value and the power of the divergence of the the size of the clusters in this limit. Our numerical results are in agreement with these predictions.


title: Fractal Dimension of Disordered Submonolayers: Determination from He Scattering Data

abstract: We propose a novel method to measure the fractal dimension of a submonolayer metal adatom system grown under conditions of limited diffusivity on a surface. The method is based on measuring the specular peak attenuation of He atoms scattered from the surface, as a function of incidence energy. The (Minkowski) fractal dimension thus obtained is that of contours of constant electron density of the adatom system. Simulation results are presented, based on experimental data. A coverage dependent fractal dimension is found from a two-decade wide scaling regime.

categories: chem-ph atom-ph cond-mat physics.chem-ph


title: Optical Theorem and the Inversion of Cross Section Data for Atom Scattering from Defects on Surfaces

abstract: The information content and properties of the cross section for atom scattering from a defect on a flat surface are investigated. Using the Sudden approximation, a simple expression is obtained that relates the cross section to the underlying atom/defect interaction potential. An approximate inversion formula is given, that determines the shape function of the defect from the scattering data. Another inversion formula approximately determines the potential due to a weak corrugation in the case of substitutional disorder. An Optical Theorem, derived in the framework of the Sudden approximation, plays a central role in deriving the equations that conveniently relate the interaction potential to the cross section. Also essential for the result is the equivalence of the operational definition for the cross section for scattering by a defect, given by Poelsema and Comsa, and the formal definition from quantum scattering theory. This equivalence is established here. The inversion result is applied to determine the shape function of an Ag atom on Pt(111) from scattering data.


title: Evidence for Non-Random Hydrophobicity Structures in Protein Chains

abstract: The question of whether proteins originate from random sequences of amino acids is addressed. A statistical analysis is performed in terms of blocked and random walk values formed by binary hydrophobic assignments of the amino acids along the protein chains. Theoretical expectations of these variables from random distributions of hydrophobicities are compared with those obtained from functional proteins. The results, which are based upon proteins in the SWISS-PROT data base, convincingly show that the amino acid sequences in proteins differ from what is expected from random sequences in a statistical significant way. By performing Fourier transforms on the random walks one obtains additional evidence for non-randomness of the distributions. We have also analyzed results from a synthetic model containing only two amino-acid types, hydrophobic and hydrophilic. With reasonable criteria on good folding properties in terms of thermodynamical and kinetic behavior, sequences that fold well are isolated. Performing the same statistical analysis on the sequences that fold well indicates similar deviations from randomness as for the functional proteins. The deviations from randomness can be interpreted as originating from anticorrelations in terms of an Ising spin model for the hydrophobicities. Our results, which differ from previous investigations using other methods, might have impact on how permissive with respect to sequence specificity the protein folding process is -- only sequences with non-random hydrophobicity distributions fold well. Other distributions give rise to energy landscapes with poor folding properties and hence did not survive the evolution.


title: On the geometry and conformation of starburst dendrimers

abstract: In this work, we propose a geometrical model for the study of conformational properties of a starburst dendrimer with the topology of a truncated Bethe lattice. A convenient embedding of the Bethe lattice in the hyperbolic plane is used to study the architecture of the dendrimer. As results, we find an upper bound for the molecular size and the density profile.


title: Direct subsurface absorption of hydrogen on Pd(111)

abstract: We summarize and discuss some of the available experimental and theoretical data important for understanding the role played by subsurface sites in dissociative chemisorption calculations for the H$_2$/Pd(111) system. Then we use a semi-empirical potential energy surface (PES) to model the interaction of a H$_2$ molecule impinging on a Pd(111) surface. The London-Eyring-Polanyi-Sato (LEPS) construction has been extended to make direct subsurface absorption possible. A 2-dimensional wave packet calculation is used to find qualitative trends in the direct subsurface absorption and to reveal the time scales involved. We suggest that a partial in-plane relaxation occurs for the slowest incoming particles, thus resulting in a higher direct subsurface absorption probability for low energies.


title: Quantum Effects in Barrier Dynamics

abstract: The dynamics near the top of a potential barrier is studied in the temperature region where quantum effects become important. The time evolution of the density matrix of a system that deviates initially from equilibrium in the vicinity of the barrier top but is in local equilibrium away from the barrier top is determined. Explicit results are given for a range of parameters where the nonequilibrium state is not affected by anharmonicities of the barrier potential except for the barrier height. In particular, for a system confined initially to one side of the barrier the relaxation to a quasi--stationary flux state is determined. The associated rate constant is evaluated and the relation to other rate formulas is discussed in detail.


title: He Scattering from Compact Clusters and from Diffusion-Limited Aggregates on Surfaces: Observable Signatures of Structure

abstract: The angular intensity distribution of He beams scattered from compact clusters and from diffusion limited aggregates, epitaxially grown on metal surfaces, is investigated theoretically. The purpose is twofold: to distinguish compact cluster structures from diffusion limited aggregates, and to find observable {\em signatures} that can characterize the compact clusters at the atomic level of detail. To simplify the collision dynamics, the study is carried out in the framework of the sudden approximation, which assumes that momentum changes perpendicular to the surface are large compared with momentum transfer due to surface corrugation. The diffusion limited aggregates on which the scattering calculations were done, were generated by kinetic Monte Carlo simulations. It is demonstrated, by focusing on the example of compact Pt Heptamers, that signatures of structure of compact clusters may indeed be extracted from the scattering distribution. These signatures enable both an experimental distinction between diffusion limited aggregates and compact clusters, and a determination of the cluster structure. The characteristics comprising the signatures are, to varying degrees, the Rainbow, Fraunhofer, specular and constructive interference peaks, all seen in the intensity distribution. It is also shown, how the distribution of adsorbate heights above the metal surface can be obtained by an analysis of the specular peak attenuation. The results contribute to establishing He scattering as a powerful tool in the investigation of surface disorder and epitaxial growth on surfaces, alongside with STM.


title: Dynamics of Proton Transfer in Mesoscopic Clusters

abstract: Proton transfer rates and mechanisms are studied in mesoscopic, liquid-state, molecular clusters. The proton transfer occurs in a proton-ion complex solvated by polar molecules comprising the cluster environment. The rates and mechanisms of the reaction are studied using both adiabatic and non-adiabatic molecular dynamics. For large molecular clusters, the proton-ion complex resides primarily on the surface of the cluster or one layer of solvent molecules inside the surface. The proton transfer occurs as the complex undergoes orientational fluctuations on the cluster surface or penetrates one solvent layer into the cluster leading to solvent configurations that favor the transfer. For smaller clusters the complex resides mostly on the surface of the cluster and proton transfer is observed only when the complex penetrates the cluster and solvent configurations that favor the proton transfer are achieved. Quantitative information on the cluster reaction rate constants is also presented.


title: Binary Assignments of Amino Acids from Pattern Conservation

abstract: We develop a simple optimization procedure for assigning binary values to the amino acids. The binary values are determined by a maximization of the degree of pattern conservation in groups of closely related protein sequences. The maximization is carried out at fixed composition. For compositions approximately corresponding to an equipartition of the residues, the optimal encoding is found to be strongly correlated with hydrophobicity. The stability of the procedure is demonstrated. Our calculations are based upon sequences in the SWISS-PROT database.


title: Molecular dynamics study of solvation effects on acid dissociation in aprotic media

abstract: Acid ionization in aprotic media is studied using Molecular Dynamics techniques. In particular, models for HCl ionization in acetonitrile and dimethylsulfoxide are investigated. The proton is treated quantum mechanically using Feynman path integral methods and the remaining molecules are treated classically. Quantum effects are shown to be essential for the proper treatment of the ionization. The potential of mean force is computed as a function of the ion pair separation and the local solvent structure is examined. The computed dissociation constants in both solvents differ by several orders of magnitude which are in reasonable agreement with experimental results. Solvent separated ion pairs are found to exist in dimethylsulfoxide but not in acetonitrile. Dissociation mechanisms in small clusters are also investigated. Solvent separated ion pairs persist even in aggregates composed of rather few molecules, for instance, as few as thirty molecules. For smaller clusters or for large ion pair separations cluster finite-size effects come into play in a significant fashion.


title: Liquid Crystal Analogue of Abrikosov Vortex Flow in Superconductors

abstract: We extend the correspondence between the Renn-Lubensky Twist-Grain-Boundary-A phase in chiral liquid crystals and the Abrikosov mixed state in superconductors to dynamical aspects. We find that for a TGB sample with free boundaries, an external electric field applied along the helical axis induces a uniform translational motion of the grain boundary system - an analogue of the well-known mixed state flux flow. Likewise, an analogue of the mixed state Nernst effect is found. In much the same way in which the flux flow carries intercore electric fields generating Joule heat in an otherwise dissipation-free system, the grain boundary flow carries along polarized charges, resulting in a finite electric conductivity in a ferroelectric.


title: Lyapunov instability of fluids composed of rigid diatomic molecules

abstract: We study the Lyapunov instability of a two-dimensional fluid composed of rigid diatomic molecules, with two interaction sites each, and interacting with a WCA site-site potential. We compute full spectra of Lyapunov exponents for such a molecular system. These exponents characterize the rate at which neighboring trajectories diverge or converge exponentially in phase space. Quam. These exponents characterize the rate at which neighboring trajectories diverge or converge exponentially in phase space. Qualitative different degrees of freedom -- such as rotation and translation -- affect the Lyapunov spectrum differently. We study this phenomenon by systematically varying the molecular shape and the density. We define and evaluate ``rotation numbers'' measuring the time averaged modulus of the angular velocities for vectors connecting perturbed satellite trajectories with an unperturbed reference trajectory in phase space. For reasons of comparison, various time correlation functions for translation and rotation are computed. The relative dynamics of perturbed trajectories is also studied in certain subspaces of the phase space associated with center-of-mass and orientational molecular motion.


title: Exact and Semiclassical Density Matrix of a Particle Moving in a Barrier Potential with Bound States

abstract: We present a barrier potential with bound states that is exactly solvable and determine the eigenfunctions and eigenvalues of the Hamiltonian. The equilibrium density matrix of a particle moving at temperature T in this nonlinear barrier potential field is determined. The exact density matrix is compared with the result of the path integral approach in the semiclassical approximation. For opaque barriers the simple semiclassical approximation is found to be sufficient at high temperatures while at low temperatures the fluctuation paths may have a caustic depending on temperature and endpoints. Near the caustics the divergence of the simple semiclassical approximation of the density matrix is removed by a nonlinear fluctuation potential. For opaque barriers the improved semiclassical approximation is again in agreement with the exact result. In particular, bound states and the form of resonance states are described accurately by the semiclassical approach.


title: Protein dynamics with off-lattice Monte Carlo moves

abstract: A Monte Carlo method for dynamics simulation of all-atom protein models is introduced, to reach long times not accessible to conventional molecular dynamics. The considered degrees of freedom are the dihedrals at C$_\alpha$-atoms. Two Monte Carlo moves are used: single rotations about torsion axes, and cooperative rotations in windows of amide planes, changing the conformation globally and locally, respectively. For local moves Jacobians are used to obtain an unbiased distribution of dihedrals. A molecular dynamics energy function adapted to the protein model is employed. A polypeptide is folded into native-like structures by local but not by global moves.


title: Scaling and Scale Breaking in Polyelectrolyte

abstract: We consider the thermodynamics of a uniformly charged polyelectrolyte with harmonic bonds. For such a system there is at high temperatures an approximate scaling of global properties like the end-to-end distance and the interaction energy with the chain-length divided by the temperature. This scaling is broken at low temperatures by the ultraviolet divergence of the Coulomb potential. By introducing a renormalization of the strength of the nearest- neighbour interaction the scaling is restored, making possible an efficient blocking method for emulating very large polyelectrolytes using small systems. The high temperature behaviour is well reproduced by the analytical high-$T$ expansions even for fairly low temperatures and system sizes. In addition, results from low-$T$ expansions, where the coefficients have been computed numerically, are presented. These results approximate well the corresponding Monte Carlo results at realistic temperatures. A corresponding analysis of screened chains is performed. The situation here is complicated by the appearance of an additional parameter, the screening length. A window is found in parameter space, where scaling holds for the end-to-end distance. This window corresponds to situations where the range of the potential interpolates between the bond length and the size of the chain. This scaling behaviour, which is verified by Monte Carlo results, is consistent with Flory scaling. Also for the screened chain a blocking approach can be devised, that performs well for low temperatures, whereas the low-$T$ expansion is inaccurate at realistic temperatures.


title: Computational Study of the Structure and Thermodynamic Properties of Ammonium Chloride Clusters Using a Parallel J-Walking Approach

abstract: The thermodynamic and structural properties of (NH$_4$Cl)$_n$ clusters, n=3-10 are studied. Using the method of simulated annealing, the geometries of several isomers for each cluster size are examined. Jump-walking Monte Carlo simulations are then used to compute the constant-volume heat capacity for each cluster size over a wide temperature range. To carry out these simulations a new parallel algorithm is developed using the Parallel Virtual Machine (PVM) software package. Features of the cluster potential energy surfaces, such as energy differences among isomers and rotational barriers of the ammonium ions, are found to play important roles in determining the shape of the heat capacity curves.


title: Correlation Energy Estimators based on M{\o}ller-Plesset Perturbation Theory

abstract: Some methods for the convergence acceleration of the M{\o}ller-Plesset perturbation series for the correlation energy are discussed. The order-by-order summation is less effective than the Feenberg series. The latter is obtained by renormalizing the unperturbed Hamilton operator by a constant factor that is optimized for the third order energy. In the fifth order case, the Feenberg series can be improved by order-dependent optimization of the parameter. Alternatively, one may use Pad{\'e} approximants or a further method based on effective characteristic polynomials to accelerate the convergence of the perturbation series. Numerical evidence is presented that, besides the Feenberg-type approaches, suitable Pad{\'e} approximants, and also the effective second order characteristic polynomial, are excellent tools for correlation energy estimation.

categories: chem-ph atom-ph hep-lat math.NA physics.chem-ph quant-ph


title: A Symmetry Adapted Approach to Molecular Spectroscopy: The Anharmonic Oscillator Symmetry Model

abstract: We apply the Anharmonic Oscillator Symmetry Model to the description of vibrational excitations in ${\cal D}_{3h}$ and ${\cal T}_d$ molecules. A systematic procedure can be used to establish the relation between the algebraic and configuration space formulations, by means of which new interactions are found in the algebraic model, leading to reliable spectroscopic predictions. We illustrate the method for the case of ${\cal D}_{3h}$-triatomic molecules and the ${\cal T}_d$ Be-cluster.


title: On the relation between algebraic and configuration space calculations of molecular vibrations

abstract: The relation between algebraic and traditional calculations of molecular vibrations is investigated. An explicit connection between interactions in configuration space and the corresponding algebraic interactions is established.


title: Bond angle distribution in amorphous germania and silica

abstract: The distribution of Ge-O-Ge and Si-O-Si bond angles alpha in amorphous germania and silica is re-determined on the basis of diffraction experiments. The bond angle alpha joining adjacent tetrahedra is the central parameter of any continuous random network description (CRN) of these glasses. New high energy photon diffraction experiments on amorphous germania (at photon energies of 97 and 149 keV) are presented, covering the momentum transfer 0.6-33.5 AA^{-1}. In photon diffraction experiments on GeO2 the contribution of the OO pairs is very small. To obtain a similar information for amorphous SiO2, high energy photon diffraction experiments have been combined with neutron diffraction data on amorphous silica in order to eliminate the OO- partial structure factor. With this technique it is shown that the Si-O-Si angle distribution is fairly narrow (sigma=7.5 degree) and in fact comparable in width to the Ge-O-Ge angle distribution (sigma=8.3 degree), a result which differs from current opinion. The narrower distribution found in this study are in much better agreement to the determinations based on 29Si-MAS-NMR. Among the various models relating the chemical shift to the bond angle, best agreement is found with those models based on the secant model. Sharp components in the bond angle distribution can be excluded within the reached real space resolution of 0.09 AA.


title: On the energy spectrum of the C60-fullerene anion

abstract: It is shown that the C60-fullerene anion has a hydrogen-like spectrum.


title: From Heisenberg matrix mechanics to EBK quantization: theory and first applications

abstract: Despite the seminal connection between classical multiply-periodic motion and Heisenberg matrix mechanics and the massive amount of work done on the associated problem of semiclassical (EBK) quantization of bound states, we show that there are, nevertheless, a number of previously unexploited aspects of this relationship that bear on the quantum-classical correspondence. In particular, we emphasize a quantum variational principle that implies the classical variational principle for invariant tori. We also expose the more indirect connection between commutation relations and quantization of action variables. With the help of several standard models with one or two degrees of freedom, we then illustrate how the methods of Heisenberg matrix mechanics described in this paper may be used to obtain quantum solutions with a modest increase in effort compared to semiclassical calculations. We also describe and apply a method for obtaining leading quantum corrections to EBK results. Finally, we suggest several new or modified applications of EBK quantization.

categories: chem-ph chao-dyn nlin.CD nucl-th physics.chem-ph quant-ph


title: Motional diminishing of optical activity: a novel method for studying molecular dynamics in liquids and plastic crystals

abstract: Molecular dynamics calculations and optical spectroscopy measurements of weakly active infrared modes are reported. The results are qualitatively understood in terms of the "motional diminishing" of IR lines, a process analogous to the motional narrowing of a nuclear magnetic resonance (NMR) signal. In molecular solids or liquids where the appropriate intramolecular resonances are observable, motional diminishing can be used to study the fluctuations of the intermolecular interactions having time scales of 1psec to 100psec.


title: Van der Waals clusters in the ultra-quantum limit: a Monte Carlo study

abstract: Bosonic van der Waals clusters of sizes three, four and five are studied by diffusion quantum Monte-Carlo techniques. In particular we study the unbinding transition, the ultra-quantum limit where the ground state ceases to exist as a bound state. We discuss the quality of trial wave functions used in the calculations, the critical behavior in the vicinity of the unbinding transition, and simple improvements of the diffusion Monte Carlo algorithm.


title: Boson-realization model for the vibrational spectra of tetrahedral molecules

abstract: An algebraic model of Boson-realization is proposed to study the vibrational spectra of a tetrahedral molecule, where ten sets of boson creation and annihilation operators are used to construct the Hamiltonian with $T_{d}$ symmetry. There are two schemes in our model. The first scheme provides an eight-parameter fit to the published experimental vibrational eigenvalues of methane with a root-mean-square deviation 11.61 $cm^{-1}$. The second scheme, where the bending oscillators are assumed to be harmonic and the interactions between the bending vibrations are neglected, provided a five-parameter fit with a root-mean-square deviation 12.42 $cm^{-1}$


title: Solvent mode participation in the non-radiative relaxation of the hydrated electron

abstract: Non-adiabatic molecular dynamics simulations are used to analyze the role of different solvent degrees of freedom in the non-radiative relaxation of the first excited state of the hydrated electron. The relaxation occurs through a multi-mode coupling between the adiabatic electronic states. The process cannot be described by a single mode promotion model frequently used in the ``large molecule'' limit of gas phase theories. Solvent librations and vibrations, and the H2O asymmetric stretch in particular, are found to be the most effective promotors of the electronic transition. Dissipation of the released energy to the solvent proceeds on two time scales: a fast 10-20 fs heating of the first solvation shell, where most of the energy is accepted by the librational degrees of freedom, and a several hundred femtosecond global reconstruction of t he solvent as the first shell transfers its excess energy to the rest of the molecules. The implications of our use of a semiclassical approximation as the criterion for good promoting and accepting modes are discussed.


title: A Criterion That Determines Fast Folding of Proteins: A Model Study

abstract: We consider the statistical mechanics of a full set of two-dimensional protein-like heteropolymers, whose thermodynamics is characterized by the coil-to-globular ($T_\theta$) and the folding ($T_f$) transition temperatures. For our model, the typical time scale for reaching the unique native conformation is shown to scale as $\tau_f\sim F(M)\exp(\sigma/\sigma_0)$, where $\sigma=1-T_f/T_\theta$, $M$ is the number of residues, and $F(M)$ scales algebraically with $M$. We argue that $T_f$ scales linearly with the inverse of entropy of low energy non-native states, whereas $T_\theta$ is almost independent of it. As $\sigma\rightarrow 0$, non-productive intermediates decrease, and the initial rapid collapse of the protein leads to structures resembling the native state. Based solely on {\it accessible} information, $\sigma$ can be used to predict sequences that fold rapidly.


title: A general algebraic model for molecular vibrational spectroscopy

abstract: We introduce the Anharmonic Oscillator Symmetry Model to describe vibrational excitations in molecular systems exhibiting high degree of symmetry. A systematic procedure is proposed to establish the relation between the algebraic and configuration space formulations, leading to new interactions in the algebraic model. This approach incorporates the full power of group theoretical techniques and provides reliable spectroscopic predictions. We illustrate the method for the case of ${\cal D}_{3h}$-triatomic molecules.


title: Ion-Chain Dynamics in Polymer Electrolytes

abstract: Representing polyether-salt systems by chains of interacting coordination shells, defined by the cation and by its nearest ligands, we derive the interaction potential between closest shells -- the inter-shells potential -- in terms of two-electron polarization effects. Values are presented for monovalent-based crystalline poly(ethylene oxide), PEO, electrolytes. For the eutectic composition $\text{PEO}_{12} \text{EuBr}_3$, the inter-shells energy is evaluated also by relating the empirical value of the nearest-ligands local-field potential with the variation of the $\text{Eu}^{3+}$ concentration. Both methods give the same results.


title: Theory for the Ultrafast Structural Response of optically excited small clusters: Time-dependence of the Ionization Potential

abstract: Combining an electronic theory with molecular dynamics simulations we present results for the ultrafast structural changes in small clusters. We determine the time scale for the change from the linear to a triangular structure after the photodetachment process Ag$_3^- \rightarrow {\rm Ag}_3$. We show that the time-dependent change of the ionization potential reflects in detail the internal degrees of freedom, in particular coherent and incoherent motion, and that it is sensitive to the initial temperature. We compare with experiment and point out the general significance of our results.


title: Symmetry-adapted algebraic description of stretching and bending vibrations of ozone

abstract: The vibrational excitations of ozone, including both bending and stretching vibrations, are studied in the framework of a symmetry-adapted algebraic approach. This method is based on the isomorphism between the U(2) algebra and the one-dimensional Morse oscillator, and the introduction of point group symmetry techniques. The use of symmetry-adapted interactions, which in the harmonic limit have a clear physical interpretation, makes it possible to systematically include higher order terms and anharmonicities. A least-square fit to all published experimental levels (up to ten quanta) of 16 O_3 and 18 O_3 yields a r.m.s. deviation of 2.5 and 1.0 1/cm, respectively.


title: Identification of Amino Acid Sequences with Good Folding Properties in an Off-Lattice Model

abstract: Folding properties of a two-dimensional toy protein model containing only two amino-acid types, hydrophobic and hydrophilic, respectively, are analyzed. An efficient Monte Carlo procedure is employed to ensure that the ground states are found. The thermodynamic properties are found to be strongly sequence dependent in contrast to the kinetic ones. Hence, criteria for good folders are defined entirely in terms of thermodynamic fluctuations. With these criteria sequence patterns that fold well are isolated. For 300 chains with 20 randomly chosen binary residues approximately 10% meet these criteria. Also, an analysis is performed by means of statistical and artificial neural network methods from which it is concluded that the folding properties can be predicted to a certain degree given the binary numbers characterizing the sequences.


title: Extended Two--Particle Green's Functions and Optical Potentials for Two--Particle Scattering by Many--Body Targets

abstract: An extension of the fermionic particle--particle propagator is presented, that possesses similar algebraic properties to the single--particle Green's function. In particular, this extended two--particle Green's function satisfies Dyson's equation and its self energy has the same analytic structure as the the self energy of the single--particle Green's function. For the case of a system interacting by one--particle potentials only, the two--particle self energy takes on a particularly simple form, just like the common self energy does. The new two--particle self energy also serves as a well behaved optical potential for the elastic scattering of a two--particle projectile by a many--body target. Due to its analytic structure, the two--particle self energy avoids divergences that appear with effective potentials derived by other means.


title: On the Stability of Endohedral Rare Gas Fullerenes

abstract: The stability of Ne@C$_{60}$ and He@C$_{60}$ is discussed in the context of a spherical model where the carbon atoms are smeared out into a uniform shell. The electronic properties of the sixty $\pi$ electrons together with those of the central atom are treated in the Thomas-Fermi approximation. Simple electrostatic reasoning elucidates the nature of the radial stability of the complex. A method to include non-spherical corrections is outlined. Possible bonding topologies of the central atom and the C$_{60}$ cage are discussed, as well as the relevance of these topologies to incipient central atom distortions.

categories: chem-ph atom-ph cond-mat mtrl-th physics.chem-ph


title: Charge transfer excitons in C60-dimers and polymers

abstract: Charge-transfer (CT) exciton effects are investigated for the optical absorption spectra of crosslinked C60 systems by using the intermediate exciton theory. We consider the C60-dimers, and the two (and three) molecule systems of the C60-polymers. We use a tight-binding model with long-range Coulomb interactions among electrons, and the model is treated by the Hartree-Fock approximation followed by the single-excitation configuration interaction method. We discuss the variations in the optical spectra by changing the conjugation parameter between molecules. We find that the total CT-component increases in smaller conjugations, and saturates at the intermediate conjugations. It decreases in the large conjugations. We also find that the CT-components of the doped systems are smaller than those of the neutral systems, indicating that the electron-hole distance becomes shorter in the doped C60-polymers.


title: A local anharmonic treatment of vibrations of methane

abstract: The stretching and bending vibrations of methane are studied in a local anharmonic model of molecular vibrations. The use of symmetry-adapted operators reduces the eigenvalue problem to block diagonal form. For the 44 observed energies we obtain a fit with a standard deviation of 0.81 cm$^{-1}$ (and a r.m.s. deviation of 1.16 cm$^{-1}$).


title: Stability and Teller's Theorem: A Thomas-Fermi Theory of Fullerenes

abstract: We study {{\rm C}$_{60}$} with the use of Thomas-Fermi theory. A spherical shell model is invoked to treat the nuclear potential, where the nuclear and core charges are smeared out into a shell of constant surface charge density. The valence electron distribution and the electrostatic potential are efficiently computed by integration of the Thomas-Fermi equation, subject to the shell boundary conditions. Total energy is numerically calculated over a range of shell radii, and the mechanical stability of the model is explored, with attention to the constraints of Teller's theorem. The calculated equilibrium radius of the shell is in good agreement with experiment.


title: On a Thomas-Fermi model of "hollow" atom

abstract: A Thomas-Fermi model of a spherical shell of positive charge is investigated, under various boundary conditions. The electron distribution and the ionization charge are given particular attention.


title: Monte Carlo Optimization of Trial Wave Functions in Quantum Mechanics and Statistical Mechanics

abstract: This review covers applications of quantum Monte Carlo methods to quantum mechanical problems in the study of electronic and atomic structure, as well as applications to statistical mechanical problems both of static and dynamic nature. The common thread in all these applications is optimization of many-parameter trial states, which is done by minimization of the variance of the local or, more generally for arbitrary eigenvalue problems, minimization of the variance of the configurational eigenvalue.

categories: chem-ph atom-ph cond-mat physics.chem-ph quant-ph


title: A Dynamical Theory of Electron Transfer: Crossover from Weak to Strong Electronic Coupling

abstract: We present a real-time path integral theory for the rate of electron transfer reactions. Using graph theoretic techniques, the dynamics is expressed in a formally exact way as a set of integral equations. With a simple approximation for the self-energy, the rate can then be computed analytically to all orders in the electronic coupling matrix element. We present results for the crossover region between weak (nonadiabatic) and strong (adiabatic) electronic coupling and show that this theory provides a rigorous justification for the salient features of the rate expected within conventional electron transfer theory. Nonetheless, we find distinct characteristics of quantum behavior even in the strongly adiabatic limit where classical rate theory is conventionally thought to be applicable. To our knowledge, this theory is the first systematic dynamical treatment of the full crossover region.


title: $\Delta I=4$ and $\Delta I=8$ bifurcations in rotational bands of diatomic molecules

abstract: It is shown that the recently observed $\Delta I=4$ bifurcation seen in superdeformed nuclear bands is also occurring in rotational bands of diatomic molecules. In addition, signs of a $\Delta I=8$ bifurcation, of the same order of magnitude as the $\Delta I=4$ one, are observed both in superdeformed nuclear bands and rotational bands of diatomic molecules.


title: Relation between a Screened Polyelectrolyte and a Field Theory

abstract: There exists a generic relationship between the thermodynamics of a continuous polymer with a generic self-interaction and the two-point function of an interacting field-theory. In addition, the (2N)-point function of the resulting field theory is similarly related to a system of N interacting polymers. In the present paper, this relation is explored for the special case of a polyelectrolyte, characterized by a screened Coulomb pair potential. The corresponding field theory can be recast in a particularly simple form, corresponding to a quantum-mechanical particle, self-interacting via the emission and absorption of a massive scalar field. This is particularly useful in a perturbative treatment: perturbative expansions for the polyelectrolyte can be simply derived from the loop expansion for the related field-theoretical two-point function, for which established computational methods exist.


title: Applications of the Optimized Effective Potential Method of Density Functional Theory to Atomic and Molecular Systems

abstract: Using the optimized effective potential method in conjunction with the semi-analytical approximation due to Krieger, Li and Iafrate, we have performed fully self-consistent exact exchange-only density-functional calculations for diatomic molecules with a fully numerical basis-set-free molecular code. The results are very similar to the ones obtained with the Hartree Fock approach. Furthermore we present results for ground states of positive atomic ions including correlation contributions in the approximation of Colle and Salvetti. It is found that the scheme performs significantly better than conventional Kohn-Sham calculations.


title: A Variational Approach to the Structure and Thermodynamics of Linear Polyelectrolytes with Coulomb and Screened Coulomb Interactions

abstract: A variational approach, based on a discrete representation of the chain, is used to calculate free energy and conformational properties in polyelectrolytes. The true bond and Coulomb potentials are approximated by a trial isotropic harmonic energy containing force constants between {\em all}monomer-pairs as variational parameters. By a judicious choice of representation and the use of incremental matrix inversion, an efficient and fast-convergent iterative algorithm is constructed, that optimizes the free energy. The computational demand scales as $N^3$ rather than $N^4$ as expected in a more naive approach. The method has the additional advantage that in contrast to Monte Carlo calculations the entropy is easily computed. An analysis of the high and low temperature limits is given. Also, the variational formulation is shown to respect the appropriate virial identities.The accuracy of the approximations introduced are tested against Monte Carlo simulations for problem sizes ranging from $N=20$ to 1024. Very good accuracy is obtained for chains with unscreened Coulomb interactions. The addition of salt is described through a screened Coulomb interaction, for which the accuracy in a certain parameter range turns out to be inferior to the unscreened case. The reason is that the harmonic variational Ansatz becomes less efficient with shorter range interactions. As a by-product a very efficient Monte Carlo algorithm was developed for comparisons, providing high statistics data for very large sizes -- 2048 monomers. The Monte Carlo results are also used to examine scaling properties, based on low-$T$ approximations to end-end and monomer-monomer separations. It is argued that the former increases faster than linearly with the number of bonds.


title: A Computational Study of Thirteen-atom Ar-Kr Cluster Heat Capacities

abstract: Heat capacity curves as functions of temperature were calculated using Monte Carlo methods for the series of Ar_{13-n}Kr_n clusters (0 <= n <= 13). The clusters were modeled classically using pairwise additive Lennard-Jones potentials. J-walking (or jump-walking) was used to overcome convergence difficulties due to quasiergodicity present in the solid-liquid transition regions, as well as in the very low temperature regions where heat capacity anomalies arising from permutational isomers were observed. Substantial discrepancies between the J-walking results and the results obtained using standard Metropolis Monte Carlo methods were found. Results obtained using the atom-exchange method, another Monte Carlo variant designed for multi-component systems, were mostly similar to the J-walker results. Quench studies were also done to investigate the clusters' potential energy surfaces; in each case, the lowest energy isomer had an icosahedral-like symmetry typical of homogeneous thirteen-atom rare gas clusters, with an Ar atom being the central atom.


title: An Alternative Conception of Tree-Adjoining Derivation

abstract: The precise formulation of derivation for tree-adjoining grammars has important ramifications for a wide variety of uses of the formalism, from syntactic analysis to semantic interpretation and statistical language modeling. We argue that the definition of tree-adjoining derivation must be reformulated in order to manifest the proper linguistic dependencies in derivations. The particular proposal is both precisely characterizable through a definition of TAG derivations as equivalence classes of ordered derivation trees, and computationally operational, by virtue of a compilation to linear indexed grammars together with an efficient algorithm for recognition and parsing according to the compiled grammar.


title: Lessons from a Restricted Turing Test

abstract: We report on the recent Loebner prize competition inspired by Turing's test of intelligent behavior. The presentation covers the structure of the competition and the outcome of its first instantiation in an actual event, and an analysis of the purpose, design, and appropriateness of such a competition. We argue that the competition has no clear purpose, that its design prevents any useful outcome, and that such a competition is inappropriate given the current level of technology. We then speculate as to suitable alternatives to the Loebner prize.


title: Restricting the Weak-Generative Capacity of Synchronous Tree-Adjoining Grammars

abstract: The formalism of synchronous tree-adjoining grammars, a variant of standard tree-adjoining grammars (TAG), was intended to allow the use of TAGs for language transduction in addition to language specification. In previous work, the definition of the transduction relation defined by a synchronous TAG was given by appeal to an iterative rewriting process. The rewriting definition of derivation is problematic in that it greatly extends the expressivity of the formalism and makes the design of parsing algorithms difficult if not impossible. We introduce a simple, natural definition of synchronous tree-adjoining derivation, based on isomorphisms between standard tree-adjoining derivations, that avoids the expressivity and implementability problems of the original rewriting definition. The decrease in expressivity, which would otherwise make the method unusable, is offset by the incorporation of an alternative definition of standard tree-adjoining derivation, previously proposed for completely separate reasons, thereby making it practical to entertain using the natural definition of synchronous derivation. Nonetheless, some remaining problematic cases call for yet more flexibility in the definition; the isomorphism requirement may have to be relaxed. It remains for future research to tune the exact requirements on the allowable mappings.


title: An Empirically Motivated Reinterpretation of Dependency Grammar

abstract: Dependency grammar is usually interpreted as equivalent to a strict form of X--bar theory that forbids the stacking of nodes of the same bar level (e.g., N' immediately dominating N' with the same head). But adequate accounts of _one_--anaphora and of the semantics of multiple modifiers require such stacking and accordingly argue against dependency grammar. Dependency grammar can be salvaged by reinterpreting its claims about phrase structure, so that modifiers map onto binary--branching X--bar trees rather than ``flat'' ones.


title: Memoization in Constraint Logic Programming

abstract: This paper shows how to apply memoization (caching of subgoals and associated answer substitutions) in a constraint logic programming setting. The research is is motivated by the desire to apply constraint logic programming (CLP) to problems in natural language processing that involve (constraint) interleaving or coroutining, such as GB and HPSG parsing.


title: SPANISH 1992 (S92): corpus-based analysis of present-day Spanish for medical purposes

abstract: S92 research was begun in 1987 to analyze word frequencies in present-day Spanish for making speech pathology evaluation tools. 500 2,000-word samples of children, adolescents and adults' language were input between 1988-1991, calculations done in 1992; statistical and Lewandowski analyses were carried out in 1993.


title: Constraint-Based Categorial Grammar

abstract: We propose a generalization of Categorial Grammar in which lexical categories are defined by means of recursive constraints. In particular, the introduction of relational constraints allows one to capture the effects of (recursive) lexical rules in a computationally attractive manner. We illustrate the linguistic merits of the new approach by showing how it accounts for the syntax of Dutch cross-serial dependencies and the position and scope of adjuncts in such constructions. Delayed evaluation is used to process grammars containing recursive constraints.


title: Principles and Implementation of Deductive Parsing

abstract: We present a system for generating parsers based directly on the metaphor of parsing as deduction. Parsing algorithms can be represented directly as deduction systems, and a single deduction engine can interpret such deduction systems so as to implement the corresponding parser. The method generalizes easily to parsers for augmented phrase structure formalisms, such as definite-clause grammars and other logic grammar formalisms, and has been used for rapid prototyping of parsing algorithms for a variety of formalisms including variants of tree-adjoining grammars, categorial grammars, and lexicalized context-free grammars.


title: A Deductive Account of Quantification in LFG

abstract: The relationship between Lexical-Functional Grammar (LFG) functional structures (f-structures) for sentences and their semantic interpretations can be expressed directly in a fragment of linear logic in a way that explains correctly the constrained interactions between quantifier scope ambiguity and bound anaphora. The use of a deductive framework to account for the compositional properties of quantifying expressions in natural language obviates the need for additional mechanisms, such as Cooper storage, to represent the different scopes that a quantifier might take. Instead, the semantic contribution of a quantifier is recorded as an ordinary logical formula, one whose use in a proof will establish the scope of the quantifier. The properties of linear logic ensure that each quantifier is scoped exactly once. Our analysis of quantifier scope can be seen as a recasting of Pereira's analysis (Pereira, 1991), which was expressed in higher-order intuitionistic logic. But our use of LFG and linear logic provides a much more direct and computationally more flexible interpretation mechanism for at least the same range of phenomena. We have developed a preliminary Prolog implementation of the linear deductions described in this work.


title: Intensional Verbs Without Type-Raising or Lexical Ambiguity

abstract: We present an analysis of the semantic interpretation of intensional verbs such as seek that allows them to take direct objects of either individual or quantifier type, producing both de dicto and de re readings in the quantifier case, all without needing to stipulate type-raising or quantifying-in rules. This simple account follows directly from our use of logical deduction in linear logic to express the relationship between syntactic structures and meanings. While our analysis resembles current categorial approaches in important ways, it differs from them in allowing the greater type flexibility of categorial semantics while maintaining a precise connection to syntax. As a result, we are able to provide derivations for certain readings of sentences with intensional verbs and complex direct objects that are not derivable in current purely categorial accounts of the syntax-semantics interface. The analysis forms a part of our ongoing work on semantic interpretation within the framework of Lexical-Functional Grammar.


title: Adjuncts and the Processing of Lexical Rules

abstract: The standard HPSG analysis of Germanic verb clusters can not explain the observed narrow-scope readings of adjuncts in such verb clusters. We present an extension of the HPSG analysis that accounts for the systematic ambiguity of the scope of adjuncts in verb cluster constructions, by treating adjuncts as members of the subcat list. The extension uses powerful recursive lexical rules, implemented as complex constraints. We show how `delayed evaluation' techniques from constraint-logic programming can be used to process such lexical rules.


title: Similarity-Based Estimation of Word Cooccurrence Probabilities

abstract: In many applications of natural language processing it is necessary to determine the likelihood of a given word combination. For example, a speech recognizer may need to determine which of the two word combinations ``eat a peach'' and ``eat a beach'' is more likely. Statistical NLP methods determine the likelihood of a word combination according to its frequency in a training corpus. However, the nature of language is such that many word combinations are infrequent and do not occur in a given corpus. In this work we propose a method for estimating the probability of such previously unseen word combinations using available information on ``most similar'' words. We describe a probabilistic word association model based on distributional word similarity, and apply it to improving probability estimates for unseen word bigrams in a variant of Katz's back-off model. The similarity-based method yields a 20% perplexity improvement in the prediction of unseen bigrams and statistically significant reductions in speech-recognition error.


title: Temporal Relations: Reference or Discourse Coherence?

abstract: The temporal relations that hold between events described by successive utterances are often left implicit or underspecified. We address the role of two phenomena with respect to the recovery of these relations: (1) the referential properties of tense, and (2) the role of temporal constraints imposed by coherence relations. We account for several facets of the identification of temporal relations through an integration of these.


title: Some Bibliographical References on Intonation and Intonational Meaning

abstract: A by-no-means-complete collection of references for those interested in intonational meaning, with other miscellaneous references on intonation included. Additional references are welcome, and should be sent to julia@research.att.com.


title: Syntactic-Head-Driven Generation

abstract: The previously proposed semantic-head-driven generation methods run into problems if none of the daughter constituents in the syntacto-semantic rule schemata of a grammar fits the definition of a semantic head given in Shieber et al. 1990. This is the case for the semantic analysis rules of certain constraint-based semantic representations, e.g. Underspecified Discourse Representation Structures (UDRSs) (Frank/Reyle 1992). Since head-driven generation in general has its merits, we simply return to a syntactic definition of `head' and demonstrate the feasibility of syntactic-head-driven generation. In addition to its generality, a syntactic-head-driven algorithm provides a basis for a logically well-defined treatment of the movement of (syntactic) heads, for which only ad-hoc solutions existed, so far.


title: Pearl: A Probabilistic Chart Parser

abstract: This paper describes a natural language parsing algorithm for unrestricted text which uses a probability-based scoring function to select the "best" parse of a sentence. The parser, Pearl, is a time-asynchronous bottom-up chart parser with Earley-type top-down prediction which pursues the highest-scoring theory in the chart, where the score of a theory represents the extent to which the context of the sentence predicts that interpretation. This parser differs from previous attempts at stochastic parsers in that it uses a richer form of conditional probabilities based on context to predict likelihood. Pearl also provides a framework for incorporating the results of previous work in part-of-speech assignment, unknown word models, and other probabilistic models of linguistic features into one parsing tool, interleaving these techniques instead of using the traditional pipeline architecture. In preliminary tests, Pearl has been successful at resolving part-of-speech and word (in speech processing) ambiguity, determining categories for unknown words, and selecting correct parses first using a very loosely fitting covering grammar.


title: Efficiency, Robustness, and Accuracy in Picky Chart Parsing

abstract: This paper describes Picky, a probabilistic agenda-based chart parsing algorithm which uses a technique called {\em probabilistic prediction} to predict which grammar rules are likely to lead to an acceptable parse of the input. Using a suboptimal search method, Picky significantly reduces the number of edges produced by CKY-like chart parsing algorithms, while maintaining the robustness of pure bottom-up parsers and the accuracy of existing probabilistic parsers. Experiments using Picky demonstrate how probabilistic modelling can impact upon the efficiency, robustness and accuracy of a parser.


title: Towards History-based Grammars: Using Richer Models for Probabilistic Parsing

abstract: We describe a generative probabilistic model of natural language, which we call HBG, that takes advantage of detailed linguistic information to resolve ambiguity. HBG incorporates lexical, syntactic, semantic, and structural information from the parse tree into the disambiguation process in a novel way. We use a corpus of bracketed sentences, called a Treebank, in combination with decision tree building to tease out the relevant aspects of a parse tree that will determine the correct parse of a sentence. This stands in contrast to the usual approach of further grammar tailoring via the usual linguistic introspection in the hope of generating the correct parse. In head-to-head tests against one of the best existing robust probabilistic parsing models, which we call P-CFG, the HBG model significantly outperforms P-CFG, increasing the parsing accuracy rate from 60% to 75%, a 37% reduction in error.


title: A Stochastic Finite-State Word-Segmentation Algorithm for Chinese

abstract: We present a stochastic finite-state model for segmenting Chinese text into dictionary entries and productively derived words, and providing pronunciations for these words; the method incorporates a class-based model in its treatment of personal names. We also evaluate the system's performance, taking into account the fact that people often do not agree on a single segmentation.


title: Natural Language Parsing as Statistical Pattern Recognition

abstract: Traditional natural language parsers are based on rewrite rule systems developed in an arduous, time-consuming manner by grammarians. A majority of the grammarian's efforts are devoted to the disambiguation process, first hypothesizing rules which dictate constituent categories and relationships among words in ambiguous sentences, and then seeking exceptions and corrections to these rules. In this work, I propose an automatic method for acquiring a statistical parser from a set of parsed sentences which takes advantage of some initial linguistic input, but avoids the pitfalls of the iterative and seemingly endless grammar development process. Based on distributionally-derived and linguistically-based features of language, this parser acquires a set of statistical decision trees which assign a probability distribution on the space of parse trees given the input sentence. These decision trees take advantage of significant amount of contextual information, potentially including all of the lexical information in the sentence, to produce highly accurate statistical models of the disambiguation process. By basing the disambiguation criteria selection on entropy reduction rather than human intuition, this parser development method is able to consider more sentences than a human grammarian can when making individual disambiguation rules. In experiments between a parser, acquired using this statistical framework, and a grammarian's rule-based parser, developed over a ten-year period, both using the same training material and test sentences, the decision tree parser significantly outperformed the grammar-based parser on the accuracy measure which the grammarian was trying to maximize, achieving an accuracy of 78% compared to the grammar-based parser's 69%.


title: Common Topics and Coherent Situations: Interpreting Ellipsis in the Context of Discourse Inference

abstract: It is claimed that a variety of facts concerning ellipsis, event reference, and interclausal coherence can be explained by two features of the linguistic form in question: (1) whether the form leaves behind an empty constituent in the syntax, and (2) whether the form is anaphoric in the semantics. It is proposed that these features interact with one of two types of discourse inference, namely {\it Common Topic} inference and {\it Coherent Situation} inference. The differing ways in which these types of inference utilize syntactic and semantic representations predicts phenomena for which it is otherwise difficult to account.