SEARCH · Engineering Papers

Results for “Transport phenomena”

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.

SourceResearch connectionFrom yearThrough year

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

At least 649 records · Page 36

All sources NASA NTRS DOE OSTI

DOE OSTI2022

(Invited) Continuum Mathematical Modeling of Water Electrolysis: A Tutorial

Widespread use of hydrogen energy is contingent on the development of reliable and economical sources of hydrogen. Electrolysis from renewably-derived low-carbon electricity is a potentially viable method of hydrogen generation. Prime among the electrolysis technologies are those utilizing ion-conducting polymers (ionomers) including proton-exchange-membrane water electrolyzer (PEMWE). However, these technologies need to exhibit increased efficiency, performance, and durability to become commercially viable. Like most electrochemical devices, PEMWEs involve multiple components (e.g., catalyst, ionomer, transport layers, membrane, plates) and multiple phases, with phenomena occurring across different time and length scales. Furthermore, it is difficult to experimentally probe many of the species and phenomena during operation. Thus, mathematical modeling at the continuum level has been an invaluable aid in exploring, understanding, and optimizing PEMWE cell and components. Furthermore, this is especially true in the highly coupled and complex physics and chemistries that occur with the membrane-electrode assembly (MEA). The physics in a typical volume-averaged non-isothermal model include multiphase transport in porous media, concentrated-solution theory, Ohm's Law, and Butler-Volmer kinetics, and ion, gas, and water transport in the ionomer.

Dizon, Arthur↗

DOE OSTI2025

Bridging interfacial properties and cell performance: A multiscale model for proton-exchange-membrane fuel cells

Here, to elucidate the impact of local interfaces on mass-transport resistance and overall cell performance of low-loaded proton-exchange-membrane fuel cells (PEMFCs), we present a multiscale modeling framework incorporating a novel modified agglomerate model. The model considers three distinct Pt-electrolyte interfaces: Pt on the carbon surface covered by either ionomer or water film and Pt inside carbon nanopores. Detailed mass-transport voltage-loss breakdowns reveal that coupled agglomerate-interface-scale mass transport dominates the mass-transport loss. The ionomer poisons the exterior-Pt surface through suppressing O 2 adsorption and intrinsic ORR activity, leading to low current-density performance. Conversely, interior-Pt interface enhances the kinetic performance but limits high current-density performance due to its low interfacial permeability. The exterior-Pt/water interface demonstrates superior kinetic performance and mass transport, though its practical implementation requires ensuring proton transport. By coupling the multiscale CL properties with ink parameters, the model identifies an optimal I to C ratio of approximately 0.5, a moderate value where the ionomer content is sufficient to guarantee proton transport without fully covering the Pt surface and forming large agglomeration, thus allowing the utilization of the Pt-water interface and avoiding high mass-transport loss. Overall, the model helps unravel limiting phenomena across different operating regimes and provides routes for optimizing performance.

Cell diagnostic↗

DOE OSTI2022

Quantum-Accurate Multiscale Modeling of Shock Hugoniots, Ramp Compression Paths, Structural and Magnetic Phase Transitions, and Transport Properties in Highly Compressed Metals

Fully characterizing high energy density (HED) phenomena using pulsed power facilities (Z machine) and coherent light sources is possible only with complementary numerical modeling for design, diagnostic development, and data interpretation. The exercise of creating numerical tests, that match experimental conditions, builds critical insight that is crucial for the development of a strong fundamental understanding of the physics behind HED phenomena and for the design of next generation pulsed power facilities. The persistence of electron correlation in HED materials arising from Coulomb interactions and the Pauli exclusion principle is one of the greatest challenges for accurate numerical modeling and has hitherto impeded our ability to model HED phenomena across multiple length and time scales at sufficient accuracy. An exemplar is a ferromagnetic material like iron, while familiar and widely used, we lack a simulation capability to characterize the interplay of structure and magnetic effects that govern material strength, kinetics of phase transitions and other transport properties. Herein we construct and demonstrate the Molecular-Spin Dynamics (MSD) simulation capability for iron from ambient to earth core conditions, all software advances are open source and presently available for broad usage. These methods are multi-scale in nature, direct comparisons between high fidelity density functional theory (DFT) and linear-scaling MSD simulations is done throughout this work, with advancements made to MSD allowing for electronic structure changes being reflected in classical dynamics. Main takeaways for the project include insight into the role of magnetic spins on mechanical properties and thermal conductivity, development of accurate interatomic potentials paired with spin Hamiltonians, and characterization of the high pressure melt boundary that is of critical importance to planetary modeling efforts.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1981

A simplified model of high pressure spray combustion

A simplified model of high-pressure spray combustion is examined. The analysis relies on a kappa-epsilon-g turbulence model in conjunction with the locally homogeneous flow (LHF) approximation of two-phase flow, which implies infinitely fast transport rates between the phases. High-pressure phenomena near the thermodynamic critical point are treated using the Redlich-Kwong equation of state. Predictions are compared with existing measurements of spray boundaries in a pressure-atomized n-pentane spray (Sauter mean diameter, approximately 30 microns) burning in stagnant air at 3, 6, and 9 MPa. The LHF model overestimates the rate of development of the flow, yielding spray lengths roughly 20% shorter than measured. Calibrated drop-life-history calculations suggest that finite interphase transport rates are the primary cause of the discrepancy.

Mao, C.-P.↗

DOE OSTI2025

Field Induced Density Wave in a Kagome Superconductor

On the kagome lattice, electrons benefit from the simultaneous presence of band topology, flat electronic bands, and Van Hove singularities, forming competing or cooperating orders. Understanding the interrelation between these distinct order parameters remains a significant challenge, leaving much of the associated physics unexplored. In the kagome superconductor KV 3 Sb 5 , which exhibits a charge density wave (CDW) state below T ≃ 78 K, we uncover an unpredicted field-induced phase transition below 6 K. The observed transition is marked by a hysteretic anomaly in the resistivity, nonlinear electrical transport, and a change in the symmetry of the electronic response as probed via the angular dependence of the magnetoresistivity. These observations surprisingly suggest the emergence of an unanticipated broken symmetry state coexisting with the original CDW. To understand this experimental observation, we developed a theoretical minimal model for the normal state inside the high-temperature parent CDW phase, where an incommensurate CDW order emerges as an instability subleading to superconductivity. Here, the incommensurate CDW emerges when superconducting fluctuations become fully suppressed by large magnetic fields. Our results suggest that, in kagome superconductors, quantum states can either coexist or are nearly degenerate in energy, indicating that these are rich platforms to expose new correlated phenomena.

Charge density waves↗

DOE OSTI2021

Enhanced Diffusive Transport in Fluctuating Porous Media

Mass transport within porous structures is a ubiquitous process in biological, geological, and technological systems. Despite the importance of these phenomena, there is no comprehensive theory that describes the complex and diverse transport behavior within porous environments. While the porous matrix itself is generally considered a static and passive participant, many porous environments are in fact dynamic, with fluctuating walls, pores that open and close, and dynamically changing cross-links. While diffusion has been measured in fluctuating structures, notably in model biological systems, it is rarely possible to isolate the effect of fluctuations because of the absence of control experiments involving an identical static counterpart, and it is generally impossible to observe the dynamics of the structure. In this paper, we present a direct comparison of the diffusion of nanoparticles of various sizes within a trackable, fluctuating porous matrix and a geometrically equivalent static matrix, in conditions spanning a range of regimes from obstructed to highly confined. The experimental system comprised a close-packed layer of colloidal spheres that were either immobilized to a planar surface or allowed to fluctuate locally, within the space defined by their nearest neighbors. Interestingly, the effective long-time diffusion coefficient was approximately 35–65% greater in the fluctuating porous matrix than in the static one (depending on the size of the nanoparticle probes), regardless of the geometric regime. This was explained by considering the enhancing effects of matrix fluctuations on the short-time diffusion coefficient and cooperative “gate-opening” motions of matrix particles and nanoparticle probes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Ultrafast low-temperature metal–insulator interface phonon dynamics and heat transport in a Pt/Gd 3 Fe 5 O 12 heterostructure

Interfacial thermal and acoustic phenomena have an important role in quantum science and technology, including in spintronic and spincaloritronic materials and devices. Simultaneous measurements of the low-temperature thermal and acoustic properties of a metal/insulator heterostructure reveal distinct dynamics in the characteristic phonon frequency ranges of acoustic and thermal transport. The measurements probed a heterostructure consisting of a thin film of Pt on the ferrimagnetic insulator gadolinium iron garnet (Gd 3 Fe 5 O 12 , GdIG) grown epitaxially on a gadolinium gallium garnet substrate. Ultrafast structural dynamics within the Pt layer were tracked using time-resolved ultrafast x-ray diffraction and analyzed to probe interfacial acoustic and thermal properties. The rapid heating of the Pt layer by a 400 nm wavelength femtosecond-duration optical pulse produced transient structural changes that provided the stimulus for these measurements. Rapid heating produced a broadband acoustic pulse that was partially reflected by the Pt/GdIG interface. Temporal frequencies up to 740 GHz, corresponding to angular frequencies of several THz, were detected in a wavelet analysis of the acoustic oscillations of the strain in the Pt layer. The structural results were analyzed to determine (i) the acoustic damping coefficient and phonon mean free path in Pt at frequencies of hundreds of GHz and (ii) the Grüneisen anharmonicity parameter. The thermal conductance of the Pt/GdIG interface was tracked using the slower, tens-of-picosecond-scale, dynamics of the initial cooling of the heated Pt layer. Analysis using a model based on the Boltzmann transport equation shows that the phonon transmission is lower at the phonon frequencies relevant to thermal transport than for subterahertz regime acoustics.

Acoustic phenomena↗

NASA NTRS2000

Exploiting the Temperature Dependence of Magnetic Susceptibility to Control Convection in Fundamental Studies of Solidification Phenomena

It is well known that convection is a dominant mass transport mechanism when materials are solidified on Earth's surface. This convection is caused by gradients in density (and therefore gravitational force) that are brought about by gradients in temperature, composition or both. Diffusion of solute is therefore dwarfed by convection and the study of fundamental parameters, such as dendrite tip shape and growth velocity in the absence of convection is nearly impossible. Significant experimental work has therefore been carried out in orbiting laboratories with the intent of minimizing convection by minimizing gravity. One of the best known experiments of this kind is the Isothermal Dendritic Growth Experiment (IDGE), supported by NASA. Naturally such experiments are costly and one objective of the present investigation is to develop an experimental method whereby convection can be- halted, in solidification and other experiments, on the surface. A second objective is to use the method to minimize convection resulting from the residual accelerations suffered by experiments in microgravity. The method to be used to minimize convection relies on the dependence of the magnetic susceptibility of a fluid on temperature or composition (whichever is driving convection). All materials experience a force when placed in a magnetic field gradient. The direction and magnitude of that force depend on the magnetic susceptibility of the material. Consequently the force will vary if the susceptibility varies with temperature or composition. With a magnetic field gradient in the right direction (typically upward) and of the right magnitude, this variation in the magnetic force can be made to exactly cancel the variation in the gravitational force. Expressed another way, normal buoyancy is exactly countered by a "magnetic buoyancy". To demonstrate the principle, a solution of MnC12 in water has been used. First the variation of the susceptibility of this paramagnetic solution with temperature and concentration was measured. Then a "cell", containing this solution and 50mm long by 15mm high by 155mm wide, was placed in a superconducting magnet at Marshall Space Flight Center. The magnetic field was measured at various positions within the bore of the magnet using a Hall effect probe. In this way, a position was found where the magnetic field gradient was predominantly upward; the magnitude of the gradient could then be adjusted by adjusting the current of the magnet. The ends of the cell consisted of machined copper blocks maintained at controlled temperatures by circulating water from constant temperature baths. The walls of the cell were of rectangular section glass tubing so that the cell contents could be seen. Velocities arising from thermal gradients within the cell were measured by particle image velocimetry (PIV). Particles used for this purpose were silver-coated hollow glass spheres of micrometers diameter and nearly the same density as the solution. A central vertical plane of the cell was illuminated by a laser beam passing through a cylindrical lens. Digital images of the particles were captured on a CCD camera and fed to a computer so that frame-to-frame movements of particles traveling with the fluid were captured. These images were employed to compute velocity maps using commercial PIV software. In a typical experiment the cold end of the cell was maintained at 10C and the warm end at 30 C. With no current in the magnet, i.e.- with natural -convection allowed to occur, the fluid was observed to circulate with an average speed of approximately 0.3 millimeters per second. It was visually apparent that this circulation was diminished as the current was increased. At currents of approximately 20A the flow was halted, to within the precision of the PIV measurements. At yet higher currents the convection was reversed with the hotter solution sinking and the cooler solution rising. At 40A this reversed convection had speeds averaging 0.43 millimeters per second. The measurements of susceptibility and density allow an estimate of the field gradient necessary to halt convection in the experiment. That estimate was 7.8T (squared) per meter and the convection was observed to halt in the magnet at a current giving 7.21T (squared) per meter from the magnetic field measurements. Calculations of the flow have been carried out using the computational fluid dynamics software FLUENT and show good agreement with the measurements.

Seybert, C. D.↗

NASA NTRS2001

Exploiting the Temperature Dependence of Magnetic Susceptibility to Control Convective in Fundamental Studies of Solidification Phenomena

Seybert, C.↗

DOE OSTI2020

Unification of nonequilibrium molecular dynamics and the mode-resolved phonon Boltzmann equation for thermal transport simulations

Nano-size confinement induces many intriguing non-Fourier heat conduction phenomena, such as nonlinear temperature gradients, temperature jumps near the contacts, and size-dependent thermal conductivity. Over the past decades, these effects have been studied and interpreted by nonequilibrium molecular dynamics (NEMD) and phonon Boltzmann transport equation (BTE) simulations separately, but no theory that unifies these two methods has ever been established. In this work, we unify these methods using a quantitative mode-level comparison and demonstrate that they are equivalent for various thermostats. We show that different thermostats result in different non-Fourier thermal transport characteristics due to the different mode-level phonon excitations inside the thermostats, which explains the different size-dependent thermal conductivities calculated using different reservoirs, even though they give the same bulk thermal conductivity. Specifically, the Langevin thermostat behaves like a thermalizing boundary in phonon BTE and provides mode-level thermal-equilibrium phonon outlets, while the Nose-Hoover chain thermostat and velocity rescaling method behave like biased reservoirs, which provide a spatially uniform heat generation and mode-level nonequilibrium phonon outlets. These findings explain why different experimental measurement methods can yield different size-dependent thermal conductivity. They also indicate that the thermal conductivity of materials can be tuned for various applications by specifically designing thermostats. Finally, the unification of NEMD and phonon BTE will largely facilitate the study of thermal transport in complex systems in the future by, e.g., replacing computationally unaffordable first-principles NEMD simulations with computationally less expensive spectral BTE simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Quantum maximum entropy closure for small flavor coherence

Quantum angular moment transport schemes are an important avenue toward describing neutrino flavor mixing phenomena in dense astrophysical environments such as supernovae and merging neutron stars. Successful implementation will require new closure relations that go beyond those used in classical transport. In this paper, we derive the first analytic expression for a quantum M1 closure, valid in the limit of small flavor coherence, based on the maximum entropy principle. We verify that the resulting closure relation has the appropriate limits and characteristic speeds in the diffusive and free-streaming regimes. We then use this new closure in a moment linear stability analysis to search for fast flavor instabilities in a binary neutron star merger simulation and find better results as compared with previously designed, ad hoc , semiclassical closures.

astrophysical & cosmological simulations↗

NASA NTRS2003

Effect of Alignment on Transport Properties of Carbon Nanotube/Metallic Junctions

Ballistic and spin coherent transport in single walled carbon nanotubes (SWCNT) are predicted to enable high sensitivity single-nanotube devices for strain and magnetic field sensing. Based upon these phenomena, electron beam lithography procedures have been developed to study the transport properties of purified HiPCO single walled carbon nanotubes for development into sensory materials for nondestructive evaluation. Purified nanotubes are dispersed in solvent suspension and then deposited on the device substrate before metallic contacts are defined and deposited through electron beam lithography. This procedure produces randomly dispersed ropes, typically 2 - 20 nm in diameter, of single walled carbon nanotubes. Transport and scanning probe microscopy studies have shown a good correlation between the junction resistance and tube density, alignment, and contact quality. In order to improve transport properties of the junctions a technique has been developed to align and concentrate nanotubes at specific locations on the substrate surface. Lithographic techniques are used to define local areas where high frequency electric fields are to be concentrated. Application of the fields while the substrate is exposed to nanotube-containing solution results in nanotube arrays aligned with the electric field lines. A second electron beam lithography layer is then used to deposit metallic contacts across the aligned tubes. Experimental measurements are presented showing the increased tube alignment and improvement in the transport properties of the junctions.

Wincheski, Buzz↗

DOE OSTI2020

Transient Studies of Nonequilibrium Transport in Two-Dimensional Semiconductors

The overarching objective of research under DOE support has been to investigate emergent carrier phenomena in different two-dimensional (2D) systems, under conditions where they are driven far from equilibrium. The description of nonequilibrium transport has long represented one of the most enduring challenges in condensed-matter physics, in large part due to the need to invoke the role of complex many-body phenomena (electron-phonon and electron-electron coupling) in this regime. In this project, we have explored aspects of this problem in a variety of different two-dimensional systems. These include the modulation-doped two-dimensional electrons gas, or 2DEG, in GaAs/AlGaAs heterostructures; the electron and hole systems in graphene and its related hetero-structures, and; the strongly-confined electron layers in transition-metal dichalcogenide (MoS 2 and WS 2 ) and trichalcogenide (TiS 3 ) materials. Significant understanding has been achieved in a number of different areas, as we elaborate upon in subsequent sections of this proposal: Nanoscale spectroscopy of defect-related barriers in MoS2 Negative differential conductance due to inter-valley scattering in WS 2 under high fields Incipient localization due to scattering from frozen-phonon disorder in quantum point contacts High-field transient transport in graphene/hexagonal BN heterostructures Universal voltage scaling of quantum corrections in graphene under nonequilibrium A defining feature of much of this research has been an approach to problem solving in which novel experimental techniques (such as transient electrical pulsing and scanning-probe microscopy) are combined with the development (in collaboration with colleagues from condensed-matter theory) of sophisticated models, a strategy that we propose to continue in our future research. In recent years, we have developed an increasing capacity for handling atomically-thin materials, and their van der Waals systems, a capacity that is reflected in the balance of topics listed above.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS2021

Topical: Boiling, Condensation and Two-Phase Flows in Microgravity

The white paper addresses both fundamental and technological questions in the science of phase-change fluid flows aimed at understanding thermal transport in flow boiling and condensation and the use of these processes in space and terrestrial thermal systems. While flow boiling and condensation are widely employed in Earth-based thermal equipment to reduce the cost, size, and power required to transfer heat, a lack of data and predictive models on the impact of gravity on these phenomena hinders the use of two-phase heat transport in a low-gravity environment. The white paper presents research tasks in long-duration microgravity recommended to expand our knowledge on phase-changing flows. Due to elimination of the masking gravity effects, the proposed research will provide a unique opportunity to bring new insights into the role of capillary and hydrodynamic forces in phase-changing flows. The potential impact of the proposed research ranges from the space power, thermal and cryogenic systems to various terrestrial applications in the petrochemical, pharmaceutical, biochemical, nuclear, and metallurgical industries.

Phase-change fluid flow↗

DOE OSTI2021

Growth, characterization, and Chern insulator state in MnBi 2⁢ Te 4 via the chemical vapor transport method

As the first intrinsic antiferromagnetic topological insulator, MnBi 2 Te 4 has provided a platform to investigate the interplay of band topology and magnetism as well as the emergent phenomena arising from such an interplay. Here we report the chemical-vapor-transport (CVT) growth and characterization of MnBi 2 Te 4 , as well as the observation of the field-induced quantized Hall conductance in 6-layer devices. Through comparative studies between our CVT-grown and flux-grown MnBi 2 Te 4 via magnetic, transport, scanning tunneling microscopy, and angle-resolved photoemission spectroscopy measurements, we find that CVT-grown MnBi 2 Te 4 is marked with higher Mn occupancy on the Mn site, slightly higher Mn Bi antisites, smaller carrier concentration, and a Fermi level closer to the Dirac point. Furthermore, a 6-layer device made from the CVT-grown sample shows by far the highest mobility of 2500 cm 2 V s in MnBi 2 Te 4 devices with the quantized Hall conductance appearing at 1.8 K and 8 T. Furthermore, our study provides a route to obtain high-quality single crystals of MnBi 2 Te 4 that are promising to make superior devices and realize emergent phenomena, such as the layer Hall effect and quantized anomalous Hall effect, etc.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Electrostatic limits on charge selectivity in ultrahigh charge density polymer membranes

Ion-exchange membranes (IEMs) with ultrahigh charge densities offer the promise of enhanced ion transport for electrochemical technologies, yet the fundamental limits of this membrane design strategy are not yet understood. In this work, we present a systematic study of bis(1-vinyl-3-imidazolium) cross-linked polymer membranes with fixed charge contents ranging from 6 to 9 mol/L[dry polymer], synthesized to isolate the effects of charge density at constant membrane water content. While ionic conductivity increases monotonically with increasing charge density, a sharp decline in charge selectivity is observed for the most densely charged membranes, defying conventional expectations. Structural, thermal, and mechanical analyses reveal a critical onset of network disruption and anomalous ion partitioning behavior in these densely charged membranes. With proper consideration of cross-linker geometry, these results are interpreted using Manning's counter-ion condensation theory, which suggests that closely packed charged cross-linkers amplify inter-chain electrostatic interactions and trigger excess counter-ion condensation. Our findings suggest that, beyond a threshold spacing between fixed charges, increased functionalization may harm rather than help charge selectivity of IEMs.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Singular Hall Response from a Correlated Ferromagnetic Flat Nodal‐Line Semimetal

Abstract Topological quantum phases are largely understood in weakly correlated systems, which have identified various quantum phenomena, such as the spin Hall effect, protected transport of helical fermions, and topological superconductivity. Robust ferromagnetic order in correlated topological materials particularly attracts attention, as it can provide a versatile platform for novel quantum devices. Here, a singular Hall response arising from a unique band structure of flat topological nodal lines in combination with electron correlation in a van der Waals ferromagnetic semimetal, Fe 3 GaTe 2 , with a high Curie temperature ofT c = 347 K is reported. High anomalous Hall conductivity violating the conventional scaling, resistivity upturn at low temperature, and a large Sommerfeld coefficient are observed in Fe 3 GaTe 2 , which implies heavy fermion features in this ferromagnetic topological material. The scanning tunneling microscopy, circular dichroism in angle‐resolved photoemission spectroscopy, and theoretical calculations support the original electronic features of the material. Thus, low‐dimensional Fe 3 GaTe 2 with electronic correlation, topology, and room‐temperature ferromagnetic order appears to be a promising candidate for robust quantum devices.

Chemistry↗

DOE OSTI2022

Nonadiabatic Molecular Dynamics Simulations Based on Time-Dependent Density Functional Tight-Binding Method

Nonadiabatic excited-state molecular dynamics underpin many photophysical and photochemical phenomena, such as exciton dynamics, charge separation and transport. In this work, we present an efficient nonadiabatic molecular dynamic (NAMD) simulation method based on time-dependent density functional tight-binding (TDDFTB) theory. Specifically, the adiabatic electronic structure, an essential NAMD input, is described at the TDDFTB level. Here, the nonadiabatic effects originating from the coupled motions of electrons and nuclei are treated by the trajectory surface hopping algorithm. To improve the computational efficiency, nonadiabatic couplings between excited states within the TDDFTB method are derived and implemented using an analytical approach. Further, the time-dependent nonadiabatic coupling scalars are calculated based on the overlap between molecular orbitals rather than the Slater determinants to speed up the simulations. In addition, the electronic decoherence scheme and a state reassigned unavoided crossings algorithm, which has been implemented in the NEXMD software, are used to improve the accuracy of the simulated dynamics and handle trivial unavoided crossings. Finally, the photoinduced nonadiabatic dynamics of a benzene molecule are simulated to demonstrate our implementation. The results for excited state NAMD simulations of benzene molecule based on TDDFTB method compare well that obtained with numerically expensive time-dependent density functional theory. The proposed methodology provides an attractive theoretical simulation tool for predicting the photophysical and photochemical properties of complex materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗