Abstract: The nanoscale photothermal effect and the optofluidic convection around plasmonic nanoparticles drive the application of such nanoparticles in micro-environment. In this work, heat transfer and fluid flow around Au nanospheres and nanorods in water medium under continuous and pulsed wave laser irradiance was investigated using an FEM based numerical framework. Au nanospheres of a wide range of diameter: 40 nm = Diameter (D) = 180 nm and relatively large nanorods (diameter: 50 nm) with varying aspect ratio (1 = Aspect ratio (A) = 5) and orientation (0 degrees = ? = 90 degrees, ? = 0 degrees, 90 degrees) with respect to the incident EM radiation were investigated for continuous wave (CW) and pulsed wave laser. It was found that although nanorods can attain much higher temperature than nanospheres, orientation of a nanorod is an important factor to be carefully considered in applications. In micro-scale spherical and hemispherical confinements (diameter < 14.4 p.m), the convective velocity fields around nanoparticles is in the order of 10-9 m/s, with only a weak effect of the slip or no-slip boundary condition on the confining walls. Importantly, the size of the confinement has a strong effect leading to an order of magnitude stronger convection for 14.4 p.m (diameter) spherical confinement as compared to 3.6 p.m confinement. Additionally close proximity of the nanoparticles to the confining walls strongly reduces (by an order of magnitude) the convective currents. The results reported herein provides important insights for the use of photothermal nanoparticles in microscale confined space (e.g. cellular environment) for applications such as optical tweezers, photoporation, etc.

Abstract: Non-equilibrium inductively coupled plasmas (ICPs) operating in hydrogen are of significant interest for applications including large-area materials processing. Increasing control of spatial gas heating, which drives the formation of neutral species density gradients and the rate of gas-temperature-dependent reactions, is critical. In this study, we use 2D fluid-kinetic simulations with the Hybrid Plasma Equipment Model to investigate the spatially resolved production of atomic hydrogen in a low-pressure planar ICP operating in pure hydrogen (10-20 Pa or 0.075-0.15 Torr, 300 W). The reaction set incorporates self-consistent calculation of the spatially resolved gas temperature and 14 vibrationally excited states. We find that the formation of neutral-gas density gradients, which result from spatially non-uniform electrical power deposition at constant pressure, can drive significant variations in the vibrational distribution function and density of atomic hydrogen when gas heating is spatially resolved. This highlights the significance of spatial gas heating on the production of reactive species in relatively high-power-density plasma processing sources.

Abstract: Plasma medicine has attracted tremendous interest in a variety of medical conditions, ranging from wound healing to antimicrobial applications, even in cancer treatment, through the interactions of cold atmospheric plasma (CAP) and various biological tissues directly or indirectly. The underlying mechanisms of CAP treatment are still poorly understood although the oxidative effects of CAP with amino acids, peptides, and proteins have been explored experimentally. In this study, machine learning (ML) technology is introduced to efficiently unveil the interaction mechanisms of amino acids and reactive oxygen species (ROS) in seconds based on the data obtained from the reactive molecular dynamics (MD) simulations, which are performed to probe the interaction of five types of amino acids with various ROS on the timescale of hundreds of picoseconds but with the huge computational load of several days. The oxidative reactions typically start with H-abstraction, and the details of the breaking and formation of chemical bonds are revealed; the modification types, such as nitrosylation, hydroxylation, and carbonylation, can be observed. The dose effects of ROS are also investigated by varying the number of ROS in the simulation box, indicating agreement with the experimental observation. To overcome the limits of timescales and the size of molecular systems in reactive MD simulations, a deep neural network (DNN) with five hidden layers is constructed according to the reaction data and employed to predict the type of oxidative modification and the probability of occurrence only in seconds as the dose of ROS varies. The well-trained DNN can effectively and accurately predict the oxidative processes and productions, which greatly improves the computational efficiency by almost ten orders of magnitude compared with the reactive MD simulation. This study shows the great potential of ML technology to efficiently unveil the underpinning mechanisms in plasma medicine based on the data from reactive MD simulations or experimental measurements. In this study, since reactive molecular dynamics simulation can currently only describe interactions between a few hundred atoms in a few hundred picoseconds, deep neural networks (DNN) are introduced to enhance the simulation results by predicting more data efficiently. image

Abstract: Promising biomedical uses for nonthermal plasma (NTP) in the fields of regenerative medicine, cancer therapy, and vaccine delivery involve the noninvasive application of uniform nonequilibrium plasma (including dielectric barrier discharge plasma) to living skin. Whereas most investigations have focused on achieving desired therapeutic outcomes, fewer studies have examined the mechanisms and pathways by which epithelial cells respond to NTP exposure. Using a transwell apical-basolateral-chambered system to culture the human keratinocyte HaCaT cell line, in vitro experiments were performed to demonstrate the effects of nanosecond-pulsed dielectric barrier discharge (nsDBD) plasma on polarized epithelial cell viability, monolayer permeability, intracellular oxidative stress, and the release of adenosine triphosphate (ATP). Application of nsDBD plasma at 60 Hz or below had minimal or no effect on HaCaT monolayer viability or permeability. nsDBD plasma exposure did, however, result in frequency-dependent reductions in intracellular glutathione (indicating direct induction of oxidative stress by nsDBD plasma) and increased extracellular ATP concentrations in the ba-solateral (subepithelial) media, which are indicators of cellular stress and an NTP-induced inflammatory response. These studies provide new insights into nsDBD plasma-induced inflammation and local innate immune responses initiated by polarized epithelial tissues.

Abstract: Studies using xenograft mouse models have shown that plasma applied to the skin overlying tumors results in tumor shrinkage. Plasma is considered a nonpenetrating treatment; however, these studies demonstrate plasma effects that occur beyond the postulated depth of physical penetration of plasma components. The present study examines the propagation of plasma effects through a tissue model using three-dimensional, cell-laden extracellular matrices (ECMs). These ECMs are used as barriers against direct plasma penetration. By placing them onto a monolayer of target cancer cells to create an in-vitro analog to in-vivo studies, we distinguished between cellular effects from direct plasma exposure and cellular effects due to cell-to-cell signaling stimulated by plasma. We show that nanosecond-pulsed dielectric barrier discharge plasma treatment applied atop an acellular barrier impedes the externalization of calreticulin (CRT) in the target cells. In contrast, when a barrier is populated with cells, CRT externalization is restored. Thus, we demonstrate that plasma components stimulate signaling among cells embedded in the barrier to transfer plasma effects to the target cells.

Abstract: Many geometries of plant organs can be described by the Gielis equation, a polar coordinate equation extended from the superellipse equation, . Here, r is the polar radius corresponding to the polar angle φ; m is a positive integer that determines the number of angles of the Gielis curve when φ ∈ [0 to 2π); and the rest of the symbols are parameters to be estimated. The pentagonal radial symmetry of calyxes and corolla tubes in top view is a common feature in the flowers of many eudicots. However, prior studies have not tested whether the Gielis equation can depict the shapes of corolla tubes. We sampled randomly 366 flowers of Vinca major L., among which 360 had five petals and pentagonal corolla tubes, and six had four petals and quadrangular corolla tubes. We extracted the planar coordinates of the outer rims of corolla tubes (in top view) (ORCTs), and then fitted the data with two simplified versions of the Gielis equation with k = 1 and m = 5: (Model 1), and (Model 2). The adjusted root mean square error (RMSEadj) was used to evaluate the goodness of fit of each model. In addition, to test whether ORCTs are radially symmetrical, we correlated the estimates of n2 and n3 in Model 1 on a log-log scale. The results validated the two simplified Gielis equations. The RMSEadj values for all corolla tubes were smaller than 0.05 for both models. The numerical values of n2 and n3 were demonstrated to be statistically equal based on the regression analysis, which suggested that the ORCTs of V. major are radially symmetrical. It suggests that Model 1 can be replaced by the simpler Model 2 for fitting the ORCT in this species. This work indicates that the pentagonal or quadrangular corolla tubes (in top view) can both be modeled by the Gielis equation and demonstrates that the pentagonal or quadrangular corolla tubes of plants tend to form radial symmetrical geometries during their development and growth.

Abstract: Measuring the inequality of leaf area distribution per plant (ILAD) can provide a useful tool for quantifying the influences of intra- and interspecific competition, foraging behavior of herbivores, and environmental stress on plants’ above-ground architectural structures and survival strategies. Despite its importance, there has been limited research on this issue. This paper aims to fill this gap by comparing four inequality indices to measure ILAD, using indices for quantifying household income that are commonly used in economics, including the Gini index (which is based on the Lorenz curve), the coefficient of variation, the Theil index, and the mean log deviation index. We measured the area of all leaves for 240 individual plants of the species Shibataea chinensis Nakai, a drought-tolerant landscape plant found in southern China. A three-parameter performance equation was fitted to observations of the cumulative proportion of leaf area vs. the cumulative proportion of leaves per plant to calculate the Gini index for each individual specimen of S. chinensis. The performance equation was demonstrated to be valid in describing the rotated and right shifted Lorenz curve, given that >96% of root-mean-square error values were smaller than 0.004 for 240 individual plants. By examining the correlation between any of the six possible pairs of indices among the Gini index, the coefficient of variation, the Theil index, and the mean log deviation index, the data show that these indices are closely related and can be used interchangeably to quantify ILAD.

Abstract: Climate models suggest that the persistence of summer precipitation regimes (PRs) is on the rise, characterized by both longer dry and longer wet durations. These PR changes may alter plant biochemical composition and thereby their economic and ecological characteristics. However, impacts of PR persistence have primarily been studied at the community level, largely ignoring the biochemistry of individual species. Here, we analyzed biochemical components of four grassland species with varying sensitivity to PR persistence (Holcus lanatus, Phleum pratense, Lychnis flos-cuculi, Plantago lanceolata) along a range of increasingly persistent PRs (longer consecutive dry and wet periods) in a mesocosm experiment. The more persistent PRs decreased nonstructural sugars, whereas they increased lignin in all species, possibly reducing plant quality. The most sensitive species Lychnis seemed less capable of altering its biochemical composition in response to altered PRs, which may partly explain its higher sensitivity. The more tolerant species may have a more robust and dynamic biochemical network, which buffers the effects of changes in individual biochemical components on biomass. We conclude that the biochemical composition changes are important determinants for plant performance under increasingly persistent precipitation regimes.

Abstract: Superimposition of dual frequencies (DFs) is one of the methods used for controlling plasma distribution in an inductively coupled plasma (ICP) source. The effects of a superimposed DF on the argon plasma characteristics have been investigated using a two-dimensional self-consistent fluid model. When both currents are fixed at 6A, the plasma density drops with decrease in one of the source frequencies due to less efficient heating and the plasma uniformity improves significantly. Moreover, for ICP operated with superimposed DFs (i.e., 4.52MHz/13.56MHz and 2.26MHz/13.56MHz), the current source exhibits the same period as the low frequency (LF) component, and the plasma density is higher than that obtained at a single frequency (i.e., 4.52 and 2.26MHz) with the same total current of 12A. However, at superimposed current frequencies of 6.78MHz/13.56MHz, the plasma density is lower than that obtained at a single frequency of 6.78MHz due to the weaker negative azimuthal electric field between two positive maxima during one period of 6.78MHz. When the superimposed DF ICP operates at 2.26 and 13.56MHz, the rapid oscillations of the induced electric field become weaker during one period of 2.26MHz as the current ratio of 2.26MHz/13.56MHz rises from 24A/7 A to 30A/1 A, and the plasma density drops with the current ratio due to weakened electron heating. The uniformity of plasma increases due to sufficient diffusion under the low-density condition.

Abstract: We study the electronic spectral properties at zero temperature of the one-dimensional (1D) version of the degenerate two-orbital Kanamori-Hubbard model, one of the well-established frameworks to study transition metal compounds, using state-of-the-art numerical techniques based on the density matrix renormalization group. While the system is Mott insulating for the half-filled case, as expected for an interacting 1D system, we find interesting and rich structures in the single-particle density of states (DOS) for the hole-doped system. In particular, we find the existence of in-gap states which are pulled down to lower energies from the upper Hubbard band with increasing the interorbital Coulomb interaction V. We analyze the composition of the DOS by projecting it onto different local excitations, and we observe that for large dopings these in-gap excitations are formed mainly by interorbital holon-doublon (HD) states and their energies follow approximately the HD states in the atomic limit. We observe that the Hund interaction J increases the width of the in-gap band, as expected from the two-particle fluctuations in the Hamiltonian. The observation of a finite density of states within the gap between the Hubbard bands for this extended 1D model indicates that these systems present a rich excitation spectra which could help us understand the microscopic physics behind multiorbital compounds.