Statistics of copolymer sequences obtained by polymeranalogous reactions of homopolymer globules

For the first time the macrokinetics of polymeranalogous reaction in a homopolymer globule is considered theoretically. Using mathematical apparatus of the theory of markovian stochastic processes we examined the statistical characteristics of the chemical structure of the emerging copolymers. It is shown that when the reaction is kinetically-controlled the sequence distribution of monomeric units in copolymer chains follows Bernoullian statistics. In the opposite regime, when the reaction is diffusion-controlled it is proved that the resulting primary structure is described by the Levy flight distribution. Such copolymers with nonmarkovian character of the alteration of units in macromolecules can not be synthesized in the bulk or in concentrated solutions by means of standard synthetic polymerization procedures. Received 14 June 2002 Published online: 16 April 2003 RID="a" ID="a"e-mail: khokhlov@polly.phys.msu.ru     

The lifetime of the soliton in the improved Davydov model at the biological temperature 300 K for protein molecules

We study the effects of quantum fluctuations and thermal perturbations on the lifetime of the soliton in the improved Davydov model proposed by us with two-quanta and with an added interaction. By using quantum perturbation theory, we compute the soliton lifetime for a wide ranges of parameter values relevant for protein molecules. The lifetime of the new soliton at the biological temperature 300 K is of the order of 10^-10 second or τ/τ≥ 500 for parameters appropriate to α-helical protein molecules. This shows clearly that the new soliton in the improved model is a viable mechanism for the bio-energy transport in the α-helix region of proteins. Received 7 January 1999 and Received in final form 16 August 2000     

Drift of a polymer chain in a porous medium --A Monte Carlo study

We investigate the drift of an end-labeled telehelic polymer chain in a frozen disordered medium under the action of a constant force applied to the one end of the macromolecule by means of an off-lattice bead spring Monte Carlo model. The length of the polymers N is varied in the range 8 < N < 128, and the obstacle concentration in the medium C is varied from zero up to the percolation threshold C≈ 0.75. For field intensities below a C-dependent critical field strength B _c, where jamming effects become dominant, we find that the conformational properties of the drifting chains can be interpreted as described by a scaling theory based on Pincus blobs. The variation of drag velocity with C in this interval of field intensities is qualitatively described by the law of Mackie-Meares. The threshold field intensity B _c itself is found to decrease linearly with C. Received 20 August 2001 and Received in final form 19 November 2001     

Helix coil mixing in twist-storing polymers

Twist-storing polymers respond with elastic energy penalty to coherent or random twisting along the local chain axis away from its equilibrium, which can be straight (as in “ribbons”) or helical (as in DNA and other biopolymers). Here we study the equilibrium conformation of such polymers, focusing on the thermodynamic balance between twist and writhe, resulting from the competition between the random coil entropy and the potential energy stored in superhelical portions of the polymer chain. Two macroscopic variables characterise such a chain, the end-to-end distance R and the link number Lk, which is a topological invariant of a given polymer with clamped ends. We find that with increasing link number Lk, the chain accommodates its excess twist in growing plectonemes, unless forced out of this state by stretching its end-to-end distance R. We calculate the force-extension relation, which exhibits crossovers between different deformation regimes. Received 16 November 2000 and Received in final form 6 February 2001     

The electrostatic persistence length of polymers beyond the OSF limit

We use large-scale Monte Carlo simulations to test scaling theories for the electrostatic persistence length l _e of isolated, uniformly charged polymers with Debye-Hückel intrachain interactions in the limit where the screening length κ^-1 exceeds the intrinsic persistence length of the chains. Our simulations cover a significantly larger part of the parameter space than previous studies. We observe no significant deviations from the prediction l _e∝κ^-2 by Khokhlov and Khachaturian which is based on applying the Odijk-Skolnick-Fixman theories of electrostatic bending rigidity and electrostatically excluded volume to the stretched de Gennes-Pincus-Velasco-Brochard polyelectrolyte blob chain. A linear or sublinear dependence of the persistence length on the screening length can be ruled out. We show that previous results pointing into this direction are due to a combination of excluded-volume and finite chain length effects. The paper emphasizes the role of scaling arguments in the development of useful representations for experimental and simulation data. Received 12 February 2002     

Elasticity theory of a twisted stack of plates

We present an elastic model of B-form DNA as a stack of thin, rigid plates or base pairs that are not permitted to deform. The symmetry of DNA and the constraint of plate rigidity limit the number of bulk elastic constants contributing to a macroscopic elasticity theory of DNA to four. We derive an effective twist-stretch energy in terms of the macroscopic stretch along and relative excess twist about the DNA molecular axis. In addition to the bulk stretch and twist moduli found previously, we obtain a twist-stretch modulus with the following remarkable properties: 1) it vanishes when the radius of the helical curve following the geometric center of each plate is zero, 2) it vanishes with the elastic constant K₂₃ that couples compression normal to the plates to a shear strain, if the plates are perpendicular to the molecular axis, and 3) it is nonzero if the plates are tilted relative to the molecular axis. This implies that a laminated helical structure carved out of an isotropic elastic medium will not twist in response to a stretching force, but an isotropic material will twist if it is bent into the shape of a helix. Received: 4 July 1997 / Received in final form: 16 October 1997 / Accepted: 21 October 1997     

Rayleigh instability of charged droplets in the presence of counterions

The classical instability of a charged spherical droplet is reconsidered in the presence of counterions. An ensemble of such droplets is studied within a simplified cell model. Screening of the electric field by the counterions is found to increase the equilibrium droplet size. Furthermore, if the ions can enter the droplet, a first-order phase transition occurs upon increasing Bjerrum length, surface tension or droplet density, leading to a phase separation. Simple scaling properties of the free energy give the shape of the phase boundary and show the system to be scale-invariant there. Pearl-necklace structures of hydrophobic polyelectrolytes are discussed as an application. Received 30 August 2001     

Adsorption of a semiflexible wormlike chain

Adsorption of ideal polymers with stiff backbone onto a flat surface is considered theoretically. Both scaling approach and quantitative theory are developed. We predict a self-similar monomer concentration profile c(x) ∼ x ^-4/3 near the surface (when the distance to the surface x is much smaller than the chain persistence length l /2). The typical conformation of a weakly adsorbed chain can be viewed as a sequence of alternating flat (2-dimensional) trains of wormlike short loops (flat blobs) and coil-like (3-dimensional) loops forming a triple-layer structure: contact layer (x < Δ) of adsorbed fragments virtually laid on the surface, proximal layer (Δ < x < l) of flat blobs, and more dilute distal corona layer (x > l). Here Δ defines the range of monomer/surface attraction, Δ ≪ l. The adsorption transition is continuous. However, its relative width is small (T ^* is the adsorption temperature, ΔT is the relevant temperature interval): ∼ , i.e. a discontinuous transition in the limit Δ/l↦ 0. Received 10 October 2002 and Received in final form 22 November 2002 RID="a" ID="a"Permanent address: Physics Department, Moscow State University, Moscow 119992, Russia. e-mail: semenov@polly.phys.msu.ru     

Effects of properties of the surfactant on its aggregate behavior

The effects of the solvent-surfactant interaction, chain length and stiffness of surfactants on the formation of aggregates and the aggregation degree of surfactants in the two-dimension solution have been investigated using discontinuous molecular-dynamics simulations. When the tail-water repulsion increases or the head-water attraction decreases, the aggregation degree increases. Increasing the chain length and raising the stiffness of surfactants lead to the increment of the aggregation degree. Received 25 July 2001 and Received in final form 7 May 2002     

Measuring the implantation depth of silver clusters in graphite

Scanning tunnelling microscopy (STM) and molecular dynamics (MD) simulations have been used to investigate the implantation of Ag₇ ^- clusters into the graphite surface. An experimental measure of the implantation depth of individual clusters is gained via thermal oxidation of the bombarded graphite surfaces. This process results in etching of the cluster-induced defects to form etch pits which grow laterally whilst retaining the depth of the implanted cluster. STM imaging of the etch pits reveals the distribution of implantation depths for deposition energies of 2 keV and 5 keV. Molecular dynamics simulations for clusters of 5 keV energy show that the implantation depth for Ag₇ ^- is largely independent of the impact site on the graphite surface and the cluster orientation. The implantation depth found by MD lies at the upper edge of the experimental depth distribution. Received 30 November 2000     

Weak violation of universality for polyelectrolyte chains: Variational theory and simulations

A variational approach is considered to calculate the free energy and the conformational properties of a polyelectrolyte chain in d dimensions. We consider in detail the case of pure Coulombic interactions between the monomers, when screening is not present, in order to compute the end-to-end distance and the asymptotic properties of the chain as a function of the polymer chain length N. We find R≃N ^ν(log N)^γ, where ν = and λ is the exponent which characterizes the long-range interaction U∝ 1/r ^λ. The exponent γ is shown to be non-universal, depending on the strength of the Coulomb interaction. We check our findings by a direct numerical minimization of the variational energy for chains of increasing size 2⁴ < N < 2¹⁵. The electrostatic blob picture, expected for small enough values of the interaction strength, is quantitatively described by the variational approach. We perform a Monte Carlo simulation for chains of length 2⁴ < N < 2¹⁰. The non-universal behavior of the exponent γ previously derived within the variational method is also confirmed by the simulation results. Non-universal behavior is found for a polyelectrolyte chain in d = 3 dimension. Particular attention is devoted to the homopolymer chain problem, when short-range contact interactions are present. Received 8 August 2000 and Received in final form 19 December 2000     

Force-extension behavior of folding polymers

The elastic response of flexible polymers made of elements which can be either folded or unfolded, having different lengths in these two states, is discussed. These situations are common for biopolymers as a result of folding interactions intrinsic to the monomers, or as a result of binding of other smaller molecules along the polymer length. Using simple flexible-chain models, we show that even when the energy ε associated with maintaining the folded state is comparable to k _B T, the elastic response of such a chain can mimic usual polymer linear elasticity, but with a force scale enhanced above that expected from the flexibility of the chain backbone. We discuss recent experiments on single-stranded DNA, chromatin fiber and double-stranded DNA with proteins weakly absorbed along its length which show this effect. Effects of polymer semiflexiblity and torsional stiffness relevant to experiments on proteins binding to dsDNA are analyzed. We finally discuss the competition between electrostatic self-repulsion and folding interactions responsible for the complex elastic response of single-stranded DNA. Received 7 August 2002 and Received in final form 7 March 2003 / Published online: 15 April 2003 RID="a" ID="a"e-mail: jmarko@uic.edu     

Optical investigation of the quasi two-dimensional monophosphate tungsten bronzes K x P 4 W 8 O 32 (x = 0 - 1.57)

The potassium doped monophosphate tungsten bronzes K_xP₄W₈O₃₂ are two-dimensional metals which show a metal-to-metal transition at a critical temperature which depends on the doping level. The metal-to-metal transition is accompanied by the formation of a commensurate charge density wave with wave vector (π/b,0) which is independent of the doping level. Undoped P₄W₈O₃₂, on the other hand, has two metal-to-metal transitions which are connected to the formation of incommensurate charge density waves. We measured the infrared reflectivity of the series K_xP₄W₈O₃₂ (x = 0 - 1.57) in the spectral range from 100 to 10 000 cm^-1 for room temperature and well below the critical temperature. Polarization-dependent infrared spectra find a two-dimensional behavior in the normal and the charge density wave state and show signatures of hybridization between one- and two-dimensional conduction bands. In undoped P₄W₈O₃₂ the essentials of the charge density wave state can be understood from the nesting vectors of the calculated Fermi surface and two gaps are observed in the infrared spectra. The gap sizes are a factor of about 2.5 bigger than the predictions from mean-field theory in the weak-coupling limit which suggests medium- or strong electron-phonon coupling. For potassium doped K_xP₄W₈O₃₂ one gap is observed in the charge density wave state. The energetics of the charge density formation may be dominated by the energy required for the lattice modulation. Received 27 April 2001 and Received in final form 21 September 2001     

Influence of the helium autoionization structures on the single/double ionization signal

Calculations of intense field (around 10¹⁶ W/cm²) single- and double-ionization processes in helium at XUV wavelengths are presented. The laser wavelength is chosen near the | 2s2p ¹ P autoionization structure and the dynamics are explored. Single and double ionization yields, as well as the photoelectron energy spectrum for photon energies around the autoionization structure are calculated. In the case of a pulse of few femtoseconds duration, no significant enhancement of the double ionization yield has been found in tuning the photon frequency around the peak of the resonance. It is also shown that in the case of a long pulse (and hence narrow compared with the relevant autoionization width), the branching ratio of double to single ionization yield can be relatively enhanced by tuning to the absorption minimum of the resonance. Received 19 February 2002 / Received in final form 2 May 2002 Published online 19 July 2002     

Delocalization effects and charge reorganizations induced by repulsive interactions in strongly disordered chains

We study the delocalization effect of a short-range repulsive interaction on the ground state of a finite density of spinless fermions in strongly disordered one dimensional lattices. The density matrix renormalization group method is used to explore the charge density and the sensitivity of the ground state energy with respect to the boundary condition (the persistent current) for a wide range of parameters (carrier density, interaction and disorder). Analytical approaches are developed and allow to understand some mechanisms and limiting conditions. For weak interaction strength, one has a Fermi glass of Anderson localized states, while in the opposite limit of strong interaction, one has a correlated array of charges (Mott insulator). In the two cases, the system is strongly insulating and the ground state energy is essentially invariant under a twist of the boundary conditions. Reducing the interaction strength from large to intermediate values, the quantum melting of the solid array gives rise to a more homogeneous distribution of charges, and the ground state energy changes when the boundary conditions are twisted. In individual chains, this melting occurs by abrupt steps located at sample-dependent values of the interaction where an (avoided) level crossing between the ground state and the first excitation can be observed. Important charge reorganizations take place at the avoided crossings and the persistent currents are strongly enhanced around the corresponding interaction value. These large delocalization effects become smeared and reduced after ensemble averaging. They mainly characterize half filling and strong disorder, but they persist away of this optimal condition. Received 5 July 2000 and Received in final form 8 November 2000     

Experimental and theoretical analysis of the 5 pnpJ =0 e , 1 e , 2 e autoionizing spectrum of Sr

The even parity 5 pnpJ = 0, 1 and 2 doubly excited autoionizing states of strontium were investigated both experimentally and theoretically. Sr atoms in an atomic beam were excited through the two-step Isolated Core Excitation (ICE) scheme 5s^{2 1} S ₀ -λ ₁ → 5 sn'p ¹ P ₁ ( n' = 12-16)-λ ₂ → [5 p _3/2 np ] _J . The final ICE transition probes the [5 p _3/2 np ] _J resonances. However, the [5 p _1/2 np ] _J series below the 5 p _1/2 threshold were excited also due to their mixing with the [5 p _3/2 np ] _J perturbers. An extended energy region was covered below and above the 5 p _1/2 ionization limit by saturating the central ICE lobe and recording as many as possible “red” and “blue” secondary lobes. J identification was achieved by using mutually parallel and perpendicular linear polarizations of the laser beams. The ICE spectra were compared to those obtained by employing a two-step excitation scheme using the bound 4 d 5 p ¹ P ₁ valence state as an intermediate one. Final identification for very complex structures was achieved after comparison with theoretical energy level positions and excitation profiles produced by the R-matrix method combined with the multichannel quantum defect theory (MQDT) method. The agreement between theoretical and observed structures is quite satisfactory. Received 31 May 2000     

The multiphoton ionization rate and the energy shift of atoms interacting with weak dichromatic fields with commensurate frequencies are simple functions of the phase difference

By implementing a time-independent, nonperturbative many-electron, many-photon theory (MEMPT), cycle-averaged complex eigenvalues were obtained for the He atom, whose real part gives the field-induced energy shift, Δ(ω ₁, F ₁;ω ₂, F ₂,ϕ), and the imaginary part is the multiphoton ionization rate, Γ(ω ₁, F ₁;ω ₂, F ₂,ϕ), where ω is the frequency, F is the field strength and ϕ is the phase difference. Through analysis and computation we show that, provided the intensities are weak, the dependence of Γ(ω ₁, F ₁;ω ₂, F ₂,ϕ) on ϕ is simple. Specifically, for odd harmonics, Γ varies linearly with cos(ϕ) whilst for even harmonics it varies linearly with cos(2ϕ). In addition, this dependence on ϕ holds for Δ(ω ₁, F ₁;ω ₂, F ₂,ϕ) as well. These relations may turn out to be applicable to other atomic systems as well, and to provide a definition of the weak field regime in the dichromatic case. When the combination of (ω ₁, F ₁) and (ω ₂, F ₂) is such that higher powers of cos(ϕ) and cos(2ϕ) become important, these rules break down and we reach the strong field regime. The herein reported results refer to Γ(ω ₁, F ₁;ω ₂, F ₂,ϕ) and Δ(ω ₁, F ₁;ω ₂, F ₂,ϕ) for He irradiated by a dichromatic ac-field consisting of the fundamental wavelength λ = 248 nm and its 2nd, 3rd and 4th higher harmonics. The intensities are in the range 1.0×10¹²-3.5×10¹⁴ W/cm², with the intensity of the harmonics being 1-2 orders of magnitude smaller. The calculations incorporated systematically electronic structure and electron correlation effects in the discrete and in the continuous spectrum, for ¹S, ¹P, ¹D, ¹F, ¹G, and ¹H two-electron states of even and odd parity. Received 9 July 2000 and Received in final form 2 November 2000     

Variational calculations of ionized-donor-bound excitons in GaAs-Al x Ga 1 - x As quantum wells

Using a two-parameter wave function, we calculate variationally the binding energy of an exciton bound to an ionized donor impurity (D^+,X) in GaAs-Al_xGa_1-xAs quantum wells for the values of the well width from 10 to 300 ?, when the dopant is located in the center of the well and at the edge of the well. The theoretical results confirm that the previous experimental speculation proposed by Reynolds et al. [Phys. Rev. B 40, 6210 (1989)] is the binding energy of D^+,X for the dopant at the edge of the well. In addition, we also calculate the center-of-mass wave function of the exciton and the average interparticle distances. The results are discussed in detail. Received 17 July 2000 and Received in final form 13 November 2000     

The dynamics of a THz Rydberg wavepacket

An optically excited Rydberg wavepacket can be generated by exciting the electron from a low-lying state to a coherent superposition of high-lying states with a short broadband optical pulse. A special kind of Rydberg wavepacket is generated in the case of a interaction of a weak THz half cycle pulse with a stationary Rydberg state, called the THz wavepacket. This THz wavepacket is a coherent superposition of the initial Rydberg state and its neighbouring states. We have investigated the time evolution of THz wavepackets by measuring the impact of two in time delayed half cycle pulses ( ≈ 200 V cm^-1) on the population of a stationary (n = 40) Rydberg state in rubidium. The first half cycle pulse creates the THz wavepacket and the second half cycle pulse probes the dynamics of the THz wavepacket. We support our experimental data by numerically solving the Schr?dinger equation and with a semi-classical picture. Whereas an optically excited wavepacket is initially localized, a THz wavepacket is initially delocalized and becomes localized after half a revival time. Received 23 August 2000 and Received in final form 27 March 2001     

Coherent transport through intramolecular junction of single-wall carbon nanotubes

We have constructed four types single-wall carbon nanotube intramolecular junctions (IMJs) of (5,5)/(8,0), (5,5)/(10,0), (5,5)/(9,0)A, and (5,5)/(9,0)B along a common axis, and calculated their electronic and transport properties using a tight binding-based Green's function approach that is particular suitable for realistic calculation of electronic transport property in extended system. Our results show that quasi-localized states can appear in the metal/semiconductor heterojunctions ((5,5)/(8,0) and (5,5)/(10,0)junctions), which is desirable for the design of a quantum device; and the conductance of M-M IMJs is very sensitive to the connectivity of the matching tubes, certain configurations of connection completely stop the flow of electron, while others permit the transmission of the current through the interface. These results may have implications for the device assembly and manipulation process of all carbon nanotubes-based microelectronic elements. Received 14 January 2003 / Received in final form 25 February 2003 Published online 4 June 2003 RID="a" ID="a"e-mail: lfyzz@yahoo.com.cn