Showing papers in "arXiv: Chemical Physics in 2004"

Development of Novel Density Functionals for Thermochemical Kinetics

Abstract: A new density functional theory (DFT) exchange-correlation functional for the exploration of reaction mechanisms is proposed. This new functional, denoted BMK (Boese-Martin for Kinetics), has an accuracy in the 2 kcal/mol range for transition state barriers but, unlike previous attempts at such a functional, this improved accuracy does not come at the expense of equilibrium properties. This makes it a general-purpose functional whose domain of applicability has been extended to transition states, rather than a specialized functional for kinetics. The improvement in BMK rests on the inclusion of the kinetic energy density together with a large value of the exact exchange mixing coefficient. For this functional, the kinetic energy density appears to correct `back' the excess exact exchange mixing for ground-state properties, possibly simulating variable exchange.

Current-Voltage Relations for Electrochemical Thin Films

Abstract: The dc response of an electrochemical thin film, such as the separator in a micro-battery, is analyzed by solving the Poisson-Nernst-Planck equations, subject to boundary conditions appropriate for an electrolytic/galvanic cell. The model system consists of a binary electrolyte between parallel-plate electrodes, each possessing a compact Stern layer, which mediates Faradaic reactions with nonlinear Butler-Volmer kinetics. Analytical results are obtained by matched asymptotic expansions in the limit of thin double layers and compared with full numerical solutions. The analysis shows that (i) decreasing the system size relative to the Debye screening length decreases the voltage of the cell and allows currents higher than the classical diffusion-limited current; (ii) finite reaction rates lead to the important possibility of a reaction-limited current; (iii) the Stern-layer capacitance is critical for allowing the cell to achieve currents above the reaction-limited current; and (iv) all polarographic (current-voltage) curves tend to the same limit as reaction kinetics become fast. Dimensional analysis, however, shows that ``fast'' reactions tend to become ``slow'' with decreasing system size, so the nonlinear effects of surface polarization may dominate the dc response of thin films.

Stacked clusters of polycyclic aromatic hydrocarbon molecules

Abstract: Clusters of polycyclic aromatic hydrocarbon (PAH) molecules are modelled using explicit all-atom potentials using a rigid body approximation. The PAH's considered range from pyrene (C10H8) to circumcoronene (C54H18), and clusters containing between 2 and 32 molecules are investigated. In addition to the usual repulsion-dispersion interactions, electrostatic point-charge interactions are incorporated, as obtained from density functional theory calculations. The general electrostatic distribution in neutral or singly charged PAH's is reproduced well using a fluctuating charges analysis, which provides an adequate description of the multipolar distribution. Global optimization is performed using a variety of methods, including basin-hopping and parallel tempering Monte Carlo. We find evidence that stacking the PAH molecules generally yields the most stable motif. A structural transition between one-dimensional stacks and three-dimensional shapes built from mutiple stacks is observed at larger sizes, and the threshold for this transition increases with the size of the monomer. Larger aggregates seem to evolve toward the packing observed for benzene in bulk.Difficulties met in optimizing these clusters are analysed in terms of the strong anisotropy of the molecules. We also discuss segregation in heterogeneous clusters and vibrational properties in the context of astrophysical observations.

Hole Scavenging and Photo-Stimulated Recombination of Electron-Hole Pairs in Aqueous TiO2 Nanoparticles

Abstract: It is shown that 532 nm and 1064 nm laser photoexcitation of trapped electrons generated by 355 nm photolysis of aqueous titania (TiO2) nanoparticles causes rapid photobleaching of their absorbance band in the visible and near IR. This photobleaching occurs within the duration of the laser pulse (3 ns FWHM); it is caused by photoinduced electron detrapping followed by rapid recombination of the resulting free electron and a trapped hole. The quantum yield for the electron photobleaching is ca. 0.28 for 532 nm and ca. 0.024 for 1064 nm photoexcitation. Complete separation of the spectral contributions from trapped electron and hole is demonstrated using glycerol as a selective hole scavenger. When glycerol is added to the solution, some light-absorbing holes are scavenged promptly within the duration of the 355 nm photoexcitation pulse, some are scavenged at a slower rate over the first 200 ns after the 355 nm pulse, and the rest are not scavenged, even at high concentration of the scavenger (> 10 vol %). A reaction with chemi- and physi- sorbed glycerol would account for the prompt and the slow hole decay, respectively. The implications of these results are discussed.

On the efficiency of exchange in parallel tempering Monte Carlo simulations

Abstract: We introduce the concept of effective fraction, defined as the expected probability that a configuration from the lowest index replica successfully reaches the highest index replica during a replica exchange Monte Carlo simulation. We then argue that the effective fraction represents an adequate measure of the quality of the sampling technique, as far as swapping is concerned. Under the hypothesis that the correlation between successive exchanges is negligible, we propose a technique for the computation of the effective fraction, a technique that relies solely on the values of the acceptance probabilities obtained at the end of the simulation. The effective fraction is then utilized for the study of the efficiency of a popular swapping scheme in the context of parallel tempering in the canonical ensemble. For large dimensional oscillators, we show that the swapping probability that minimizes the computational effort is 38.74%. By studying the parallel tempering swapping efficiency for a 13-atom Lennard-Jones cluster, we argue that the value of 38.74% remains roughly the optimal probability for most systems with continuous distributions that are likely to be encountered in practice.

Electron Photodetachment from Aqueous Anions. I. Quantum Yields for Generation of Hydrated Electron by 193 and 248 nm Laser Photoexcitation of Miscellaneous Inorganic Anions

Abstract: Time resolved transient absorption spectroscopy has been used to determine quantum yields for electron photodetachment in 193 nm and (where possible) 248 nm laser excitation of miscellaneous aqueous anions, including hexacyanoferrate(II), sulfate, sulfite, halide anions (Cl-, Br-, and I-), pseudohalide anions (OH-, HS-, CNS-), and several common inorganic anions for which no quantum yields have been reported heretofore: SO3=, NO2-, NO3-, ClO3- and ClO4-. Molar extinction coefficients for these anions and photoproducts of electron detachment from these anions at the excitation wavelengths were also determined. These results are discussed in the context of recent ultrafast kinetic studies and compared with the previous data obtained by product analyses. We suggest using electron photodetachment from the aqueous halide and pseudohalide anions as actinometric standard for time-resolved studies of aqueous photosystems in the UV.

Path-integral Monte-Carlo simulations for electronic dynamics on molecular chains: I. Sequential hopping and super exchange

Abstract: An improved real-time quantum Monte Carlo procedure is presented and applied to describe the electronic transfer dynamics along molecular chains. The model consists of discrete electronic sites coupled to a thermal environment which is integrated out exactly within the path integral formulation. The approach is numerically exact and its results reduce to known analytical findings (Marcus theory, golden rule) in proper limits. Special attention is paid to the role of superexchange and sequential hopping at lower temperatures in symmetric donor-bridge-acceptor systems. In contrast to previous approximate studies, superexchange turns out to play a significant role only for extremely high lying bridges where the transfer is basically frozen or for extremely low temperatures where for weaker dissipation a description in terms of rate constants is no longer feasible. For bridges with increasing length an algebraic decrease of the yield is found for short as well as for longer bridges. The approach can be extended to electronic systems with more complicated topologies including impurities and in presence of external time dependent forces.

A Liouville equation for systems which exchange particles with reservoirs: transport through a nano-device

Abstract: A Redfield-like Liouville equation for an open system that couples to one or more leads and exchanges particles with them is derived. The equation is presented for a general case. A case study of time-dependent transport through a single quantum level for varying electrostatic and chemical potentials in the leads is presented. For the case of varying electrostatic potentials the proposed equation yields, for the model study, the results of an exact solution.

A general expression for bimolecular association rates with orientational constraints

Abstract: We present a general expression for the association rate for partially diffusion-controlled reactions between spherical molecules with an asymmetric reactive patch on each surface. Reaction can occur only if the two patches are in contact and properly aligned to within specified angular tolerances. This extends and generalizes previous approaches that considered only axially symmetric patches; the earlier solutions are shown to be limiting cases of our general expression. Previous numerical results on the rate of protein-protein association with high steric specificity are in very good agreement with the value computed from our analytic expression. Using the new expression, we investigate the influence of orientational constraints on the rate constant. We find that for angular constraints of ~5-15 degrees, a typical range for example in the case of protein-protein interactions, the reaction rate is about 2 to 3 orders of magnitude higher than expected from a simple geometric model.

Electron Photodetachment from Aqueous Anions. II. Ionic Strength Effect on Geminate Recombination Dynamics and Quantum Yield for Hydrated Electron

Abstract: In concentrated solutions of NaClO4 and Na2SO4, the quantum yield for free electron generated by detachment from photoexcited anions (such as I-, OH-, ClO^4-, and [SO3]^2-) linearly decreases by 6-12% per 1 M ionic strength. In 9 M sodium perchlorate solution, this quantum yield decreases by roughly an order of magnitude. Ultrafast kinetic studies of 200 nm photon induced electron detachment from Br-, HO- and [SO3]^2- suggest that the prompt yield of thermalized electron does not change in these solutions; rather, the ionic strength effect originates in more efficient recombination of geminate pairs. Within the framework of the recently proposed mean force potential (MFP) model of charge separation dynamics in such photosystems, the observed changes are interpreted as an increase in the short-range attractive potential between the geminate partners. Association of sodium cation(s) with the electron and the parent anion is suggested as the most likely cause for the observed modification of the MFP. Electron thermalization kinetics suggest that the cation associated with the parent anion (by ion pairing and/or ionic atmosphere interaction) is passed to the detached electron in the course of the photoreaction. The precise atomic-level mechanism for the ionic strength effect is presently unclear; any further advance is likely to require the development of an adequate quantum molecular dynamics model.

Radiation Chemistry of Organic Liquids: Saturated Hydrocarbons

Abstract: In this review (124 refs), several problems in radiolysis of saturated hydrocarbons are examined. Special attention is paid to the chemistry of radical cations, high-mobility holes, excited state and spur dynamics, magnetic field and spin effects, and optically detected magnetic resonance spectroscopy.

Photo-Stimulated Electron Detrapping and the Two-State Model for Electron Transport in Nonpolar Liquids

Abstract: In common nonpolar liquids, such as saturated hydrocarbons, a dynamic equilibrium between trapped (localized) and quasifree (extended) states has been postulated for the excess electron (the two-state model). Using time-resolved dc conductivity, the effect of 1064 nm laser photoexcitation of trapped electrons on the charge transport has been observed in liquid n-hexane and methylcyclohexane. The light promotes the electron from the trap into the conduction band of the liquid, instantaneously increasing the conductivity by orders of magnitude. From the analysis of the two-pulse, two-color photoconductivity data, the residence time of the electrons in traps has been estimated as ca. 8.4 ps for n-hexane and ca. 13 ps for methylcyclohexane (at 295 K). The rate of detrapping decreases at lower temperature with an activation energy of ca. 200 meV (280-320 K); the lifetime-mobility product for quasifree electrons scales linearly with the temperature. We suggest that the properties of trapped electrons in hydrocarbon liquids can be well accounted for using the simple electron bubble (Wigner-Seiz spherical well) model. The estimated localization time of the quasifree electron is 20-50 fs; both time estimates are in good agreement with the "quasiballistic" model. This localization time is significantly lower than the value of ca. 300 fs obtained using time-domain terahertz (THz) spectroscopy for the same system [E. Knoesel et al., J. Chem. Phys. 121, 394 (2004)]. We suggest that the THz signal originates from the oscillations of electron bubbles rather than the free-electron plasma; vibrations of these bubbles may be responsible for the deviations from the Drude behavior observed below 0.4 THz. Various implications of these results are discussed.

The Clausius inequality: implications for nonequilibrium thermodynamic steady states with NEMD corroboration

Abstract: The Clausius inequality for closed systems was deduced from the Riemann integration of closed Carnot cycle loops for irreversible transitions. The corresponding inequalities for the recently developed open system Carnot cycles are derived here and their properties indicate that no new nonequilibrium entropy results from them as has been proposed over the decades. A simulation using a unique dimer molecule with a hysteresis potential - which fulfills all the dynamical laws - is used to show breakdown of the local equilibrium hypothesis. The cardinal Curie principle is shown inadequate as a paradigm to select force-flux relationships

Optimal choice of dividing surface for the computation of quantum reaction rates

Abstract: We consider the calculation of quantum mechanical rate constants for chemical reactions via algorithms that utilize short-time values of the symmetrized flux-flux correlation function. We argue that the dividing surface that makes optimal use of the short-time quantum information is the surface that minimizes the value at the origin of the symmetrized flux-flux correlation function. We also demonstrate that, in the classical limit, this quantum variational criterion produces the same dividing surface as Wigner's variational principle. Finally, we argue that the quantum variational criterion behaves in a nearly optimal fashion with respect to the minimization of the extent of re-crossing flux.

An Application of Transfer Hamiltonian Quantum Mechanics to Multi-Scale Modeling

Abstract: In quantum/classical (QM/CM) partitioning methods for multi-scale modeling, one is often forced to introduce uncontrolled phenomenological effects of the environment (CM) in the quantum (QM) domain as ab initio quantum calculations are computationally too intensive to be applied to the whole sample. We propose a method, in which two qualitatively different components of the information about the state of the CM region are incorporated into the QM calculations. First, pseudoatoms constructed to describe the chemistry of the nearest neighbor exchange interactions replace the atoms at the boundary of the CM and the QM regions. Second, the remaining effect of the CM bulk environment due to long-range Coulombic interactions is modeled in terms of dipoles. We have tested this partitioning method in a silica nanorod and a 3-membered silica ring for which ab initio quantum data for the whole system is available to assess the quality of the proposed partitioning method.

Analysis of an algebraic model for the chromophore vibrations of CF$_3$CHFI

Abstract: We extract the dynamics implicit in an algebraic fitted model Hamiltonian for the hydrogen chromophore's vibrational motion in the molecule $CF_3CHFI$. The original model has 4 degrees of freedom, three positions and one representing interbond couplings. A conserved polyad allows the reduction to 3 degrees of freedom. For most quantum states we can identify the underlying motion that when quantized gives the said state. Most of the classifications, identifications and assignments are done by visual inspection of the already available wave function semiclassically transformed from the number representation to a representation on the reduced dimension toroidal configuration space corresponding to the classical action and angle variables. The concentration of the wave function density to lower dimensional subsets centered on idealized simple lower dimensional organizing structures and the behavior of the phase along such organizing centers already reveals the atomic motion. Extremely little computational work is needed.

Calculation of the Critical Temperature for 2- and 3-Dimensional Ising Models and for 2-Dimensional Potts Models Using the Transfer Matrix Method

Abstract: A new graphical method is developed to calculate the critical temperature of 2- and 3-dimensional Ising models as well as that of the 2-dimensional Potts models. This method is based on the transfer matrix method and using the limited lattice for the calculation. The reduced internal energy per site has been accurately calculated for different 2-D Ising and Potts models using different size-limited lattices. All calculated energies intersect at a single point when plotted versus the reduced temperature. The reduced temperature at the intersection is 0.4407, 0.2746, and 0.6585 for the square, triangular, and honeycombs Ising lattices and 1.0050, 0.6309, and 1.4848 for the square, triangular, and honeycombs Potts lattices, respectively. These values are exactly the same as the critical temperatures reported in the literature, except for the honeycomb Potts lattice. For the two-dimensional Ising model, we have shown that the existence of such an intersection point is due to the duality relation. The method is then extended to the simple cubic Ising model, in which the intersection point is found to be dependent on the lattice sizes. We have found a linear relation between the lattice size and the intersection point. This relation is used to obtain the critical temperature of the unlimited simple cubic lattice. The obtained result, 0.221(2), is in a good agreement with the accurate value of 0.22165 reported by others.

Chemical Equilibrium as Balance of the Thermodynamic Forces

Abstract: The article sets forth comprehensive basics of thermodynamics of chemical equilibrium as balance of the thermodynamic forces. Based on the linear equations of irreversible thermodynamics, De Donder definition of the thermodynamic force, and Le Chatelier's principle, new thermodynamics of chemical equilibrium offers an explicit account for multiple chemical interactions within the system. Basic relations between energetic characteristics of chemical transformations and reaction extents are based on the idea of chemical equilibrium as balance between internal and external thermodynamic forces, which is presented in the form of a logistic equation, containing only one new parameter. Solutions to the basic equation define the domain of states of the chemical system, from true equilibrium to true chaos. The new theory is derived exclusively from the currently recognized ideas and covers equilibrium thermodynamics as well as non-equilibrium thermodynamics in a unique concept.

Quantum mechanical considerations on algebraic structure of central molecular chirality

Abstract: The chiral algebra of tetrahedral molecules, derived from Fischer projections, is discussed in the framework of quantum mechanics. A quantum chiral algebra is obtained whose operators, acting as rotations or inversions, commute with the Hamiltonian of the system. It is shown that energy and chirality eigenstates are strictly related through the Heisenberg relations, while chirality operators conserve parity eigenstates.

Structural, Vibrational and Mechanical Studies of Hydroxyapatite produced by wet-chemical methods

Abstract: Hydroxyapatite samples were produced by two different wet-chemical methods, and characterized by x-ray diffraction, infrared and compression strength measurements. The x-ray diffraction measurements were simulated using the Rietveld method, and structural data as lattice parameters and average crystallite size were obtained. The infrared spectra showed the presence of CO$_3^{2-}$ ions in all samples, indicating a contamination by these ions. By mixing samples produced by both methods, a bioceramic was obtained and, after sintering, samples with very high compression strengths (26--30 MPa) were obtained.

Accounting for correlations with core electrons by means of the generalized RECP: Atoms Hg and Pb and their compounds

Abstract: A way to account for atomic-type correlations between the chemically active (valence) and innermore (core) electrons in the framework of the generalized relativistic effective core potential (GRECP) method is suggested. The "correlated" GRECPs (CGRECPs) are generated for the Hg and Pb atoms. Only correlations for the external twelve and four electrons of them, correspondingly, should be treated explicitly in the subsequent calculations with these CGRECPs whereas the innermore electrons are excluded from the calculations. Results of atomic calculations with the correlated and earlier GRECP versions are compared with the corresponding all-electron Dirac-Coulomb values. Calculations with the above GRECPs and CGRECPs are also carried out for the lowest-lying states of the HgH molecule and its cation and for the ground state of the PbO molecule as compared to earlier calculations and experimental data. The accuracy for the vibrational frequencies is increased up to an order of magnitude and the errors for the bond lengths (rotational constants) are decreased in about two times when the correlated GRECPs are applied instead of earlier GRECP versions employing the same innercore-outercore-valence partitioning.

Negative capacitance and related instabilities in theoretical models of the electric double layer and membrane capacitors

Abstract: Various models leading to predictions of negative capacitance, C, are briefly reviewed. Their relation to the nature of electric control is discussed. We reconfirm that the calculated double layer capacitance can be negative under S-control - an artificial construct that requires uniform distribution of the electrode surface charge density, S. It is shown that the combined relaxation of the ionic and electronic contributions can result in C<0 even for the local statistical ionic models with strictly positive diffuse layer capacitance. In reality, however, only the total charge q (or the average surface charge density S) can be experimentally fixed in isolated cell studies (q-control). For those S where C becomes negative under S-control, the transition to q-control (i.e. relaxing the lateral charge density distribution, fixing its mean value to S) leads to instability of the uniform distribution and a transition to a non-uniform phase. As an illustration, a "membrane capacitor" model is discussed. This exactly solvable model, allowing for both uniform and inhomogeneous relaxation of the electrical double layer, helps to demonstrate both the onset and some important features of the instability. Perspectives for further development are discussed briefly.

Extraction of Thermodynamic Data from Ternary Diffusion Coefficients of Lysozyme Chloride in Water and Aqueous Na$_2$SO$_4$

Abstract: This paper presents, for ternary lysozyme-Na$_2$SO$_4$-water system, the thermodynamic data extracted from the measured values of four ternary diffusion coefficients and the Onsager reciprocal relations. The calculation for derivatives of solute chemical potentials with respect to solute molar concentrations was made using the method presented in \cite{1}. This method is applicable to systems in which the molar concentration of one solute is very small compared to that of the other, like in our case. The approach is illustrated for the lysozyme chloride-Na$_2$SO$_4$-water system at 25$^o$ C, pH 4.5 and at 0.6 mM (8.6 mg/mL) lysozyme chloride and 0.1, 0.25, 0.5, 0.65, and 0.8 M Na$_2$SO$_4$ concentrations. The calculated solute chemical potential derivatives were used to compute the protein cation charge approximately. We also compute the diffusion Onsager coefficients $(L_{ij})_o$ for each composition at pH 4.5.

Numerical Evaluation of Cuspoid and Bessoid Oscillating Integrals for Applications in Chemical Physics

Abstract: Oscillating integrals often arise in the theoretical description of phenomena in chemical physics, in particular in atomic and molecular collisions, and in spectroscopy A computer code for the numerical evaluation of the oscillatory cuspoid canonical integrals and their first-order partial derivatives is described The code uses a novel adaptive contour algorithm, which chooses a contour in the complex plane that avoids the violent oscillatory and exponential natures of the integrand and modifies its choice as necessary Applications are made to the swallowtail canonical integral and to a bessoid integral

Core excitation in Ozone localized to one of two symmetry-equivalent chemical bonds - molecular alignment through vibronic coupling

Abstract: Core excitation from terminal oxygen O$_T$ in O$_3$ is shown to be an excitation from a localized core orbital to a localized valence orbital. The valence orbital is localized to one of the two equivalent chemical bonds. We experimentally demonstrate this with the Auger Doppler effect which is observable when O$_3$ is core-excited to the highly dissociative O$_{T}$1s$^{-1}$7a$_1^1$ state. Auger electrons emitted from the atomic oxygen fragment carry information about the molecular orientation relative to the electromagnetic field vector at the moment of excitation. The data together with analytical functions for the electron-peak profiles give clear evidence that the preferred molecular orientation for excitation only depends on the orientation of one bond, not on the total molecular orientation. The localization of the valence orbital "7a$_1$" is caused by mixing of the valence orbital "5b$_2$" through vibronic coupling of anti-symmetric stretching mode with b$_2$-symmetry. To the best of our knowledge, it is the first discussion of the localization of a core excitation of O$_3$. This result explains the success of the widely used assumption of localized core excitation in adsorbates and large molecules.

An estimate of the inter-system crossing time in light-emitting polymers

Abstract: The reported enhanced singlet-exciton yields in light-emitting polymers over the statistical limit of 25% has attracted wide experimental and theoretical attention. Most theoretical estimates of the singlet-exciton yield depend crucially on estimates of the inter-system crossing rates induced by spin-orbit coupling. In this paper we use the experimentally determined phosphorescent life-time and energy of the lowest-lying triplet state, as well as calculated values of Huang-Rhys factors to estimate the spin-orbit matrix element and inter-system crossing time between the lowest-lying singlet and triplet states.

Strongly localised molecular orbitals for $\alpha$-quartz

Abstract: A previously proposed computational procedure for constructing a set of nonorthogonal strongly localised one-electron molecular orbitals (O. Danyliv, L. Kantorovich - physics/0401107) is applied to a perfect $\alpha$-quartz crystal characterised by an intermediate type of chemical bonding. The orbitals are constructed by applying various localisation methods to canonical Hartree-Fock orbitals calculated for a succession of finite molecular clusters of increased size with appropriate boundary conditions. The calculated orbitals span the same occupied Fock space as the canonical HF solutions, but have an advantage of reflecting the true chemical nature of the bonding in the system. The applicability of several localisation techniques as well as of a number of possible choices of localisation regions (structure elements) are discussed for this system in detail.

Inspecting plastic deformation of Pd by means of fractal geometry

Abstract: The influence of phase transformation-induced plastic deformation in Pd|H system on the electrode surface was investigated. Since the Pd surface is subject of severe plastic deformation during this process, the structure and roughness of the electrode surface significantly changes. Quantitative analysis of the electrode surfaces for comparative study of such changes is a valuable tool to inspect the plastic deformation induced. Fractal dimension can be used as a quantitative measure for this purpose. Since inappropriate methods may lead to significant errors, an appropriate approach was proposed for the determination of fractal dimensions in such systems. It was demonstrated that the surface roughness generated is mainly due to the plastic deformation induced, not the other side processes. The methodology is of general interest for the investigation of plastic deformations. The methodology is of general interest for the investigation of plastic deformations, and also other solid-state structural changes.

Neural Network Approach to Study of Langmuir Layers at the Water Surface

Abstract: The structure of the layers of a phthalocyanines is studied through a measurement of surface pressure-area isotherms at various initial surface concentrations of the spreaded compound. Layer characteristics are determined from experimental data by replotting the isotherms in accordance with M. Volmer s equation. Neural networks taught on approximations of the isoterms. The method is compared with the traditional least square one. Another task for NN is a forecast of main characteristics of the layers.

Theoretical calculation of the surface energy of water

Abstract: The estimation of the surface tension of water is theoretically dealt on the basis of the dipole molecular model. It is known that the experimentally determined surface tension of freshly exposed surface has a higher value than the nominal value of 73 [mN/m]. We calculated the value corresponding to the fresh surface where the reorientation of the molecules has not occurred.