BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Tapio Ala-Nissila

Tapio Ala-Nissila appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.softcond-mat.stat-mechcond-mat.mes-hallcond-mat.mtrl-sciphysics.comp-phquant-phBiological Physicsphysics.chem-phphysics.flu-dynBiomoleculesphysics.app-phphysics.opticsphysics.soc-phPopulations and Evolution

Catalog footprint

What is connected

50works
14topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

50 published item(s)

preprint2026arXiv

Density-dependent sodium-storage mechanisms in hard carbon materials

Understanding the sodium-storage mechanism in hard carbon (HC) anodes is crucial for advancing sodium-ion battery (SIB) technology. However, the intrinsic complexity of HC microstructures and their interactions with sodium remain not fully elucidated. We present a multiscale methodology that integrates grand-canonical Monte Carlo (GCMC) simulations with a machine-learning interatomic potential based on the Gaussian approximation potential (GAP) framework to investigate sodium insertion mechanisms in hard carbons with different levels of porosity, achieved by simulating structural models with densities ranging from 0.7 to 1.9 g cm$^{-3}$. Structural and thermodynamic analyses reveal the interplay between pore size and accessibility and the relative contributions of adsorption, intercalation, and pore filling to the overall storage capacity. Low-density carbons favor pore-filling, achieving extremely high capacities at near-zero voltages, whereas high-density carbons primarily store sodium through adsorption and intercalation, leading to lower but more stable capacities. Intermediate-density carbons ($1.3-1.6$ g cm$^{-3}$) provide the most balanced performance, combining moderate capacity (480 and 310 mAh g$^{-1}$), safe operating voltages, and minimal volume expansion ($<10$\%). These findings establish a direct correlation between carbon density and electrochemical behavior, providing atomic-scale insight into how hard carbon morphology governs sodium-storage. The proposed framework offers a rational design principle for optimizing HC-based SIB anodes toward high energy density and long-term cycling stability.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Driven polymer translocation into a channel: Iso-flux tension propagation theory and Langevin dynamics simulations

Iso-flux tension propagation (IFTP) theory and Langevin dynamics (LD) simulations are employed to study the dynamics of channel-driven polymer translocation in which a polymer translocates into a narrow channel and the monomers in the channel experience a driving force $f_{\rm c}$. In the high driving force limit, regardless of the channel width, IFTP theory predicts $τ\propto f_{\textrm{c}}^β$ for the translocation time, where $β=-1$ is the force scaling exponent. Moreover, LD data show that for a very narrow channel fitting only a single file of monomers, the entropic force due to the subchain inside the channel does not play a significant role in the translocation dynamics, and the force exponent $β= -1$ regardless of the force magnitude. As the channel width increases the number of possible spatial configurations of the subchain inside the channel becomes significant, and the resulting entropic force causes the force exponent to drop below unity.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Modified Poisson-Boltzmann theory for polyelectrolytes in monovalent salt solutions with finite-size ions

We present a soft-potential-enhanced Poisson-Boltzmann (SPB) theory to efficiently capture ion distributions and electrostatic potential around rodlike charged macromolecules. The SPB model is calibrated with a coarse-grained particle-based model for polyelectrolytes (PEs) in monovalent salt solutions as well as compared to a full atomistic molecular dynamics simulations with explicit solvent. We demonstrate that our modification enables the SPB theory to accurately predict monovalent ion distributions around a rodlike PE in a wide range of ion and charge distribution conditions in the weak-coupling regime. These include excess salt concentrations up to 1 M, and ion sizes ranging from small ions, such as Na or Cl, to softer and larger ions with size comparable to the PE diameter. The work provides a simple way to implement an enhancement that effectively captures the influence of ion size and species into the PB theory in the context of PEs in aqueous salt solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Neuroevolution machine learning potentials: Combining high accuracy and low cost in atomistic simulations and application to heat transport

We develop a neuroevolution-potential (NEP) framework for generating neural network based machine-learning potentials. They are trained using an evolutionary strategy for performing large-scale molecular dynamics (MD) simulations. A descriptor of the atomic environment is constructed based on Chebyshev and Legendre polynomials. The method is implemented in graphic processing units within the open-source GPUMD package, which can attain a computational speed over $10^7$ atom-step per second using one Nvidia Tesla V100. Furthermore, per-atom heat current is available in NEP, which paves the way for efficient and accurate MD simulations of heat transport in materials with strong phonon anharmonicity or spatial disorder, which usually cannot be accurately treated either with traditional empirical potentials or with perturbative methods.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Temperature-resilient anapole modes associated with TE polarization in semiconductor nanowire

Polarization-dependent scattering anisotropy of cylindrical nanowires has numerous potential applications in, for example, nanoantennas, photothermal therapy, thermophotovoltaics, catalysis, sensing, optical filters and switches. In all these applications, temperature-dependent material properties play an important role and often adversely impact performance depending on the dominance of either radiative or dissipative damping. Here, we employ numerical modeling based on Mie scattering theory to investigate and compare the temperature and polarization-dependent optical anisotropy of metallic (gold, Au) nanowires with indirect (silicon, Si) and direct (gallium arsenide, GaAs) bandgap semiconducting nanowires. Results indicate that plasmonic scattering resonances in semiconductors, within the absorption band, deteriorate with an increase in temperature whereas those occurring away from the absorption band strengthen as a result of the increase in phononic contribution. Indirect-bandgap thin ($20 \,\mathrm{nm}$) Si nanowires present low absorption efficiencies for both the transverse electric (TE, $E_{\perp}$) and magnetic (TM, $E_{\parallel}$) modes, and high scattering efficiencies for the TM mode at shorter wavelengths making them suitable as highly efficient scatterers. Temperature-resilient higher-order anapole modes with their characteristic high absorption and low scattering efficiencies are also observed in the semiconductor nanowires ($r \! = \! 125 \! - \! 130$ nm) for the TE polarization. Herein, the GaAs nanowires present $3 \! - \! 7$ times greater absorption efficiencies compared to the Si nanowires making them especially suitable for temperature-resilient applications such as scanning near-field optical microscopy (SNOM), localized heating, non-invasive sensing or detection that require strong localization of energy in the near field.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Adaptive and optimized COVID-19 vaccination strategies across geographical regions and age groups

We evaluate the efficiency of various heuristic strategies for allocating vaccines against COVID-19 and compare them to strategies found using optimal control theory. Our approach is based on a mathematical model which tracks the spread of disease among different age groups and across different geographical regions, and we introduce a method to combine age-specific contact data to geographical movement data. As a case study, we model the epidemic in the population of mainland Finland utilizing mobility data from a major telecom operator. Our approach allows to determine which geographical regions and age groups should be targeted first in order to minimize the number of deaths. In the scenarios that we test, we find that distributing vaccines demographically and in an age-descending order is not optimal for minimizing deaths and the burden of disease. Instead, more lives could potentially be saved by using strategies which emphasize high-incidence regions and distribute vaccines in parallel to multiple age groups. The level of emphasis that high-incidence regions should be given depends on the overall transmission rate in the population. This observation highlights the importance of updating the vaccination strategy when the effective reproduction number changes due to the general contact patterns changing and new virus variants entering.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Heat transport across graphene/hexagonal-BN tilted grain boundaries from phase-field crystal model and molecular dynamics simulations

We study the interfacial thermal conductance of grain boundaries (GBs) between monolayer graphene and hexagonal boron nitride (h-BN) sheets using a combined atomistic approach. First, realistic samples containing graphene/h-BN GBs with different tilt angles are generated using the phase-field crystal (PFC) model developed recently [P. Hirvonen \textit{et al.}, Phys. Rev. B \textbf{100}, 165412 (2019)] that captures slow diffusive relaxation inaccessible to molecular dynamics (MD) simulations. Then, large-scale MD simulations using the efficient GPUMD package are performed to assess heat transport and rectification properties across the GBs. We find that lattice mismatch between the graphene and h-BN sheets plays a less important role in determining the interfacial thermal conductance as compared to the tilt angle. In addition, we find no significant thermal rectification effects for these GBs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Directing Near-Infrared Photon Transport with Core@Shell Particles

Directing the propagation of near-infrared radiation is a major concern in improving the efficiency of solar cells and thermal insulators. A facile approach to scatter light in the near-infrared region without excessive heating is to embed compact layers with semiconductor particles. The directional scattering by semiconductor@oxide (core@shell) spherical particles (containing Si, InP, TiO$_2$, SiO$_2$, or ZrO$_2$) with a total radius varying from 0.1 to 4.0 μm and in an insulating medium at low volume fraction is investigated using Lorenz-Mie theory and multiscale modelling. The optical response of each layers is calculated under irradiation by the sun or a blackbody emitter at 1180 K. Reflectance efficiency factors of up to 83.7% and 63.9% are achieved for near-infrared solar and blackbody radiation in 200 μm thick compact layers with only 1% volume fraction of bare Si particles with a radius of 0.23 μm and 0.50 μm, respectively. The maximum solar and blackbody efficiency factors of layers containing InP particles was slightly less (80.2% and 60.7% for bare particles with a radius of 0.25 μm and 0.60 μm, respectively). The addition of an oxide coating modifies the surrounding dielectric environment, which improves the solar reflectance efficiency factor to over 90% provided it matches the scattering mode energies with the incident spectral density. The layers are spectrally-sensitive and can be applied as a back or front reflector for solar devices, high temperature thermal insulators, and optical filters in Gradient Heat Flux Sensors for fire safety applications.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Interpretation of apparent thermal conductivity in finite systems from equilibrium molecular dynamics simulations

We propose a way to properly interpret the apparent thermal conductivity obtained for finite systems using equilibrium molecular dynamics simulations (EMD) with fixed or open boundary conditions in the transport direction. In such systems the heat current autocorrelation function develops negative values after a correlation time which is proportional to the length of the simulation cell in the transport direction. Accordingly, the running thermal conductivity develops a maximum value at the same correlation time and eventually decays to zero. By comparing EMD with nonequilibrium molecular dynamics (NEMD) simulations, we conclude that the maximum thermal conductivity from EMD in a system with domain length 2L is equal to the thermal conductivity from NEMD in a system with domain length L. This facilitates the use of nonperiodic-boundary EMD for thermal transport in finite samples in close correspondence to NEMD.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Kinetic roughening of the urban skyline

\begin{abstract} In this Letter we follow the asymptotic spatial correlation of buildings' heights $G_{\infty}(r)$ in the whole of the Netherlands \cite{bag3d}, which comprises $\approx 10,000,000$ buildings, for the purpose of recovering its scaling with respect to space and time given respectively by the exponents $r^{2 α}$ and $t^{2β}$. This allows us to identify the universality class of the evolution of the urban skyline seen as a dynamically evolving interface. Two major classes of cities were identified based on the recovered value of $α=0.4$ and $α= 0$, which correspond respectively to the KPZ and the EW universality classes. Picking a discrete model from each of these classes and mapping it to physical rules for constructions in cities we conclude that imposed restrictions on buildings' heights are reflected in the exponent $α$ and thus have implications on how the skyline evolves.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Polymer translocation through nanopore assisted by an environment of active rods

We use a combination of computer simulations and iso-flux tension propagation (IFTP) theory to investigate translocation dynamics of a flexible linear polymer through a nanopore into an environment composed of repulsive active rods in 2D. We demonstrate that the rod activity induces a crowding effect on the polymer, leading to a time-dependent net force that facilitates translocation into the active environment. Incorporating this force into the IFTP theory for pore-driven translocation allows us to characterise translocation dynamics in detail and derive a scaling form for the average translocation time as $\tildeτ \sim \tilde{L}_{\textrm{r}}^ν / \tilde{F}_{\textrm{SP}} $, where $\tilde{L}_{\textrm{r}}$ and $\tilde{F}_{\textrm{SP}}$ are the rod length and self-propelling force acting on the rods, respectively, and $ν$ is the Flory exponent.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Pulling a DNA molecule through a nanopore embedded in an anionic membrane: tension propagation coupled to electrostatics

We consider the influence of electrostatic forces on driven translocation dynamics of a flexible polyelectrolyte being pulled through a nanopore by an external force on the head monomer. To this end, we augment the iso-flux tension propagation (IFTP) theory with electrostatics for a negatively charged biopolymer pulled through a nanopore embedded in a similarly charged anionic membrane. We show that for the realistic case such as a single-stranded DNA, the translocation dynamics at low salt where screening is weak and at finite negative membrane charge is unexpectedly accelerated despite the large repulsive electrostatic interactions between the polymer coil on the {\it cis} side and the charged membrane. This is due to the rapid release of the electrostatic potential energy of the coil during translocation, leading to an effectively attractive force that assists end-driven translocation. The speedup results in non-monotonic polymer length and membrane charge dependence of the exponent $α$ characterizing the translocation time $τ\propto N_0^α$ of the polymer with length $N_0$. In the regime of long polymers $N_0\gtrsim500$, the translocation exponent exceeds its upper limit $α=2$ previously observed for the same system without electrostatic interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

State leakage during fast decay and control of a superconducting transmon qubit

Superconducting Josephson junction qubits constitute the main current technology for many applications, including scalable quantum computers and thermal devices. Theoretical modeling of such systems is usually done within the two-level approximation. However, accurate theoretical modeling requires taking into account the influence of the higher excited states without limiting the system to the two-level qubit subspace. Here, we study the dynamics and control of a superconducting transmon using the numerically exact stochastic Liouville-von Neumann equation approach. We focus on the role of state leakage from the ideal two-level subspace for bath induced decay and single-qubit gate operations. We find significant short-time state leakage due to the strong coupling to the bath. We quantify the leakage errors in single-qubit gates and demonstrate their suppression with DRAG control for a five-level transmon in the presence of decoherence. Our results predict the limits of accuracy of the two-level approximation and possible intrinsic constraints in qubit dynamics and control for an experimentally relevant parameter set.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Thermal conductivity reduction in carbon nanotube by fullerene encapsulation: A molecular dynamics study

Single-walled carbon nanotubes (SWCNTs) in their pristine form have high thermal conductivity whose further improvement has attracted a lot of interest. Some theoretical studies have suggested that the thermal conductivity of a $(10,10)$ SWCNT is dramatically enhanced by C$_{60}$ fullerene encapsulation. However, recent experiments on SWCNT bundles show that fullerene encapsulation leads to a reduction rather than an increase in thermal conductivity. Here, we employ three different molecular dynamics methods to study the influence of C$_{60}$ encapsulation on heat transport in a $(10,10)$ SWCNT. All the three methods consistently predict a reduction of the thermal conductivity of $(10,10)$ SWCNT upon C$_{60}$ encapsulation by $20\%-30\%$, in agreement with experimental results on bundles of SWCNTs. We demonstrate that there is a simulation artifact in the Green-Kubo method which gives anomalously large thermal conductivity from artificial convection. Our results show that the C$_{60}$ molecules conduct little heat compared to the outer SWCNT and reduce the phonon mean free paths of the SWCNT by inducing extra phonon scattering. We also find that the thermal conductivity of a $(10,10)$ SWCNT monotonically decreases with increasing filling ratio of C$_{60}$ molecules.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A minimal Tersoff potential for diamond silicon with improved descriptions of elastic and phonon transport properties

Silicon is an important material and many empirical interatomic potentials have been developed for atomistic simulations of it. Among them, the Tersoff potential and its variants are the most popular ones. However, all the existing Tersoff-like potentials fail to reproduce the experimentally measured thermal conductivity of diamond silicon. Here we propose a modified Tersoff potential and develop an efficient open source code called GPUGA (graphics processing units genetic algorithm) based on the genetic algorithm and use it to fit the potential parameters against energy, virial and force data from quantum density functional theory calculations. This potential, which is implemented in the efficient open source GPUMD (graphics processing units molecular dynamics) code, gives significantly improved descriptions of the thermal conductivity and phonon dispersion of diamond silicon as compared to previous Tersoff potentials and at the same time well reproduces the elastic constants. Furthermore, we find that quantum effects on the thermal conductivity of diamond silicon at room temperature are non-negligible but small: using classical statistics underestimates the thermal conductivity by about 10\% as compared to using quantum statistics.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Influence of Boundaries and Thermostatting on Nonequilibrium Molecular Dynamics Simulations of Heat Conduction in Solids

Nonequilibrium molecular dynamics (NEMD) has been extensively used to study thermal transport at various length scales in many materials. In this method, two local thermostats at different temperatures are used to generate a nonequilibrium steady state with a constant heat flux. Conventionally, the thermal conductivity of a finite system is calculated as the ratio between the heat flux and the temperature gradient extracted from the linear part of the temperature profile away from the local thermostats. Here we show that, with a proper choice of the thermostat, the nonlinear part of the temperature profile should actually not be excluded in thermal transport calculations. We compare NEMD results against those from the atomistic Green's function method in the ballistic regime, and those from the homogeneous nonequilibrium molecular dynamics method in the ballistic-to-diffusive regime. These comparisons suggest that in all the transport regimes, one should directly calculate the thermal conductance from the temperature difference between the heat source and sink and, if needed, convert it to the thermal conductivity by multiplying it with the system length. Furthermore, we find that the Langevin thermostat outperforms the Nosé-Hoover (chain) thermostat in NEMD simulations because of its stochastic and local nature. We show that this is particularly important for studying asymmetric carbon-based nanostructures, for which the Nosé-Hoover thermostat can produce artifacts leading to unphysical thermal rectification. Our findings are important to obtain correct results from molecular dynamics simulations of nanoscale heat transport as the accuracy of the interatomic potentials is rapidly improving.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Calorimetric measurement of work for a driven harmonic oscillator

A calorimetric measurement has recently been proposed as a promising technique to measure thermodynamic quantities in a dissipative superconducting qubit. These measurements rely on the fact that the system is projected into energy eigenstates whenever energy is exchanged with the environment. This requirement imposes a restriction on the class of systems that can be measured in this way. Here we extend the calorimetric protocol to the measurement of a driven quantum harmonic oscillator. We employ a scheme based on the two-level approximation to define a new work quantity and show how its statistics relates to the standard two-measurement protocol. We find that for the average work the two-level approximation holds in the underdamped regime for short driving times and, in the overdamped regime, for any driving time. However, this approximation fails for the variance and higher moments of work at finite temperatures. Furthermore, we show how to relate the work statistics obtained through this scheme to the work statistics given by the two-measurement protocol.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Driven translocation of a semi-flexible polymer through a nanopore

We study the driven translocation of a semi-flexible polymer through a nanopore by means of a modified version of the iso-flux tension propagation theory (IFTP), and extensive molecular dynamics (MD) simulations. We show that in contrast to fully flexible chains, for semi-flexible polymers with a finite persistence length $\tilde{\ell}_p$ the {\it trans} side friction must be explicitly taken into account to properly describe the translocation process. In addition, the scaling of the end-to-end distance $R_N$ as a function of the chain length $N$ must be known. To this end, we first derive a semi-analytic scaling form for $R_N$, which reproduces the limits of a rod, an ideal chain, and an excluded volume chain in the appropriate limits. We then quantitatively characterize the nature of the {\it trans} side friction based on MD simulations of semi-flexible chains. Augmented with these two factors, the modified IFTP theory shows that there are three main regimes for the scaling of the average translocation time $τ\propto N^α$. In the stiff chain (rod) limit $N/\tilde{\ell}_p \ll 1$, {$α= 2$}, which continuously crosses over in the regime $ 1 < N/\tilde{\ell}_p < 4$ towards the ideal chain behavior with {$α= 3/2$}, which is reached in the regime $N/\tilde{\ell}_p \sim 10^2$. Finally, in the limit $N/\tilde{\ell}_p \gg 10^6$ the translocation exponent approaches its symptotic value $1+ν$, where $ν$ is the Flory exponent. Our results are in good agreement with available simulations and experimental data.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Efficient dynamical correction of the transition state theory rate estimate for a flat energy barrier

The recrossing correction to the transition state theory estimate of a thermal rate can be difficult to calculate when the energy barrier is flat. This problem arises, for example, in polymer escape if the polymer is long enough to stretch between the initial and final state energy wells while the polymer beads undergo diffusive motion back and forth over the barrier. We present an efficient method for evaluating the correction factor by constructing a sequence of hyperplanes starting at the transition state and calculating the probability that the system advances from one hyperplane to another towards the product. This is analogous to what is done in forward flux sampling except that there the hyperplane sequence starts at the initial state. The method is applied to the escape of polymers with up to 64 beads from a potential well. For high temperature, the results are compared with direct Langevin dynamics simulations as well as forward flux sampling and excellent agreement between the three rate estimates is found. The use of a sequence of hyperplanes in the evaluation of the recrossing correction speeds up the calculation by an order of magnitude as compared with the traditional approach. As the temperature is lowered, the direct Langevin dynamics simulations as well as the forward flux simulations become computationally too demanding, while the harmonic transition state theory estimate corrected for recrossings can be calculated without significant increase in the computational effort.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Efficient protocol for qubit initialization with a tunable environment

We propose an efficient qubit initialization protocol based on a dissipative environment that can be dynamically adjusted. Here the qubit is coupled to a thermal bath through a tunable harmonic oscillator. On-demand initialization is achieved by sweeping the oscillator rapidly into resonance with the qubit. This resonant coupling with the engineered environment induces fast relaxation to the ground state of the system, and a consecutive rapid sweep back to off resonance guarantees weak excess dissipation during quantum computations. We solve the corresponding quantum dynamics using a Markovian master equation for the reduced density operator of the qubit-bath system. This allows us to optimize the parameters and the initialization protocol for the qubit. Our analytical calculations show that the ground-state occupation of our system is well protected during the fast sweeps of the environmental coupling and, consequently, we obtain an estimate for the duration of our protocol by solving the transition rates between the low-energy eigenstates with the Jacobian diagonalization method. Our results suggest that the current experimental state of the art for the initialization speed of superconducting qubits at a given fidelity can be considerably improved.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Global transition path search for dislocation formation in Ge on Si(001)

Global optimization of transition paths in complex atomic scale systems is addressed in the context of misfit dislocation formation in a strained Ge film on Si(001). Such paths contain multiple intermediate minima connected by minimum energy paths on the energy surface emerging from the atomic interactions in the system. The challenge is to find which intermediate states to include and to construct a path going through these intermediates in such a way that the overall activation energy for the transition is minimal. In the numerical approach presented here, intermediate minima are constructed by heredity transformations of known minimum energy structures and by identifying local minima in minimum energy paths calculated using a modified version of the nudged elastic band method. Several mechanisms for the formation of a 90° misfit dislocation at the Ge-Si interface are identified when this method is used to construct transition paths connecting a homogeneously strained Ge film and a film containing a misfit dislocation. One of these mechanisms which has not been reported in the literature is detailed. The activation energy for this path is calculated to be 26% smaller than the activation energy for half loop formation of a full, isolated 60° dislocation. An extension of the common neighbor analysis method involving characterization of the geometrical arrangement of second nearest neighbors is used to identify and visualize the dislocations and stacking faults.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Heat flux and information backflow in cold condensed matter systems

We examine non-Markovian effects in an open quantum system from the point of view of information flow. To this end, we consider the spin-boson model with a cold reservoir, accounting for the exact time-dependent correlations between the system and the bath to study the exchange of information and heat. We use an information theoretic measure of the relevant memory effects and demonstrate that the information backflow from the reservoir to the system does not necessarily correlate with the backflow of heat. We also examine the influence of temperature and coupling strength on the loss and gain of information between the system and the bath. Finally, we discuss how additional driving changes the backflow of information, giving rise to potential applications in reservoir engineering.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Minimum energy path for the nucleation of misfit dislocations in Ge/Si(001) heteroepitaxy

A possible mechanism for the formation of a 90° misfit dislocation at the Ge/Si(001) interface through homogeneous nucleation is identified from atomic scale calculations where a minimum energy path connecting the coherent epitaxial state and a final state with a 90° misfit dislocation is found using the nudged elastic band method. The initial path is generated using a repulsive bias activation procedure in a model system including 75000 atoms. The energy along the path exhibits two maxima in the energy. The first maximum occurs as a 60° dislocation nucleates. The intermediate minimum corresponds to an extended 60° dislocation. The subsequent energy maximum occurs as a second 60° dislocation nucleates in a complementary, mirror glide plane, simultaneously starting from the surface and from the first 60° dislocation. The activation energy of the nucleation of the second dislocation is 30% lower than that of the first one showing that the formation of the second 60° dislocation is aided by the presence of the first one. The simulations represent a step towards unraveling the formation mechanism of 90° dislocations, an important issue in the design of growth procedures for strain released Ge overlayers on Si(100) surfaces, and more generally illustrate an approach that can be used to gain insight into the mechanism of complex nucleation paths of extended defects in solids.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Multiscale modeling of polycrystalline graphene: A comparison of structure and defect energies of realistic samples from phase field crystal models

We extend the phase field crystal (PFC) framework to quantitative modeling of polycrystalline graphene. PFC modeling is a powerful multiscale method for finding the ground state configurations of large realistic samples that can be further used to study their mechanical, thermal or electronic properties. By fitting to quantum-mechanical density functional theory (DFT) calculations, we show that the PFC approach is able to predict realistic formation energies and defect structures of grain boundaries. We provide an in-depth comparison of the formation energies between PFC, DFT and molecular dynamics (MD) calculations. The DFT and MD calculations are initialized using atomic configurations extracted from PFC ground states. Finally, we use the PFC approach to explicitly construct large realistic polycrystalline samples and characterize their properties using MD relaxation to demonstrate their quality.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Quantum jump model for a system with a finite-size environment

Measuring the thermodynamic properties of open quantum systems poses a major challenge. A calorimetric detection has been proposed as a feasible experimental scheme to measure work and fluctuation relations in open quantum systems. However, the detection requires a finite size for the environment, which influences the system dynamics. This process cannot be modeled with the standard stochastic approaches. We develop a quantum jump model suitable for systems coupled to a finite-size environment. With the method we study the common fluctuation relations and prove that they are satisfied.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Perfect quantum excitation energy transport via single edge perturbation in a complete network

We consider quantum excitation energy transport (EET) in a network of two-state nodes in the Markovian approximation by employing the Lindblad formulation. We find that EET from an initial site, where the excitation is inserted to the sink, is generally inefficient due to the inhibition of transport by localization of the excitation wave packet in a symmetric, fully-connected network. We demonstrate that the EET efficiency can be significantly increased up to %100 by perturbing hopping transport between the initial node and the one connected directly to the sink, while the rate of energy transport is highest at a finite value of the hopping parameter. We also show that prohibiting hopping between the other nodes which are not directly linked to the sink does not improve the efficiency. We show that external dephasing noise in the network plays a constructive role for EET in the presence of localization in the network, while in the absence of localization it reduces the efficiency of EET.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Polymer escape from a confining potential

The rate of escape of polymers from a two-dimensionally confining potential well has been evaluated using self-avoiding as well as ideal chain representations of varying length, up to 80 beads. Long timescale Langevin trajectories were calculated using the path integral hyperdynamics method to evaluate the escape rate. A minimum is found in the rate for self-avoiding polymers of intermediate length while the escape rate decreases monotonically with polymer length for ideal polymers. The increase in the rate for long, self-avoiding polymers is ascribed to crowding in the potential well which reduces the free energy escape barrier. An effective potential curve obtained using the centroid as an independent variable was evaluated by thermodynamic averaging and Kramers rate theory then applied to estimate the escape rate. While the qualitative features are well reproduced by this approach, it significantly overestimates the rate, especially for the longer polymers. The reason for this is illustrated by constructing a two-dimensional effective energy surface using the radius of gyration as well as the centroid as controlled variables. This shows that the description of a transition state dividing surface using only the centroid fails to confine the system to the region corresponding to the free energy barrier and this problem becomes more pronounced the longer the polymer is. A proper definition of a transition state for polymer escape needs to take into account the shape as well as the location of the polymer.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Theory of polymer translocation through a flickering nanopore under an alternating driving force

We develop a theory for polymer translocation driven by a time-dependent force through an oscillating nanopore. To this end, we extend the iso-flux tension propagation theory (IFTP) [Sarabadani \textit{et al., J. Chem. Phys.}, 2014, \textbf{141}, 214907] for such a setup. We assume that the external driving force in the pore has a component oscillating in time, and the flickering pore is similarly described by an oscillating term in the pore friction. In addition to numerically solving the model, we derive analytical approximations that are in good agreement with the numerical simulations. Our results show that by controlling either the force or pore oscillations, the translocation process can be either sped up or slowed down depending on the frequency of the oscillations and the characteristic time scale of the process. We also show that while in the low and high frequency limits the translocation time $τ$ follows the established scaling relation with respect to chain length $N_0$, in the intermediate frequency regime small periodic fluctuations can have drastic effects on the dynamical scaling. The results can be easily generalized for non-periodic oscillations and elucidate the role of time dependent forces and pore oscillations in driven polymer translocation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Transition State Theory Approach to Polymer Escape from a One Dimensional Potential Well

The rate of escape of an ideal bead-spring polymer in a symmetric double-well potential is calculated using transition state theory (TST) and the results compared with direct dynamical simulations. The minimum energy path of the transitions becomes flat and the dynamics diffusive for long polymers making the Kramers-Langer estimate poor. However, TST with dynamical corrections based on short time trajectories started at the transition state gives rate constant estimates that agree within a factor of two with the molecular dynamics simulations over a wide range of bead coupling constants and polymer lengths. The computational effort required by the TST approach does not depend on the escape rate and is much smaller than that required by molecular dynamics simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Biopolymer filtration in corrugated nanochannels

We examine pressure-driven non-equilibrium transport of linear, circular and star polymers through a nanochannel containing a rectangular pit with full hydrodynamic interactions and thermal fluctuations. We demonstrate that with sufficiently small pressure differences, there is contour length-dependent entropic trapping of the polymer in the pit when the pit and the polymer sizes are compatible. This is due to competition between flow and chain relaxation in the pit, which leads to a non-monotonic dependence of the polymer mobility on its size and should aid in the design of nanofiltration devices based on the polymer size and shape.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Controlling polymer translocation and ion transport via charge correlations

We develop a correlation-corrected transport theory in order to predict ionic and polymer transport properties of membrane nanopores in physical conditions where mean-field electrostatics breaks down. The experimentally observed low KCl conductivity of open alpha-Hemolysin pores is quantitatively explained by the presence of surface polarization effects. Upon the penetration of a DNA molecule into the pore, these polarization forces combined with the electroneutrality of DNA sets a lower boundary for the ionic current, explaining the weak salt dependence of blocked pore conductivities at dilute ion concentrations. The addition of multivalent counterions into the solution results in the reversal of the polymer charge and the direction of the electroosmotic flow. With trivalent spermidine or quadrivalent spermine molecules, the charge inversion is strong enough to stop the translocation of the polymer and to reverse its motion. This mechanism can be used efficiently in translocation experiments in order to improve the accuracy of DNA sequencing by minimizing the translocation velocity of the polymer.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Iso-Flux Tension Propagation Theory of Driven Polymer Translocation: The Role of Initial Configurations

We investigate the dynamics of pore-driven polymer translocation by theoretical analysis and molecular dynamics (MD) simulations. Using the tension propagation theory within the constant flux approximation we derive an explicit equation of motion for the tension front. From this we derive a scaling relation for the average translocation time $τ$, which captures the asymptotic result $τ\propto N_0^{1+ν}$, where $N_0$ is the chain length and $ν$ is the Flory exponent. In addition, we derive the leading correction-to-scaling term to $τ$ and show that all terms of order $N_0^{2ν}$ exactly cancel out, leaving only a finite-chain length correction term due to the effective pore friction, which is linearly proportional to $N_0$. We use the model to numerically include fluctuations in the initial configuration of the polymer chain in addition to thermal noise. We show that when the {\it cis} side fluctuations are properly accounted for, the model not only reproduces previously known results but also considerably improves the estimates of the monomer waiting time distribution and the time evolution of the translocation coordinate $s(t)$, showing excellent agreement with MD simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Alteration of gas phase ion polarizabilities upon hydration in high dielectric liquids

We investigate the modification of gas phase ion polarizabilities upon solvation in polar solvents and ionic liquids. To this aim, we develop a classical electrostatic theory of charged liquids composed of solvent molecules modeled as finite size dipoles, and embedding polarizable ions that consist of Drude oscillators. In qualitative agreement with ab-initio calculations of polar solvents and ionic liquids, the hydration energy of a polarizable ion in both type of dielectric liquid is shown to favor the expansion of its electronic cloud. Namely, the ion carrying no dipole moment in the gas phase acquires a dipole moment in the liquid environment, but its electron cloud also reaches an enhanced rigidity. We find that the overall effect is an increase of the gas phase polarizability upon hydration. In the specific case of ionic liquids, it is shown that this hydration process is driven by a collective solvation mechanism where the dipole moment of a polarizable ion induced by its interaction with surrounding ions self-consistently adds to the polarization of the liquid, thereby amplifying the dielectric permittivity of the medium in a substantial way. We propose this self-consistent hydration as the underlying mechanism behind the high dielectric permittivities of ionic liquids composed of small charges with negligible gas phase dipole moment. Hydration being a correlation effect, the emerging picture indicates that electrostatic correlations cannot be neglected in polarizable liquids.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Microscopic formulation of non-local electrostatics in polar liquids embedding polarizable ions

Non-local electrostatic interactions associated with the finite solvent size and ion polarizability are investigated within the mean-field linear response theory. To this end, we introduce a field theoretic model of a polar liquid composed of linear multipole solvent molecules and embedding polarizable ions modeled as Drude oscillators. Unlike previous dipolar Poisson-Boltzmann formulations treating the solvent molecules as point dipoles, our model is able to qualitatively reproduce the non-local dielectric response behavior of polar liquids observed in Molecular Dynamics simulations and Atomic Force Microscope experiments for water solvent at charged interfaces. The present theory explains the formation of the associated interfacial hydration layers in terms of a cooperative dipolar response mechanism driven by the reaction of the solvent molecules to their own polarization field. We also incorporate into the theory the relative multipole/dipole moments of water molecules obtained from quantum mechanical calculations, and show that the multipolar contributions to the dielectric permittivity are largely dominated by the dipolar one. We find that this stems from the mutual cancellation of the first two interfacial hydration layers of opposite net charge for multipolar liquids. Within the same non-local dielectric response theory, we show that the induced ion polarizability reverses the interfacial ion density trends predicted by the Poisson-Boltzmann theory, resulting in a surface affinity of coions and exclusion of counterions. The results indicate that the consideration of the discrete charge composition of solvent molecules and ions is the key step towards a microscopic understanding of non-local electrostatic effects in polar solvents.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Polymer Escape from a Metastable Kramers potential: Path Integral Hyperdynamics Study

We study the dynamics of flexible, semiflexible, and self-avoiding polymer chains moving under a Kramers metastable potential. Due to thermal noise, the polymers, initially placed in the metastable well, can cross the potential barrier, but these events are extremely rare if the barrier is much larger than thermal energy. To speed up the slow rate processes in computer simulations, we extend the recently proposed path integral hyperdynamics method to the cases of polymers. We consider the cases where the polymers' radii of gyration are comparable to the distance between the well bottom and the barrier top. We find that, for a flexible polymer, the crossing rate ($\mathcal{R}$) monotonically decreases with chain contour length ($L$), but with the magnitude much larger than the Kramers rate in the globular limit. For a semiflexible polymer, the crossing rate decreases with $L$ but becomes nearly constant for large $L$. For a fixed $L$, the crossing rate becomes maximum at an intermediate bending stiffness. For a self-avoiding chain, the rate is a nonmonotonic function of $L$, first decreasing with $L$, and then, above certain length, increasing with $L$. These findings can be instrumental for efficient separation of biopolymers.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dipolar depletion effect on the differential capacitance of carbon based materials

The remarkably low experimental values of the capacitance data of carbon based materials in contact with water solvent needs to be explained from a microscopic theory in order to optimize the efficiency of these materials. We show that this experimental result can be explained by the dielectric screening deficiency of the electrostatic potential, which in turn results from the interfacial solvent depletion effect driven by image dipole interactions. We show this by deriving from the microscopic system Hamiltonian a non-mean-field dipolar Poisson-Boltzmann equation. This can account for the interaction of solvent molecules with their electrostatic image resulting from the dielectric discontinuity between the solvent medium and the substrate. The predictions of the extended dipolar Poisson-Boltzmann equation for the differential capacitance are compared with experimental data and good agreement is found without any fitting parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Electrostatic correlations in inhomogeneous charged fluids beyond loop expansion

Electrostatic correlation effects in inhomogeneous symmetric electrolytes are investigated within a previously developed electrostatic self-consistent (SC) theory (R.R. Netz and H. Orland, Eur. Phys.J. E 11, 301 (2003)). To this aim, we introduce two computational approaches that allow to solve the SC equations beyond the loop expansion. Both approaches can handle the case of dielectrically discontinuous boundaries where the one-loop theory is known to fail. By comparing the theoretical results obtained from these schemes with the results of the MC simulations that we ran for ions at neutral single dielectric interfaces as well as with previous MC data for charged interfaces, we first show that the weak coupling (WC) Debye-Huckel (DH) theory remains quantitatively accurate up to the bulk ion density rhob=0.01 M, whereas the SC theory exhibits a good quantitative accuracy up to rhob=0.2 M. Then, we derive from the perturbative SC scheme the one-loop theory of asymmetrically partitioned salt systems around a dielectrically homogeneous charged surface. It is shown that correlation effects originate in these systems from a competition between the salt screening loss at the interface driving the ions to the bulk region, and the interfacial counterion screening excess attracting them towards the surface. In the case of weak surface charges, the interfacial salt screening loss is the dominant effect. As a result, correlations decrease the MF density of both coions and counterions. With increasing surface charge, the surface-attractive counterion screening excess starts to dominate, and correlation effects amplify in this regime the MF density of both type of ions. We also show that at a characteristic value of the electrostatic coupling parameter, electrostatic correlations result in a charge inversion effect.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Hydrodynamic Effects on Confined Polymers

We consider the statics and dynamics of a flexible polymer confined between parallel plates both in the presence and absence of hydrodynamic interactions. The hydrodynamic interactions are described at the level of the fluctuating, compressible Navier-Stokes equation. We consider two cases: (i) confinement for both the solvent and the polymer, and (ii) confinement for the polymer only (in a 3D solvent), which is experimentally feasible, for instance, by (optical) trapping. We find a continuous transition from 2D to 3D dynamic scaling as a function of decreasing degree of confinement within the de Gennes and the weak-confinement regimes. We demonstrate that, in the presence of hydrodynamics, the polymer's center-of-mass diffusion coefficient in the direction parallel to the walls scales differently as a function of the level of confinement in cases (i) and (ii). We also find that in the commonly used Langevin dynamics description, the polymer swells more parallel to the walls than in the presence of hydrodynamics, and the planar diffusion coefficient shows scaling behavior similar to case (ii) rather than case (i). In addition, we quantify the differences in the static structure factor of the polymer between cases (i) and (ii), and between case (i) and Langevin dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Hydrodynamic forces on steady and oscillating porous particles

We derive new analytical results for the hydrodynamic force exerted on a sinusoidally oscillating porous shell and a sphere of uniform density in the Stokes limit. The coupling between the spherical particle and the solvent is done using the Debye-Bueche-Brinkman (DBB) model, i.e. by a frictional force proportional to the local velocity difference between the permeable particle and the solvent. We compare our analytical results and existing dynamic theories to Lattice-Boltzmann simulations of full Navier-Stokes equations for the oscillating porous particle. We find our analytical results to agree with simulations over a broad range of porosities and frequencies.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Influence of nanoparticle size, loading, and shape on the mechanical properties of polymer nanocomposites

We study the influence of spherical, triangular, and rod-like nanoparticles on the mechanical properties of a polymer nanocomposite (PNC), via coarse-grained molecular dynamics simulations. We focus on how the nanoparticle size, loading, mass, and shape influence the PNC's elastic modulus, stress at failure and resistance against cavity formation and growth, under external stress. We find that in the regime of strong polymer-nanoparticle interactions, the formation of a polymer network via temporary polymer-nanoparticle crosslinks has a predominant role on the PNC reinforcement. Spherical nanoparticles, whose size is comparable to that of the polymer monomers, are more effective at toughening the PNC than larger spherical particles. When comparing particles of spherical, triangular, and rod-like geometries, the rod-like nanoparticles emerge as the best PNC toughening agents.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Influence of non-universal effects on dynamical scaling in driven polymer translocation

We study the dynamics of driven polymer translocation using both molecular dynamics (MD) simulations and a theoretical model based on the non-equilibrium tension propagation on the {\it cis} side subchain. We present theoretical and numerical evidence that the non-universal behavior observed in experiments and simulations are due to finite chain length effects that persist well beyond the relevant experimental and simulation regimes. In particular, we consider the influence of the pore-polymer interactions and show that they give a major contribution to the non-universal effects. In addition, we present comparisons between the theory and MD simulations for several quantities, showing extremely good agreement in the relevant parameter regimes. Finally, we discuss the potential limitations of the present theories.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Scaling theory of driven polymer translocation

We present a theoretical argument to derive a scaling law between the mean translocation time $τ$ and the chain length $N$ for driven polymer translocation. This scaling law explicitly takes into account the pore-polymer interactions, which appear as a correction term to asymptotic scaling and are responsible for the dominant finite size effects in the process. By eliminating the correction-to-scaling term we introduce a rescaled translocation time and show, by employing both the Brownian Dynamics Tension Propagation theory [Ikonen {\it et al.}, Phys. Rev. E {\bf 85}, 051803 (2012)] and molecular dynamics simulations that the rescaled exponent reaches the asymptotic limit in a range of chain lengths that is easily accessible to simulations and experiments. The rescaling procedure can also be used to quantitatively estimate the magnitude of the pore-polymer interaction from simulations or experimental data. Finally, we also consider the case of driven translocation with hydrodynamic interactions (HIs). We show that by augmenting the BDTP theory with HIs one reaches a good agreement between the theory and previous simulation results found in the literature. Our results suggest that the scaling relation between $τ$ and $N$ is retained even in this case.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Unifying model of driven polymer translocation

We present a Brownian dynamics model of driven polymer translocation, in which non-equilibrium memory effects arising from tension propagation (TP) along the cis side subchain are incorporated as a time-dependent friction. To solve the effective friction, we develop a finite chain length TP formalism, expanding on the work of Sakaue [Sakaue, PRE 76, 021803 (2007)]. The model, solved numerically, yields results in excellent agreement with molecular dynamics simulations in a wide range of parameters. Our results show that non-equilibrium TP along the cis side subchain dominates the dynamics of driven translocation. In addition, the model explains the different scaling of translocation time w.r.t chain length observed both in experiments and simulations as a combined effect of finite chain length and pore-polymer interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Polymer Translocation Induced by a Bad Solvent

We employ 3D Langevin Dynamics simulations to study the dynamics of polymer chains translocating through a nanopore in presence of asymmetric solvent conditions. Initially a large fraction ($>$ 50%) of the chain is placed at the \textit{cis} side in a good solvent while the $trans$ segments are placed in a bad solvent that causes the chain to collapse and promotes translocation from the $cis$ to the $trans$ side. In particular, we study the ratcheting effect of a globule formed at the \textit{trans} side created by the translocated segment, and how this ratchet drives the system towards faster translocation. Unlike in the case of unbiased or externally forced translocation where the mean first passage time $\langle τ\rangle $ is often characterized by algebraic scaling as a function of the chain length $N$ with a single scaling exponent $α$, and the histogram of the mean first passage time $P(τ/\langleτ\rangle)$ exhibits scaling, we find that scaling is not well obeyed. For relatively long chains we find $\langle τ\rangle \sim N^α$ where $α\approx 1$ for $\varepsilon/k_{B}T > 1$. In this limit, we also find that translocation proceeds with a nearly constant velocity of the individual beads(monomers), which is attributed to the coiling of the globule. We provide an approximate theory assuming rotat ional motion restricted on a 2D disc to demonstrate that there is a crossover from diffusive behavior of the center of mass for short chains to a single file translocation for long chains, where the average translocation time scales linearly with the chain length $N$.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Driven polymer translocation through nanopores: slow versus fast dynamics

We investigate the dynamics of polymer translocation through nanopores under external driving by 3D Langevin Dynamics simulations, focusing on the scaling of the average translocation time $τ$ versus the length of the polymer, $τ\sim N^α$. For slow translocation, i.e., under low driving force and/or high friction, we find $α\approx 1+ν\approx 1.588$ where $ν$ denotes the Flory exponent. In contrast, $α\approx 1.37$ is observed for fast translocation due to the highly deformed chain conformation on the trans side, reflecting a pronounced non-equilibrium situation. The dependence of the translocation time on the driving force is given by $τ\sim F^{-1}$ and $τ\sim F^{-0.80}$ for slow and fast translocation, respectively. These results clarify the controversy on the magnitude of the scaling exponent $α$ for driven translocation.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Three-dimensional "Mercedes-Benz" model for water

In this paper we introduce a three-dimensional version of the Mercedes-Benz model to describe water molecules. In this model van der Waals interactions and hydrogen bonds are given explicitly through a Lennard-Jones potential and a Gaussian orientation-dependent terms, respectively. At low temperature the model freezes forming Ice-I and it reproduces the main peaks of the experimental radial distribution function of water. In addition to these structural properties, the model also captures the thermodynamical anomalies of water: the anomalous density profile, the negative thermal expansivity, the large heat capacity and the minimum in the isothermal compressibility.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Scaling exponents of Forced Polymer Translocation through a nano-pore

We investigate several scaling properties of a translocating homopolymer through a thin pore driven by an external field present inside the pore only using Langevin Dynamics (LD) simulation in three dimension (3D). Specifically motivated by several recent theoretical and numerical studies that are apparently at odds with each other, we determine the chain length dependence of the scaling exponents of the average translocation time, the average velocity of the center of mass, $<v_{CM}>$, the effective radius of gyration during the translocation process, and the scaling exponent of the translocation coordinate ($s$-coordinate) as a function of the translocation time. We further discuss the possibility that in the case of driven translocation the finite pore size and its geometry could be responsible that the veclocity scaling exponent is less than unity and discuss the dependence of the scaling exponents on the pore geometry for the range of $N$ studied here.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Polymer translocation through a nanopore under an applied external field

We investigate the dynamics of polymer translocation through a nanopore under an externally applied field using the 2D fluctuating bond model with single-segment Monte Carlo moves. We concentrate on the influence of the field strength $E$, length of the chain $N$, and length of the pore $L$ on forced translocation. As our main result, we find a crossover scaling for the translocation time $τ$ with the chain length from $τ\sim N^{2ν}$ for relatively short polymers to $τ\sim N^{1 + ν}$ for longer chains, where $ν$ is the Flory exponent. We demonstrate that this crossover is due to the change in the dependence of the translocation velocity v on the chain length. For relatively short chains $v \sim N^{- ν}$, which crosses over to $v \sim N^{- 1}$ for long polymers. The reason for this is that with increasing $N$ there is a high density of segments near the exit of the pore, which slows down the translocation process due to slow relaxation of the chain. For the case of a long nanopore for which $R_\parallel $, the radius of gyration $R_{g}$ along the pore, is smaller than the pore length, we find no clear scaling of the translocation time with the chain length. For large $N$, however, the asymptotic scaling $τ\sim N^{1 + ν}$ is recovered. In this regime, $τ$ is almost independent of $L$. We have previously found that for a polymer, which is initially placed in the middle of the pore, there is a minimum in the escape time for $R_\parallel \approx L$. We show here that this minimum persists for a weak fields $E$ such that $EL$ is less than some critical value, but vanishes for large values of $EL$.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Polymer translocation through a nanopore: a two-dimensional Monte Carlo simulation

We investigate the problem of polymer translocation through a nanopore in the absence of an external driving force. To this end, we use the two-dimensional (2D) fluctuating bond model with single-segment Monte Carlo moves. To overcome the entropic barrier without artificial restrictions, we consider a polymer which is initially placed in the middle of the pore, and study the escape time required for the polymer to completely exit the pore on either end. In particular, we examined the effect of the pore length on the escape time.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

Nucleation, growth, and scaling in slow combustion

We study the nucleation and growth of flame fronts in slow combustion. This is modeled by a set of reaction-diffusion equations for the temperature field, coupled to a background of reactants and augmented by a term describing random temperature fluctuations for ignition. We establish connections between this model and the classical theories of nucleation and growth of droplets from a metastable phase. Our results are in good argeement with theoretical predictions.

0 citations0 reviewsNo code linkedOpen access