Search NASA⌕ Search

SEARCH · Search NASA

Results for “early formulation”

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 55 records · Page 3

ASSESSMENT OF HIGH-VOLUME HARVESTED FLY ASH BLENDS FOR USE IN PRECAST CONSTRUCTION

Declining supplies of fresh fly ashes in some regions have further driven the importance of using harvested fly ash as a supplementary cementitious material and thus replacement for more energy-intensive ordinary Portland cement (OPC). Higher replacement rates of OPC with fly ash in general has traditionally been limited in precast operations due to stringent early-age strength requirements necessitated by initial prestress and lifting/handling of hardened components – often well within 24 hours after fresh concrete placement. Therefore, a series of mix designs were developed, as part of a larger framework for reassessing high-volume fly ash (HVFA) use in the context of precast construction, to demonstrate the feasibility of using larger fractions of harvested fly ash in concrete formulations suitable for precast fabrication demands. More specifically, a target minimum 24-hour compressive strength of 24.1 MPa (3500 psi) was set, in accordance with conventional precast structural design checks, while increasing the fly ash content to 40% replacement of OPC – at least a 15% increase above the traditional limit of 25%. The effect of using more sustainable Type IL cement in contrast to conventional Type III Portland cement (the traditional standard for precast use) – both in conjunction with HVFA use – was also examined. Early-age mechanical properties including compressive strength, modulus of rupture, and modulus of elasticity were evaluated within the aforementioned early-age window. Lastly, corresponding design equations were then reassessed for their applicability with respect to the high-early strength HVFA mixes examined herein.

Gombeda, Matthew J.↗

Aluminosilicate colloidal gels: from the early age to the precipitation of zeolites

Aluminosilicate hydrogels are often considered to be precursors for the crystallisation of zeolites carried out under hydrothermal conditions. The preparation of mechanically homogeneous aluminosilicate gels enables the study of these materials through bulk rheology and observation of the aging dynamics until the precipitation of crystalline zeolites. The first part of this study deals with the establishment of ternary state diagrams, in order to identify the range of chemical formulations that enable preparation of single-phase homogeneous gels. Then, by studying the viscoelastic moduli during the gelation reaction, and by yielding the gel under large deformation, we propose an empirical law considering the partial order of reaction on each chemical element, to predict the gelation time according to the chemical formulation. The scaling behavior of the elastic properties of this colloidal gel shows a transition from a strong link behavior to a weak link regime. Long term aging results in the shrinkage of the gel, accompanied by syneresis of interstitial liquid at the surface. Zeolites precipitate through crystallisation by a particle attachment mechanism, when thermodynamic equilibrium is reached. The stoichiometry of the precipitated zeolites is not only consistent with the concentration of the remaining species in the supernatant but, surprisingly, it is also very close to the partial order of the reaction of the chemical elements involved in the determination of the critical gel point. This indicates a strong correlation between the morphology of the soft amorphous gel network that is formed at an early age and those of the final solid precipitated crystals.

36 MATERIALS SCIENCE↗

Rapid early gas accretion for the inner Galactic disc: A case for short accretion timescale

Recent observations of the Milky Way and galaxies at high redshifts suggest that galaxy discs were already in place soon after the Big Bang. While the gas infall history of the Milky Way in the inner disc has long been assumed to be characterised by a short accretion timescale, this has not been directly constrained using observations. Using data for the inner regions of the Milky Way recently produced by APOGEE and Gaia and of unprecedented quantity and quality, we aim to derive strong constraints on the infall history of the inner (< 6 kpc) Galaxy (with a focus on stars between 4 and 6 kpc, which we show is an appropriate proxy for the entire inner disc). We implemented gas infall into a chemical evolution model of the Galaxy disc, and used a Schmidt–Kennicutt law to connect the infall to the star formation. We explore a number of models, and two different formulations of the infall law. In one formulation, the infall is non-parametric, and in the other the infall has an explicitly exponential form. We fit the model parameters to the time–[Si/Fe] distribution of solar vicinity stars, and the metallicity and [Si/Fe] distribution function of stars with a galactocentric radius of between 4 and 6 kpc from APOGEE. Our results point to a fast, early gas accretion, and an upper limit on the accretion timescale of around 2 Gyr in the inner disc of the Milky Way. This suggests that at least half the baryons were in place within 2–3 Gyr of the Big Bang, and that half the stars of the inner disc formed within the first 5 Gyr, during the thick disc formation phase. This implies that the stellar mass of the inner disc is dominated by the thick disc, supporting our previous work, and that the gas accretion onto the inner disc was rapid and early.

79 ASTRONOMY AND ASTROPHYSICS↗

Non-linear dynamics of jet quenching

We develop a comprehensive analytic framework for jet quenching in QCD media, based on a medium-induced parton cascade sourced by collinear virtual splittings. We show that the energy flow out of the jet cone, driven by turbulent gluon cascades, is governed by a non-linear rate equation that resums gluon splittings at arbitrary angles and is enhanced by the medium length, L. The solution of this equation sets the initial condition for a non-linear DGLAP-like evolution equation, which describes the collinear early vacuum cascade resolved by the medium at angles exceeding the medium resolution angle, θ c . For asymptotic jet energies, the medium-induced cascade displays an exponential behavior that generalizes the Poisson-like distribution of parton energy loss. This formulation enables the resummation of leading contributions in α s ln(1/R), and α s ln(R/θ c ), and powers of α s L. We briefly explore the limit of strong quenching, where analytic treatments are feasible, offering insights into the impact of parton cascades on jet quenching. These results provide guidance for future numerical simulations and analytical investigations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

United States Early Career Researchers in Collider Physics input to the European Strategy for Particle Physics Update

This document represents a contribution of the United States early career collider physics community to the 2025--2026 update to the European Strategy for Particle Physics. Preferences with regard to different future collider options and R&D priorities were assessed via a survey. The early career community was defined as anyone who is a graduate student, postdoctoral researcher, untenured faculty member, or research scientist under 40 years of age. In total, 105 participants responded to the survey between February and March 10th, 2025. Questions were formulated primarily to gauge the enthusiasm and preferences for different collider options in line with the recommendations of the United States' P5 report, relevant to the European Strategy Update.

43 PARTICLE ACCELERATORS↗

The influence of electrode crack dimensions on the durability of polymer electrolyte membrane fuel cells

Electrode cracks in polymer electrolyte membrane fuel cells (PEMFCs) are correlated with early onset failures. Here, in this work we investigate the influence of cracked gas diffusion electrodes (GDEs) on the durability of the membrane electrode assembly (MEA) using a combined chemical-mechanical accelerated stress test (AST). Electrode crack dimensions were systematically tuned using ink formulations and material selection strategies. A parameter to describe the crack width areal density (Φ CW ) was used to quantify the degree of discontinuity in the electrode surfaces. Open circuit voltage (OCV) transient analyses were used to benchmark and characterize the failure mechanisms in the MEAs as a function of the Φ CW . While smaller electrode-level cracks, on the order of microns, yielded a 28 % decrease in operating lifetime, larger cracks that propagated from a discontinuous, microporous layer (MPL) coating, decreased the operating lifetime by 56 %. This work emphasizes the need for material processing strategies that consider defect tolerances to limit membrane failures in PEMFCs.

08 HYDROGEN↗

Emerging Potassium Metal Anodes: Perspectives on Control of the Electrochemical Interfaces

Potassium metal acts as the anode in emerging potassium metal batteries (KMBs). It further serves as the counter-electrode for potassium ion battery (KIB) half-cells, with its reliable performance being critical for assessing the working electrode material. This first-of-its-kind critical review focuses on the dual challenge of controlling the potassium metal-substrate and the potassium metal-electrolyte interface so as to prevent dendrites. The discussion begins with a comparison of the physical and chemical properties of K metal anodes versus the much oft studied Li and Na metal anodes. Based on established descriptions for root causes of dendrites, the problem should be less severe for K than for Li or Na, while in fact the opposite is observed. The key reason that the K metal surface rapidly becomes dendritic in common electrolytes is its unstable solid electrolyte interphase (SEI). An unstable SEI layer is defined as being non self-passivating. No SEI is perfectly stable during cycling, and all SEI structures are heterogenous both vertically and horizontally relative to the electrolyte interface. The difference between a "stable" and an "unstable" SEI may be viewed as the relative degree to which during cycling it thickens and becomes further heterogeneous. The unstable SEI on K metal leads to a number of interrelated problems, such as low cycling Coulombic efficiency (CE), a severe impedance rise, large overpotentials, and possibly electrical shorting; all of which have been reported to occur as early as in the first ten plating/stripping cycles. Many of the traditional "interface fixes" employed for Li and Na metal anodes, such as various artificial SEIs, surface membranes, barrier layers, the secondary separators, etc. have not been attempted or optimized for the case of K. This is an important area for further exploration, with an understanding that success may come harder than with Li due to K-based SEI reactivity with both ether and ester solvents. The second critical problem with K metal anodes is that they do not thermally or electrochemically wet a standard (untreated) Cu foil current collector. Published experimental and modeling research directly highlights the weak bonding between the K atoms and a Cu surface. Existing surface treatment approaches that achieve improved K wetting are discussed, along with the general design rules for future studies. Also discussed are geometry-based methods to tune nucleation, as well dual approaches where nucleation and SEI structure are tuned through complementary schemes to achieve extended half-cell and full battery stability. We hypothesize that K metal never achieves a planar wetting morphology even at cycle one making the dendrites "baked-in". We propose that classical thin films growth models, Frank van der Merwe (F-M), Volmer Weber (V-W), and Stranski-Krastanov (S-K) can be employed to describe early stage plating behavior. It is demonstrated that island-like V-W growth is the applicable description for the natural plating behavior of K on pristine Cu. Moving forward, there are three inter-related thrusts to be pursued: First, K salt - based electrolyte formulations have to mature and become further tailored to handle the increased reactivity of a metal rather than an ion anode. Second, the K-based SEI structure needs to be further understood and ultimately tuned to be less reactive. Third, the energetics of the K metal - current collector interface must be controlled to promote planar wetting/dewetting throughout cycling.

25 ENERGY STORAGE↗

Sequential formulation of all‐way coupled finite strain thermoporomechanics for largely deformable gas hydrate deposits

We develop a numerically stable sequential formulation of thermoporomechanics for largely deformable gas hydrate deposits, extended from the fixed stress split of infinitesimal transformation. Constitutive equations are based on the total Lagrangian approach for both flow and geomechanics, including dynamic full tensor permeability and thermal conductivity updated from the deformation gradient. For space discretization, we take the cell-centered finite volume and node-based finite element method for flow and geomechanics, respectively. Then, we propose a sequential implicit method for all-way coupled thermoporomechanics, where the nonisothermal multiphase flow problem of gas hydrates is solved implicitly first and then the geomechanics problem is solved implicitly at the next step. During solution of the flow problem, we fix the rate of first Pioal total stress for numerical stability as well as apply porosity correction and entropy correction to account for geomechanical effects. We test numerical examples where flow and geomechanics parameters are based on deep oceanic gas hydrate deposits. When applying depressurization, even though the results between the infinitesimal transformation and finite strain geomechanics are similar in the early stages due to small deformation, we find differences between them in the late times as deformation becomes large. Accordingly, permeability and thermal conductivity tensors become nonisotropic full tensors although they are initially isotropic. Furthermore, we identify numerical stability of the developed sequential method from the test cases that exhibit the highly complex coupled gas hydrate systems with large deformation. Thus, the proposed sequential formulation can be applied in largely deformable gas hydrate systems.

42 ENGINEERING↗

Developing Technology Performance Level Assessments for Early-Stage Wave Energy Converter Technologies: Preprint

The advantage of using Technology Performance Level (TPL) in conjunction with Technology Readiness Level (TRL) assessments in guiding technology development trajectories to successful outcomes in less time, at less overall cost, and with less encountered risk has been well articulated in the literature. In partnership with industry and international collaborators, a TPL assessment methodology for grid-connected applications has been developed through the application of the systems engineering approach. Metrics under seven different categories have been developed, weighted based on their relative relevance, and combined to yield a composite score. The methodology has been implemented in a spreadsheet tool plus a web application specifically aimed at assessing early stage (TRL 1-3) concepts. The target use cases are (a) technology developers improving their design, to find fatal flaws early, to get feedback on current design, to identify areas of improvement that will yield the highest return on investment, (b) reviewers assessing technologies in competitions or for making funding decisions, (c) investor or project developer doing due diligence, (d) policy makers landscaping the technology domain for formulating R&D strategy. The methodology and the tools are undergoing continuous improvement based on the experience and lessons learnt from applying it to internal and external marine energy technology development projects. The methodology is also being adapted for assessing WECs servicing markets outside the continental grid - broadly categorized as Powering the Blue Economy (PBE) applications. Such applications have vastly different functional requirements entailing a modification of the methodology to account for their higher risk tolerance, reduced price sensitivities, lower power needs, different permitting protocols, etc. This paper presents the latest status of the TPL assessment methodology and tools, describes its adaptation to select PBE markets, and explores its extension to other domains where it could provide a comprehensive and holistic measure of a nascent or disruptive technology's technoeconomic performance potential.

metrics↗

Developing Technology Performance Level Assessments for Early-Stage Wave Energy Converter Technologies

The advantage of using Technology Performance Level (TPL) in conjunction with Technology Readiness Level (TRL) assessments in guiding technology development trajectories to successful outcomes in less time, at less overall cost, and with less encountered risk has been well articulated in the literature. In partnership with industry and international collaborators, a TPL assessment methodology for grid-connected applications has been developed through the application of the systems engineering approach. Metrics under seven different categories have been developed, weighted based on their relative relevance, and combined to yield a composite score. The methodology has been implemented in a spreadsheet tool plus a web application specifically aimed at assessing early stage (TRL 1-3) concepts. The target use cases are (a) technology developers improving their design, to find fatal flaws early, to get feedback on current design, to identify areas of improvement that will yield the highest return on investment, (b) reviewers assessing technologies in competitions or for making funding decisions, (c) investor or project developer doing due diligence, (d) policy makers landscaping the technology domain for formulating R&D strategy. The methodology and the tools are undergoing continuous improvement based on the experience and lessons learnt from applying it to internal and external marine energy technology development projects. The methodology is also being adapted for assessing WECs servicing markets outside the continental grid - broadly categorized as Powering the Blue Economy (PBE) applications. Such applications have vastly different functional requirements entailing a modification of the methodology to account for their higher risk tolerance, reduced price sensitivities, lower power needs, different permitting protocols, etc. This paper presents the latest status of the TPL assessment methodology and tools, describes its adaptation to select PBE markets, and explores its extension to other domains where it could provide a comprehensive and holistic measure of a nascent or disruptive technology's technoeconomic performance potential.

metrics↗

Nearly optimal state preparation for quantum simulations of lattice gauge theories

Here, we present several improvements to the recently developed ground-state preparation algorithm based on the quantum eigenvalue transformation for unitary matrices (QETU), apply this algorithm to a lattice formulation of U(1) gauge theory in (2+1) dimensions, as well as propose an alternative application of QETU, a highly efficient preparation of Gaussian distributions. The QETU technique was originally proposed as an algorithm for nearly optimal ground-state preparation and ground-state energy estimation on early fault-tolerant devices. It uses the time-evolution input model, which can potentially overcome the large overall prefactor in the asymptotic gate cost arising in similar algorithms based on the Hamiltonian input model. We present modifications to the original QETU algorithm that significantly reduce the cost for the cases of both exact and Trotterized implementation of the time evolution circuit. We use QETU to prepare the ground state of a U(1) lattice gauge theory in two spatial dimensions, explore the dependence of computational resources on the desired precision and system parameters, and discuss the applicability of our results to general lattice gauge theories. We also demonstrate how the QETU technique can be utilized for preparing Gaussian distributions and wave packets in a way which outperforms existing algorithms for as little as n q ≳ 2–5 qubits.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Microemulsions: Breakthrough Electrolytes for Redox Flow Batteries

Aqueous and non-aqueous redox flow batteries (RFBs) have limited energy and current densities, respectively, due to the nature of the electrolytes. New approaches to electrolyte design are needed to improve the performance of RFBs. In this work, we combined a highly conductive aqueous phase and an organic redox-active phase in a microemulsion to formulate a novel RFB electrolyte. As a proof-of-concept, we demonstrate an RFB using this microemulsion electrolyte with maximum current density of 17.5 mA·cm −2 with a 0.19 M posolyte and 0.09 M negolyte at a flow rate of only ∼2.5 ml·min −1 , comparable to early vanadium electrolyte RFBs at similar flow rates on a per molar basis. The novel active negolyte component is an inexpensive oil-soluble vitamin (K 3 ). By combining aqueous and organic phases, the solvent potential window and energy density may be increased without sacrificing current density and new redox couples may be accessed. Microemulsion electrolytes show great promise for improved performance and increased energy densities in aqueous RFBs but the path forward is complex. We end with discussion of areas that need work to achieve the potential of these electrolytes.

25 ENERGY STORAGE↗

Mapping Incidence and Prevalence Peak Data for SIR Modeling Applications

Infectious disease modeling and forecasting have played a key role in helping assess and respond to epidemics and pandemics. Recent work has leveraged data on disease peak infection and peak hospital incidence to fit compartmental models for the purpose of forecasting and describing the dynamics of a disease outbreak. Incorporating these data can greatly stabilize a compartmental model fit on early observations, where slight perturbations in the data may lead to model fits that forecast wildly unrealistic peak infection. We introduce a new method for incorporating historic data on the value and time of peak incidence of hospitalization into the fit for a Susceptible-Infectious-Recovered (SIR) model by formulating the relationship between an SIR model’s starting parameters and peak incidence as a system of two equations that can be solved computationally. We demonstrate how to calculate SIR parameter estimates – which describe disease dynamics such as transmission and recovery rates – using this method, and determine that there is a noticeable loss in accuracy whenever prevalence data is misspecified as incidence data. To exhibit the modeling potential, we update the Dirichlet-Beta State Space modeling framework to use hospital incidence data, as this framework was previously formulated to incorporate only data on total infections. This approach is assessed for practicality in terms of accuracy and speed of computation via simulation.

97 MATHEMATICS AND COMPUTING↗

Assessment of an E10 gasoline surrogate: Qualitative and quantitative comparisons of in-cylinder spray morphology

A minimum-component gasoline fuel surrogate that captures both chemical and physical behaviors of a full-distillate fuel is needed for high-fidelity CFD simulations. This study evaluates gasoline spray characteristics in a direct-injection spark-ignition engine under motored operation. Two fuels are compared; PACE-20, which is a 9-component surrogate formulation of RD5-87, is compared with its target fuel RD5-87, which is a full-boiling range research grade E10 gasoline. The spray morphologies of both fuels are recorded for a centrally-located direct-injection 8-hole spray subject to intake air cross-flow during the early part of the intake stroke. High-speed imaging recorded scattered light of the side and axial projections of the liquid spray. Quantitative metrics were developed and employed to facilitate comparison of spray morphologies as well as to identify the transition in spray morphology due to flash boiling. This paper builds on a previous study of RD5-87 where coolant temperature (20°C–100°C), in-cylinder pressure (40–110 kPa), engine speed (650–1950 rpm), and injection pressure (60–180 bar) were systematically changed to span operating conditions with and without flash boiling. Images of the PACE-20 morphology are selected for a sub-set of operating conditions from the previous study where distinctive morphology changes occurred. Visual inspection of the images and quantitative metrics demonstrate that the PACE-20 spray morphology is equivalent to that of the RD5-87 in most cases. The exception was for changes in the ambient-gas pressure where the flash-boiling transition occurred at ∼5 kPa higher in-cylinder pressure for PACE-20. Three empirical metrics, Merging Index, Asymmetry, and Flash Index are proposed here and they were found to be useful both as quantitative comparisons of the fuel morphologies, and for identifying the transition in spray morphology due to flash boiling.

Kim, Namho↗

Seismic Tomography 2024

Seismic tomography is the most abundant source of information about the internal structure of the Earth at scales ranging from a few meters to thousands of kilometers. It constrains the properties of active volcanoes, earthquake fault zones, deep reservoirs and storage sites, glaciers and ice sheets, or the entire globe. It contributes to outstanding societal problems related to natural hazards, resource exploration, underground storage, and many more. The recent advances in seismic tomography are being translated to nondestructive testing, medical ultrasound, and helioseismology. Nearly 50 yr after its first successful applications, this article offers a snapshot of modern seismic tomography. Focused on major challenges and particularly promising research directions, it is intended to guide both Earth science professionals and early-career scientists. The individual contributions by the coauthors provide diverse perspectives on topics that may at first seem disconnected but are closely tied together by a few coherent threads: multiparameter inversion for properties related to dynamic processes, data quality, and geographic coverage, uncertainty quantification that is useful for geologic interpretation, new formulations of tomographic inverse problems that address concrete geologic questions more directly, and the presentation and quantitative comparison of tomographic models. Finally, it remains to be seen which of these problems will be considered solved, solved to some extent, or practically unsolvable over the next decade.

58 GEOSCIENCES↗

Metal Independent Correlations for Site-specific Binding Energies of Relevant Catalytic Intermediates

Establishing energy correlations among different metals can accelerate the discovery of efficient and cost-effective catalysts for complex reactions. Using a recently introduced coordination-based model, we can predict site specific metal binding energies (∆EM) that can be used as a descriptor for chemical reactions. In this study, we have examined a range of metals including Ag, Au, Co, Cu, Ir, Ni, Os, Pd, Pt, Rh, and Ru and found linear correlations between predicted ∆EM and adsorption energies of CH and O (∆ECH and ∆EO) at various coordination environment for all the considered metals. Interestingly, all the metals correlate with one another under specific surface site coordination, indicating that different metals are interrelated in a particular coordination environment. Furthermore, we have tested and verified for PtPd and PtIr-based alloys that they follow a similar behavior. Moreover, we have expanded the metal space by taking some early transition metals along with few s-block metals and shown a cyclic behavior of the adsorbate binding energy (∆EA) versus ∆EM. Therefore, ∆ECH and ∆EO can be efficiently interpolated between metals, alloys and intermetallics based on information related to one metal only. This simplifies the process of screening new metal catalyst formulations and their reaction energies.

Density function theory↗

Evaluation of Spray and Combustion Models for Simulating Dilute Combustion in a Direct-Injection Spark-Ignition Engine

Dilute combustion in spark-ignition engines has the potential to improve thermal efficiency by mitigating knock and by reducing throttling and wall heat losses. However, ignition and combustion processes can become unstable for dilute operation due to a lowered laminar flame speed, resulting in excessive cycle-to-cycle variability (CCV) of the combustion process. To compensate for the slower combustion in less reactive mixtures, a modified intake port geometry can be employed to generate a strong tumble flow in the cylinder and elevate turbulence levels around the spark plug, thereby promoting a faster transition to turbulent deflagration. Consequently, optimizing combustion chamber geometry and operating strategy is crucial to maximizing the benefits of using dilute combustion with enhanced in-cylinder turbulence across a wide range of operating conditions. Computational fluid dynamics (CFD) simulations can be utilized for virtual engine optimization tasks, but this would require the models to be truly predictive regarding the impact of changes to the engine design and operational parameters.In this study, multicycle large-eddy simulations (LES) are performed for a direct-injection spark-ignition engine to investigate the model performance in predicting engine combustion characteristics with respect to changes in the intake configuration. A tumble plate that blocks the lower part of the intake port inlet is used to vary the tumble. A set of CFD models that have been recently developed are employed, which takes into account the drag of nonspherical droplets, flash-boiling behavior of liquid sprays, spray-wall interaction, surrogate formulation of a research-grade E10 gasoline, and fast chemical kinetic solvers. Simulation results are compared to experimental engine data in terms of cylinder pressure, apparent heat release rate, mass fraction burned timing, and flame images. It is found that LES employing the state-of-the-art CFD models are capable of properly predicting the spray processes and reproducing the measured mean cylinder pressure for the case with the tumble plate. On the other hand, the LES over-predicts the combustion rate during the early combustion stage and under-estimates the CCV, and these discrepancies become larger when the tumble plate is removed.

computational fluid dynamics simulation↗

The emergence and transmission dynamics of HIV-1 CRF07_BC in Mainland China

A total of 1155 partial pol gene sequences of human immunodeficiency virus (HIV)-1 CRF07_BC were sampled between 1997 and 2015, spanning 13 provinces in Mainland China and risk groups [heterosexual, injecting drug users (IDU), and men who have sex with men (MSM)] to investigate the evolution, adaptation, spatiotemporal and risk group dynamics, migration patterns, and protein structure of HIV-1 CRF07_BC. Due to the unequal distribution of sequences across time, location, and risk group in the complete dataset (‘full1155’), subsampling methods were used. Maximum-likelihood and Bayesian phylogenetic analysis as well as discrete trait analysis of geographical location and risk group were carried out. To study mutations of a cluster of HIV-1 CRF07_BC (CRF07-1), we performed a comparative analysis of this cluster to the other CRF07_BC sequences (‘backbone_295’) and mapped the mutations observed in the respective protein structure. Our findings showed that HIV-1 CRF07_BC most likely originated among IDU in Yunnan Province between October 1992 to July 1993 [95 per cent hightest posterior density (HPD): May 1989–August 1995] and that IDU in Yunnan Province and MSM in Guangdong Province likely served as the viral sources during the early and more recent spread in Mainland China. We also revealed that HIV-1 CRF07-1 has been spreading for roughly 20 years and continues to cause local transmission in Mainland China and worldwide. Overall, our study sheds light on the dynamics of HIV-1 CRF07_BC distribution patterns in Mainland China. Our research may also be useful in formulating public health policies aimed at controlling acquired immune deficiency syndrome in Mainland China and globally.

60 APPLIED LIFE SCIENCES↗