Source author record

Mikko Karttunen

Mikko Karttunen appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

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

20 published item(s)

preprint2022arXiv

Epithelial Tissue Growth Dynamics: Universal or Not?

Universality of interfacial roughness in growing epithelial tissue has remained a controversial issue. Kardar-Parisi-Zhang (KPZ) and Molecular Beam Epitaxy (MBE) universality classes have been reported among other behaviors including total lack of universality. Here, we utilize a kinetic division model for deformable cells to investigate cell-colony scaling. With seemingly minor model changes, it can reproduce both KPZ- and MBE-like scaling in configurations that mimic the respective experiments. This result neutralizes the apparent scaling controversy. It can be speculated that this diversity in growth behavior is beneficial for efficient evolution and versatile growth dynamics.

preprint2022arXiv

Machine Learning-Driven Process of Alumina Ceramics Laser Machining

Laser machining is a highly flexible non-contact manufacturing technique that has been employed widely across academia and industry. Due to nonlinear interactions between light and matter, simulation methods are extremely crucial, as they help enhance the machining quality by offering comprehension of the inter-relationships between the laser processing parameters. On the other hand, experimental processing parameter optimization recommends a systematic, and consequently time-consuming, investigation over the available processing parameter space. An intelligent strategy is to employ machine learning (ML) techniques to capture the relationship between picosecond laser machining parameters for finding proper parameter combinations to create the desired cuts on industrial-grade alumina ceramic with deep, smooth and defect-free patterns. Laser parameters such as beam amplitude and frequency, scanner passing speed and the number of passes over the surface, as well as the vertical distance of the scanner from the sample surface, are used for predicting the depth, top width, and bottom width of the engraved channels using ML models. Owing to the complex correlation between laser parameters, it is shown that Neural Networks (NN) are the most efficient in predicting the outputs. Equipped with an ML model that captures the interconnection between laser parameters and the engraved channel dimensions, one can predict the required input parameters to achieve a target channel geometry. This strategy significantly reduces the cost and effort of experimental laser machining during the development phase, without compromising accuracy or performance. The developed techniques can be applied to a wide range of ceramic laser machining processes.

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.

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.

preprint2020arXiv

Manipulation of Diatomic Molecules with 'Oriented External Electric Fields': Linear Correlations in Atomic Properties Lead to Non-Linear Molecular Responses

'Oriented external electric fields (OEEFs)' have been shown to have great potential in being able to provide unprecedented control of chemical reactions, catalysis and selectivity with applications ranging from H2 storage to molecular machines. We report a theoretical study of the atomic origins of molecular changes due to OEEFs; understanding the characteristics of OEEF-induced couplings between atomic and molecular properties is an important step toward comprehensive understanding of the effects of strong external fields on molecular structure, stability, and reactivity. We focus on the atomic and molecular (bond) properties of a set of homo- (H2, N2, O2, F2, and Cl2) and hetero-diatomic (HF, HCl, CO, and NO) molecules under intense external electric fields in the context of quantum theory of atoms in molecules (QTAIM). It is shown that atomic properties (atomic charges and energies, and localization index) correlate linearly with the field strengths, but molecular properties (bond length, electron density at bond critical point, bond length, and electron delocalization index) exhibit non-linear responses to the imposed fields. In particular, the changes in the electron density distribution alter the shapes and locations of the zero-flux surfaces, atomic volumes, atomic electron population, and localization/delocalization indices. At the molecular level, the topography and topology of the molecular electrostatic potential undergo dramatic changes. The external fields also perturb the covalent-polar-ionic characteristic of the studied chemical bonds, hallmarking the impact of electric fields on the stability and reactivity of chemical compounds. The findings are well-rationalized within the framework of the quantum theory of atoms in molecules and form a coherent conceptual understanding of these effects in prototypical molecules such as diatomics.

preprint2016arXiv

Ab initio calculations of optical properties of silver clusters: Cross-over from molecular to nanoscale behavior

Electronic and optical properties of silver clusters were calculated using two different \textit{ab initio} approaches: 1) based on all-electron full-potential linearized-augmented plane-wave method and 2) local basis function pseudopotential approach. Agreement is found between the two methods for small and intermediate sized clusters for which the former method is limited due to its all-electron formulation. The latter, due to non-periodic boundary conditions, is the more natural approach to simulate small clusters. The effect of cluster size is then explored using the local basis function approach. We find that as the cluster size increases, the electronic structure undergoes a transition from molecular behavior to nanoparticle behavior at a cluster size of 140 atoms (diameter $\sim 1.7$\,nm). Above this cluster size the step-like electronic structure, evident as several features in the imaginary part of the polarizability of all clusters smaller than Ag$_\mathrm{147}$, gives way to a dominant plasmon peak localized at wavelengths 350\,nm$\leλ\le$ 600\,nm. It is, thus, at this length-scale that the conduction electrons' collective oscillations that are responsible for plasmonic resonances begin to dominate the opto-electronic properties of silver nanoclusters.

preprint2016arXiv

Getting excited: Challenges in quantum-classical studies of excitons in polymeric systems

A combination of classical molecular dynamics (MM/MD) and quantum chemical calculations based on the density functional theory (DFT) was performed to describe conformational properties of diphenylethyne (DPE), methylated-DPE and poly para phenylene ethynylene (PPE). DFT calculations were employed to improve and develop force field parameters for MM/MD simulations. Many-body Green's functions theory within the GW approximation and the Bethe-Salpeter equation were utilized to describe excited states of the systems. Reliability of the excitation energies based on the MM/MD conformations was examined and compared to the excitation energies from DFT conformations. The results show an overall agreement between the optical excitations based on MM/MD conformations and DFT conformations. This allows for calculation of excitation energies based on MM/MD conformations.

preprint2016arXiv

Molecular Dynamics simulations and Kelvin Probe Force microscopy to study of cholesterol-induced electrostatic nanodomains in complex lipid mixtures

The molecular arrangement of lipids and proteins within biomembranes and monolayers gives rise to complex film morphologies as well as regions of distinct electrical surface potential, topographical and electrostatic nanoscale domains. To probe these nanodomains in soft matter is a challenging task both experimentally and theoretically. This work addresses the effects of cholesterol, lipid composition, lipid charge, and lipid phase on the monolayer structure and the electrical surface potential distribution. Atomic Force Microscopy (AFM) was used to resolve topographical nanodomains and Kelvin Probe Force Microscopy (KPFM) to resolve electrical surface potential of these nanodomains in lipid monolayers. Model monolayers composed of dipalmitoylphosphatidylcholine (DPPC), 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), 1,2-dioleoyl-sn-glycero-3-[phospho-rac-(3-lysyl(1-glycerol))] (DOPG), sphingomyelin, and cholesterol were studied. It is shown that cholesterol changes nanoscale domain formation, affecting both topography and electrical surface potential. The molecular basis for differences in electrical surface potential was addressed with atomistic molecular dynamics (MD). MD simulations qualitatively match the experimental results, with 100s of mV difference in electrostatic potential between liquid-disordered bilayer (Ld, less cholesterol and lower chain order) and a liquid-ordered bilayer (Lo, more cholesterol and higher chain order). Importantly, the difference in electrostatic properties between Lo and Ld phases suggests a new mechanism by which membrane composition couples to membrane function.

preprint2016arXiv

Solvent effects on optical excitations of poly para phenylene ethynylene studied by QM/MM simulations based on Many-Body Green's Functions Theory

Electronic excitations in dilute solutions of poly para phenylene ethynylene (poly-PPE) are studied using a QM/MM approach combining many-body Green's functions theory within the $GW$ approximation and the Bethe-Salpeter equation with polarizable force field models. Oligomers up to a length of 7.5\,nm (10 repeat units) functionalized with nonyl side chains are solvated in toluene and water, respectively. After equilibration using atomistic molecular dynamics (MD), the system is partitioned into a quantum region (backbone) embedded into a classical (side chains and solvent) environment. Optical absorption properties are calculated solving the coupled QM/MM system self-consistently and special attention is paid to the effects of solvents. The model allows to differentiate the influence of oligomer conformation induced by the solvation from electronic effects related to local electric fields and polarization. It is found that the electronic environment contributions are negligible compared to the conformational dynamics of the conjugated PPE. An analysis of the electron-hole wave function reveals a sensitivity of energy and localization characteristics of the excited states to bends in the global conformation of the oligomer rather than to the relative of phenyl rings along the backbone.

preprint2016arXiv

The good, the bad and the user in soft matter simulations

Molecular dynamics (MD) simulations have become popular in materials science, biochemistry, biophysics and several other fields. Improvements in computational resources, in quality of force field parameters and algorithms have yielded significant improvements in performance and reliability. On the other hand, no method of research is error free. In this review, we discuss a few examples of errors and artifacts due to various sources and discuss how to avoid them. Besides bringing attention to artifacts and proper practices in simulations, we also aim to provide the reader with a starting point to explore these issues further. In particular, we hope that the discussion encourages researchers to check software, parameters, protocols and, most importantly, their own practices in order to minimize the possibility of errors. The focus here is on practical issues.

preprint2015arXiv

Crossovers in supercooled solvation water: Effects of hydrophilic & hydrophobic interactions

We use systematic 8 ns ab initio molecular dynamics (AIMD) to study the structure and dynamics of water in bulk, and close to both hydrophobic and hydrophilic (carbonyl) groups of tetramethylurea (TMU). We observe crossovers in the dynamical behavior around the hydrophobic group at $T_X = 256 \pm 4$K and another one at $265 \pm 5$K related to the relative strength of water-water and water-carbonyl hydrogen bonds (HBs). For bulk water, the temperature of the apparent divergence in relaxation times is located at $T_c = 210 \pm 10$K. To better identify the effects arising from the hydrophilic carbonyl group, systems of water with a methane molecule were used as references. These findings are related to the structural and thermodynamic transitions reported for proteins in solution and may also play a role in the mechanism of cold denaturation of proteins.

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.

preprint2014arXiv

Micelle fragmentation and wetting in confined flow

We use coarse-grained molecular-dynamics (MD) simulations to investigate the structural and dynamical properties of micelles under non-equilibrium Poiseuille flow in a nano-confined geometry. The effects of flow, confinement, and the wetting properties of die-channel walls on spherical sodium dodecyl sulfate (SDS) micelles are explored when the micelle is forced through a die-channel slightly smaller than its equilibrium size. Inside the channel, the micelle may fragment into smaller micelles. In addition to the flow rate, the wettability of the channel surfaces dictates whether the micelle fragments and determines the size of the daughter micelles: The overall behavior is determined by the subtle balance between hydrodynamic forces, micelle-wall interactions and self-assembly forces.

preprint2014arXiv

Modeling the Behavior of Confined Colloidal Particles Under Shear Flow

We investigate the behavior of colloidal suspensions with different volume fractions confined between parallel walls under a range of steady shears. We model the particles using molecular dynamics (MD) with full hydrodynamic interactions implemented through the use of a lattice-Boltzmann (LB) fluid. A quasi-2d ordering occurs in systems characterized by a coexistence of coupled layers with different densities, order, and granular temperature. We present a phase diagram in terms of shear and volume fraction for each layer, and demonstrate that particle exchange between layers is required for entering the disordered phase.

preprint2013arXiv

Phase Field Crystal Model for Magneto-Elasticity in Isotropic Ferromagnetic Solids

A new isotropic magneto-elastic phase field crystal (PFC) model to study the relation between morphological structure and magnetic properties of pure ferromagnetic solids is introduced. Analytic calculations were used to determine the phase diagram and obtain the relationship between elastic strains and magnetization. Time dependent numerical simulations were used to demonstrate the effect of grain boundaries on the formation of magnetic domains. It was shown that the grain boundaries act as nucleating sites for domains of reverse magnetization. Finally, we derive a relation for coercivity versus grain mis-orientation in the isotropic limit.

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.

preprint2012arXiv

Multiphase density functional theory parameterization of the Gupta potential for silver and gold

The ground state energies of Ag and Au in the face-centered cubic (FCC), body-centered cubic (BCC), simple cubic (SC) and the hypothetical diamond-like phase, and dimer were calculated as a function of bond length using density functional theory (DFT). These energies were then used to parameterize the many-body Gupta potential for Ag and Au. This parameterization over several phases of Ag and Au was performed to guarantee transferability of the potentials and to make them appropriate for studies of related nanostructures. Depending on the structure, the energetics of the surface atoms play a crucial role in determining the details of the nanostructure. The accuracy of the parameters was tested by performing a 2 ns MD simulation of a cluster of 55 Ag atoms -- a well studied cluster of Ag, the most stable structure being the icosahedral one. Within this time scale, the initial FCC lattice was found to transform to the icosahedral structure at room temperature. The new set of parameters for Ag was then used in a temperature dependent atom-by-atom deposition of Ag nanoclusters of up to 1000 atoms. We find a deposition temperature of 500 $\pm$50 K where low energy clusters are generated, suggesting an optimal annealing temperature of 500 K for Ag cluster synthesis.

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.

preprint2003arXiv

Crumpling of a stiff tethered membrane

first-principles numerical simulation model for crumpling of a stiff tethered membrane is introduced. In our model membranes, wrinkles, ridge formation, ridge collapse, as well as the initiation of stiffness divergence, are observed. The ratio of the amplitude and wave length of the wrinkles, and the scaling exponent of the stiffness divergence, are consistent with both theory and experiment. We observe that close to the stiffness divergence there appears a crossover beyond which the elastic behavior of a tethered membrane becomes similar to that of dry granular media. This suggests that ridge formation in membranes and force-chain network formation in granular packings are different manifestations of a single phenomenon.

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.