Search NASA⌕ Search

SEARCH · Search NASA

Results for “interaction potential”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Symmetrized Drude Oscillator Force Fields Improve Numerical Performance of Polarizable Molecular Dynamics

Drude oscillator potentials are a popular and computationally efficient class of polarizable models that represent each polarizable atom as a positively charged Drude core harmonically bound to a negatively charged Drude shell. Here we show that existing force fields that place all non-Coulomb forces on the Drude core and none on the shell inadvertently couple the dipole to non-Coulombic forces. This introduces errors where interactions with neutral particles can erroneously induce atomic polarization, leading to spurious polarizations in the absence of an electric field, exacerbating violations of equipartition in the employed Carr–Parinello scheme. A suitable symmetrization of the interaction potential that correctly splits the force between the Drude core and shell can correct this shortcoming, improving the stability and numerical performance of Drude oscillator-based simulations. The symmetrization procedure is straightforward and only requires the rescaling of a few force field parameters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Twisted bilayer graphene as topological heavy fermion: II. Analytical approximations of the model parameters

The recently-introduced topological heavy fermion (THF) model [1] of twisted bilayer graphene (TBG) aims to reconcile the quantum-dot-like electronic structure of the latter observed by scanning tunneling microscopy, with its electron delocalization seen in transport measurements. The THF model achieves this by coupling localized (heavy) fermions with anomalous conduction electrons. Originally, the parameters of the THF model were obtained numerically from the Bistritzer–Macdonald (BM) model of TBG [1]. In this work, we derive analytical expressions for the THF model parameters as a function of the twist angle, the ratio between the tunneling amplitudes at the AA and AB regions (w 0 /w 1 ), and the screening length of the interaction potential. By numerically computing the THF model parameters across an extensive experimentally-relevant parameter space, we show that the resulting approximations are remarkably good, i.e., within the 30% relative error for almost the entire parameter space. At the single-particle level, the THF model accurately captures the energy spectrum of the BM model over a large phase space of angles and tunneling amplitude ratios. When interactions are included, we also show that the THF description of TBG is good around the magic angle for realistic values of the tunneling amplitude ratios (0.6 ≤ w 0 /w 1 ≤ 1.0), for which the hybridization between the localized and conduction fermions γ is smaller than the onsite repulsion of the heavy fermions U 1 (i.e., |γ| < U 1 ).

Approximation methods↗
DOE OSTI2020

Predictive Model of Charge Mobilities in Organic Semiconductor Small Molecules with Force-Matched Potentials

Charge mobility of crystalline organic semiconductors (OSC) is limited by local dynamic disorder. Recently, the charge mobility for several high mobility OSCs, including TIPS-pentacene, were accurately predicted from a density functional theory (DFT) simulation constrained by the crystal structure and the inelastic neutron scattering spectrum, which provide direct measures of the structure and the dynamic disorder in the length scale and energy range of interest. However, the computational expense required for calculating all of the atomic and molecular forces is prohibitive. In this work, we demonstrate the use of density functional tight binding (DFTB), a semiempirical quantum mechanical method that is 2 to 3 orders of magnitude more efficient than DFT. We show that force matching a many-body interaction potential to DFT derived forces yields highly accurate DFTB models capable of reproducing the low-frequency intricacies of experimental inelastic neutron scattering (INS) spectra and accurately predicting charge mobility. We subsequently predicted charge mobilities from our DFTB model of a number of previously unstudied structural analogues to TIPS-pentacene using dynamic disorder from DFTB and transient localization theory. The approach we establish here could provide a truly rapid simulation pathway for accurate materials properties prediction, in our vision applied to new OSCs with tailored properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Accurate calculation of many-body energies in water clusters using a classical geometry-dependent induction model

We incorporate geometry-dependent distributed multipole and polarizability surfaces into an induction model that is used to describe the 3- and 4-body terms of the interaction between water molecules. The expansion is carried out up to hexadecapole with the multipoles distributed on the atom sites. Dipole-dipole, dipole-quadrupole, and quadrupole-quadrupole distributed polarizabilities are used to represent the response of the multipoles to an electric field. We compare the model against two large databases consisting of 43,844 3-body terms and 3,603 4-body terms obtained from high level ab initio calculations previously used to fit the MB-pol and q-AQUA interaction potentials. The classical induction model with no adjustable parameters reproduces the ab-initio 3- and 4-body terms contained in these teo Databases with a Root-Mean-Square-Error (RMSE) of 0.104/0.058 and a Mean-Absolute-Error (MAE) of 0.054/0.026 kcal/mol, respectively, results that are on a par with those obtained 1 by fitting the same data using tens of thousands of Permutationally Invariant Polynomials (PIPs). This demonstrates the accuracy of this physically motivated model in describing the 3- and 4-body terms in the interactions between water molecules with no adjustable parameters. The triple-dipole-dispersion energy was included in the 3-body energy and was found to be small but not quite negligible. The model represents a practical, efficient and transferable approach for obtaining accurate non-additive interactions for multi-component systems without the need of performing tens of thousands of high level electronic structure calculations and fitting them with tens of thousands of PIPs.

molecular analysis water, aqueous↗
DOE OSTI2025

Algae Asphalt to Enhance Pavement Sustainability and Performance at Subzero Temperatures

This paper evaluates the potential of algae-derived biobinders as sustainable alternatives for pavement construction. It specifically examines the physicochemical and rheological properties of biomodified binders and their potential to offset carbon emissions when used as partial replacements for conventional petroleum-based asphalt binders. Biosequestration of CO 2 using microalgal cell factories is a promising way of recycling CO 2 into biomass via photosynthesis. Our study demonstrates that incorporating algae-derived binders into asphalt can significantly reduce carbon emissions. Each 1% increase in algae-based biobinder leads to an approximate 4.5% decrease in net carbon emissions. This indicates that a blend containing about 22% biobinder has the potential to achieve carbon neutrality. Blends with higher proportions may even result in net-negative emissions, highlighting a promising strategy for environmentally responsible road construction. In terms of performance, the study shows that certain algae-derived biobinders significantly enhance the cracking resistance of asphalt, particularly under subzero temperatures, by improving its stress-relief capacity. A key contribution of this work is the introduction of polarizability as a novel molecular-level parameter for assessing the compatibility of algae-derived bio-oils with asphalt. By capturing the electronic responsiveness of bio-oil molecules, polarizability serves as a predictive indicator of their interaction potential with asphalt components, providing a new dimension for evaluating the binder performance at the molecular scale. Among the tested materials, the biobinder derived from Haematococcus pluvialis demonstrated particularly strong improvements in resistance to permanent deformation under repeated loading conditions analogous to traffic-induced stress, as well as enhanced resistance to moisture-induced damage. In conclusion, these findings advance the chemistry-driven design of biomass-based binders and highlight a promising pathway toward the development of low-carbon, high-performance, and sustainable infrastructure materials.

Algae↗
DOE OSTI2020

Composite Analysis Solid Waste Release Data Reduction of Activity Flux from Waste Sites to the Vadose Zone for Baseline Assessment

This environmental calculation file (ECF) describes the data reduction methodology for the solid waste release datasets generated by the Composite Analysis Solid Waste Release (CASWR) model (CP-62766, Model Package Report: Composite Analysis Solid Waste Release Model [CASWR Model]) and documents the resulting reduced datasets for the updated Composite Analysis (CA) baseline assessment. This ECF is in support of the updated Hanford Site CA. The purpose of the updated Hanford Site CA is to provide an estimate of the cumulative radiological impacts from active and planned low-level radioactive waste disposals and other potentially interacting radioactive waste disposal sources that will remain following Hanford Site closure. The vadose zone modeling activities associated with the CA includes solid waste forms as contaminant sources. To support the modeling of solid waste transport from the vadose zone to the groundwater at the Hanford Site, the CASWR model was designed to generate deterministic radionuclide release rates for the Hanford Site Central Plateau solid waste disposal sites. The complete CASWR model-generated dataset for the CA Baseline Assessment consists of 2,378,831 sets of time and release rates (1,073,869 sets for the 200 East Area and 1,304,962 sets for the 200 West Area), representing a course of the simulation from year 1953 to year 12,069, 16 radiological contaminants of potential concern (COPCs), and 237 solid waste sources as documented in ECFHANFORD- 19-0112, Solid Waste Release Calculations for the Composite Analysis Baseline Assessment. For several of the STOMP vadose zone models, the number of release rate data pairs required for release at all waste sites would exceed the limitations imposed by the STOMP modelling software package used to simulate contaminant transport through the vadose zone. As a result, a data reduction of the CASWR dataset is required prior to its use as a source input for the CA vadose zone modeling activities.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Impacts of abiotic and biotic factors on tundra productivity near Utqiaġvik, Alaska

Abstract Earlier snowmelt, warmer temperatures and herbivory are among the factors that influence high-latitude tundra productivity near the town of Utqiaġvik in northern Alaska. However, our understanding of the potential interactions between these factors is limited. MODIS observations provide cover fractions of vegetation, snow, standing water, and soil, and fractional absorption of photosynthetically active radiation by canopy chlorophyll (fAPAR chl ) per pixel. Here, we evaluated a recent time-period (2001–2014) that the tundra experienced large interannual variability in vegetation productivity metrics (i.e. fAPAR chl and APAR chl ), which was explainable by both abiotic and biotic factors. We found earlier snowmelt to increase soil and vegetation cover, and productivity in June, while warmer temperatures significantly increased monthly productivity. However, abiotic factors failed to explain stark decreases in productivity during August of 2008, which coincided with a severe lemming outbreak. MODIS observations found this tundra ecosystem to completely recover two years later, resulting in elevated productivity. This study highlights the potential roles of both climate and herbivory in modulating the interannual variability of remotely retrieved plant productivity metrics in Arctic coastal tundra ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Similarity Metric for Data Optimization and Efficient Training of Reactive Machine Learning Force Fields for Hydrocarbon Radiolysis

Radiolysis is a common approach to sterilize polymers, chemically modify them for upcycling, and accelerate their decomposition for recycling purposes. Reactive molecular dynamics (MD) simulations provide a powerful tool to generate atomic-level trajectories of the reactive processes and quantify radiolytic chemical degradation pathways. For this, machine learning (ML) surrogate models for reactive force fields with quantum mechanical accuracy are now widely used, which require ML training data sets that can provide information on atomic environments for target chemical systems. However, radiolysis chemistry can be highly complex and diverse, which poses significant challenges for generating training data to parametrize ML models. In this regard, we developed a method for optimizing the training data set using a cosine similarity metric to help guide training set selection for radiolysis of polyethylene, a model hydrocarbon polymer, as well as to enhance the transferability of our reactive ML force field (MLFF) to a variety of molecular and polymeric systems. Our approach performs atom-by-atom comparisons between local atomic environments to pinpoint important data points associated with rare and localized events, such as radiolysis damage within structures. We apply this approach to train the Chebyshev Interaction Model for Efficient Simulation (ChIMES) MLFF model, which expresses the atomic interaction potentials in terms of linear combinations of many-body Chebyshev polynomials. We first show that our method can reduce our training set size by ∼70% while improving overall accuracy compared to more standard MD model fitting approaches. We then validate our optimum model against diverse hydrocarbon simulation data, including simple alkanes and systems with unsaturated carbon bonds, over a wide range of thermodynamic conditions. Finally, we use our ChIMES model to perform MD simulations of radiolytic damage with large-scale systems that help avoid system size effects. Overall, our approach yields an MD force field that retains most of the accuracy of the underlying quantum method while yielding many orders of improvement in computational efficiency. In conclusion, our efforts will have impact on future hydrocarbon polymer radiolysis studies, where the chemical details of the polymer–radiation interactions can have a strong effect on the resulting products observed in experiments.

Hydrocarbons↗
DOE OSTI2022

Driven superconducting vortex dynamics in systems with twofold anisotropy in the presence of pinning

Abstract We examine the dynamics of superconducting vortices with twofold anisotropic interaction potentials driven over random pinning, and compare the behavior under drives applied along the hard and the soft anisotropy directions. As the driving force increases, the number of topological defects reaches a maximum near the depinning threshold, and then decreases as the vortices form one-dimensional (1D) chains. This coincides with a transition from a pinned nematic to a moving smectic aligned with the soft anisotropy direction. The system is generally more ordered when the drive is applied along the soft direction of the anisotropy. For driving along the hard direction, there is a critical value of the twofold anisotropy above which the system remains aligned with the soft direction. Hysteretic behavior appears upon cycling the driving force, with 1D vortex chains persisting during the decreasing leg below the threshold for chain formation for increasing drive. More anisotropic systems have a greater amount of structural disorder in the moving state. For lower anisotropy, the system forms a moving smectic-A state, while at higher anisotropy, a moving nematic state appears instead.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Covalent labeling of the Arabidopsis plasma membrane H + ‐ ATPase reveals 3D conformational changes involving the C‐terminal regulatory domain

The plasma membrane proton pump is the primary energy transducing, electrogenic ion pump of the plasma membrane in plants and fungi. Compared to its fungal counterpart, the plant plasma membrane proton pump's regulatory C‐terminal domain (CTD) contains an additional regulatory segment that links multiple sensory pathways regulating plant cell length through phosphorylation and recruitment of regulatory 14‐3‐3 proteins. However, a complete structural model of a plant proton pump is lacking. Here, we performed covalent labeling with mass spectrometric analysis (CL‐MS) on the Arabidopsis pump AHA2 to identify potential interactions between the CTD and the catalytic domains. Our results suggest that autoinhibition in the plant enzyme is much more structurally complex than in the fungal enzyme.

Blackburn, Matthew R. [Department of Biochemistry ↗
DOE OSTI2020

Irrigation, carbon amelioration, nitrogen, and stover removal effects on continuous corn

Corn (Zea mays L.) residue or stover is harvested as supplemental feed for livestock and is a potential feedstock for cellulosic biofuels. Limited information is available on corn stover removal effects on grain yield under different irrigation rates, nitrogen (N) fertilizer rates, and practices to maintain soil carbon (C) and minimize soil erosion. Here, we evaluated potential interactions between irrigation rate (full or limited), C amelioration practices (cover crop, surface-applied manure, or no amelioration practice), fertilizer N management (125 or 200 kg N ha –1 ), and stover removal (residue removal or no residue removal), on no-till continuous corn grain yield located on a silt loam soil in south-central Nebraska (2011–2018). Stover removal increased (P = .0017) grain yield by 1.02 Mg ha –1 and grain N uptake compared with no stover removal. Manure increased corn grain yield by 899 kg ha –1 compared to either the control or cover crop with stover removal while C amelioration practices did not affect grain yield under stover retention. Grain N concentration was higher (P < .0001) with stover removal (13.2 g kg –1 ) than residue retained (12.7 g kg –1 ). Partial factor productivity (grain yield/N rate) was highest for manure treatments with residue removal. The farmgate N balance (N applied – N removed) from stover removal was negative, yet grain yield was not affected after eight growing seasons. Grain yield increased with stover removal and manure application resulting in a cost-effective C amelioration practice whereas incorporating a winter cover crop resulted in similar grain yield as the control.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Analysis of Dynamical Peculiarities in Nanoalloys at Subsystems Level: Dynamical Degrees of Freedom, Temperature Differences, and the Chameleon Effect

Abstract A novel analysis of the dynamical behavior of nanoalloy systems, as represented by model Ni/Al 13‐atom clusters, over a broad range of energies that cover the stage‐wise transition of the systems from their solid‐like to liquid‐like state is presented. Conceptually, the analysis is rooted in partitioning the systems into judiciously chosen subsystems and characterizing the latter in terms of subsystem‐specific dynamical descriptors that include dynamical degrees of freedom , root‐mean‐square bond‐length fluctuation , and element‐specific subsystem temperature . The analysis reveals a host of intriguing new peculiarities in the dynamical behavior of the Ni/Al 13‐mers, among which are what we call the chameleon effect and the difference in the temperatures of the Ni and Al subsystems at high energies, a difference that strongly depends on the cluster composition and also changes with energy. These do not have an analog in pure Ni 13 and Al 13 and are explained in terms of the coupled effects of the difference between the masses of the Ni and Al atoms (the mass effect) and of the difference in the anharmonicity of the overall interaction potential as experienced by the Ni and Al subsystems of the clusters (the potential effect).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Measurement of temperature of a dusty plasma from its configuration

A new concept called “configurational temperature” is introduced in the context of dusty plasma, where the temperature of the dust particles submerged in the plasma can be measured directly from the positional information of the individual dust particles and the interaction potential between the dust grains. This method does not require the velocity information of individual particles, which is a key parameter to measure the dust temperature in the conventional method. The technique is initially tested using two-dimensional (2D) OpenMP parallel molecular dynamics and Monte Carlo simulation and then compared with the temperature evaluated from experimental data. The experiment have been carried out in the Dusty Plasma Experimental (DPEx) device, where a 2D stationary plasma crystal of melamine formaldehyde particles is formed in the cathode sheath of a DC glow discharge argon plasma. The kinetic temperature of the dust is calculated using the standard particle image velocimetry technique at different pressures. An extended simulation result for the three-dimensional case is also presented, which can be employed for the temperature measurement of a three-dimensional dust crystal in laboratory devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

How Fisheries Biologists Can Facilitate the Clean Energy Transition

Abstract Fisheries and aquatic biologists play a critical role in creating environmentally protective hydropower flow requirements that govern flow timing, frequency, magnitude, and rate of change. Hydropower's role in the U.S. electrical grid is expected to evolve in response to increased wind and solar generation as hydropower will be called upon to quickly ramp up and down in response to changes in wind and solar generation. For this reason, hydropower is expected to have increased value as fossil generation is phased out, even as rapid flow fluctuations linked with hydropower flexibility may strand fish, alter habitat, and create unsafe recreational conditions. We face a new challenge in facilitating the renewable energy transition—designing environmental flow requirements that protect against the impacts of flow fluctuations while allowing adequate hydropower flexibility to support a stable grid. In this paper, we discuss hydropower environmental flow requirements, operational flexibility, and electrical grid stability, their potential interactions, and opportunities to align environmental and power system needs to support healthy ecosystems, multiple water uses, and decarbonization of the electric grid.

13 HYDRO ENERGY↗
DOE OSTI2025

Integrated structural model of the palladin–actin complex using XL ‐ MS , docking, NMR , and SAXS

Abstract Palladin is an actin‐binding protein that accelerates actin polymerization and is linked to the metastasis of several types of cancer. Previously, three lysine residues in an immunoglobulin‐like domain of palladin have been identified as essential for actin binding. However, it is still unknown where palladin binds to F‐actin. Evidence that palladin binds to the sides of actin filaments to facilitate branching is supported by our previous study showing that palladin was able to compensate for Arp2/3 in the formation of Listeria actin comet tails. Here, we used chemical crosslinking to covalently link palladin and F‐actin residues based on spatial proximity. Samples were then enzymatically digested, separated by liquid chromatography, and analyzed by tandem mass spectrometry. Peptides containing the crosslinks and specific residues involved were then identified for input to the HADDOCK docking server to model the most likely binding conformation. Small‐angle x‐ray scattering was used to provide further insight into palladin flexibility and the binding interface, and NMR spectra identified potential interactions between palladin's Ig domains. Our final structural model of the F‐actin:palladin complex revealed how palladin interacts with and stabilizes F‐actin at the interface between two actin monomers. Three actin residues that were identified in this study also appear commonly in the actin‐binding interface with other proteins such as myotilin, myosin, and tropomodulin. An accurate structural representation of the complex between palladin and actin extends our understanding of palladin's role in promoting cancer metastasis through the regulation of actin dynamics.

Sargent, Rachel [Department of Chemistry and Bioch↗
DOE OSTI2022

$T\overline{T}$ deformations and the width of fundamental particles

We provide a simple geometric meaning for deformations of so-called $T\overline{T}$ type in relativistic and non-relativistic systems. Deformations by the cross products of energy and momentum currents in integrable quantum field theories are known to modify the thermodynamic Bethe ansatz equations by a “CDD factor”. In turn, CDD factors may be interpreted as additional, fixed shifts incurred in scattering processes: a finite width added to the fundamental particles (or, if negative, to the free space between them). We suggest that this physical effect is a universal way of understanding $T\overline{T}$ deformations, both in classical and quantum systems. We first show this in non-relativistic systems, with particle conservation and translation invariance, using the deformation formed out of the densities and currents of particles and momentum. This holds at the level of the equations of motion, and for any interaction potential, integrable or not. We then argue, and show by similar techniques in free relativistic particle systems, that $T\overline{T}$ deformations of relativistic systems produce the equivalent phenomenon, accounting for length contractions. We also show that, in both the relativistic and non-relativistic cases, the width of particles is equivalent to a state-dependent change of metric, where the distance function discounts the particles’ widths, or counts the additional free space. This generalises and explains the known field-dependent coordinate change describing $T\overline{T}$ deformations. The results connect such deformations with generalised hydrodynamics, where the relations between scattering shifts, widths of particles and state-dependent changes of metric have been established.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Noncommutative gauge symmetry in the fractional quantum Hall effect

Abstract We show that a system of particles on the lowest Landau level can be coupled to a probe U(1) gauge field$$ \mathcal{A} $$ A μ in such a way that the theory is invariant under a noncommutative U(1) gauge symmetry. While the temporal component$$ \mathcal{A} $$ A 0 of the probe field is coupled to the projected density operator, the spatial components$$ \mathcal{A} $$ A i are best interpreted as quantum displacements, which distort the interaction potential between the particles. We develop a Seiberg-Witten-type map from the noncommutative U(1) gauge symmetry to a simpler version, which we call “baby noncommutative” gauge symmetry, where the Moyal brackets are replaced by the Poisson brackets. The latter symmetry group is isomorphic to the group of volume preserving diffeomorphisms. By using this map, we resolve the apparent contradiction between the noncommutative gauge symmetry, on the one hand, and the particle-hole symmetry of the half-filled Landau level and the presence of the mixed Chern-Simons terms in the effective Lagrangian of the fractional quantum Hall states, on the other hand. We outline the general procedure which can be used to write down effective field theories which respect the noncommutative U(1) symmetry.

Physics↗
DOE OSTI2022

Tea consumption and risk of bladder cancer in the Bladder Cancer Epidemiology and Nutritional Determinants (BLEND) Study: Pooled analysis of 12 international cohort studies

Tea has been shown to be associated with reduced risk of several diseases including cardiovascular diseases, stroke, metabolic syndrome, and obesity. However, the results on the relationship between tea consumption and bladder cancer are conflicting. This research aimed to assess the association between tea consumption and risk of bladder cancer using a pooled analysis of prospective cohort data. Individual data from 532,949 participants in 12 cohort studies, were pooled for analyses. Cox regression models stratified by study centre was used to estimate hazard ratios (HR) and corresponding 95% CIs. Fractional polynomial regression models were used to examine the dose–response relationship. A higher level of tea consumption was associated with lower risk of bladder cancer incidence (compared with no tea consumption: HR = 0.87, 95% C.I. = 0.77–0.98 for low consumption; HR = 0.86, 95% C.I. = 0.77–0.96 for moderate consumption; HR = 0.84, 95% C.I. = 0.75–0.95 for high consumption). When stratified by sex and smoking status, this reduced risk was statistically significant among men and current and former smokers. In addition, dose–response analyses showed a lower bladder cancer risk with increment of 100 ml of tea consumption per day (HR-increment = 0.97; 95% CI = 0.96–0.98). A similar inverse association was found among males, current and former smokers while never smokers and females showed non-significant results, suggesting potential sex-dependent effect. Higher consumption of tea is associated with reduced risk of bladder cancer with potential interaction with sex and smoking status. Further studies are needed to clarify the mechanisms for a protective effect of tea (e.g. inhibition of the survival and proliferation of cancer cells and anti-inflammatory mechanisms) and its interaction with smoking and sex.

60 APPLIED LIFE SCIENCES↗