Single Energy Research Articles

Computational Assessment of Clinical Drugs against SARS-CoV-2: Foreseeing Molecular Mechanisms and Potent Mpro Inhibitors.

The emergence of new SARS-CoV-2 variants of concern (VOC) is a propulsion for accelerated potential therapeutic discovery. SARS-CoV-2's main protease (Mpro), essential for host cell viral replication, is a pre-eminent druggable protein target. Here, we perform extensive drug re-profiling of the comprehensive Excelra database, which compiles various under-trial drug candidates for COVID-19 treatment. For mechanistic understanding, the most promising screened-out molecules with targets are subjected to molecular dynamics (MD) simulations. Post-MD analyses demonstrate Darunavir, Ponatinib, and Tomivosertib forming a stable complex with Mpro, characterized by less fluctuation of Cα atoms, smooth and stable root-mean-square deviation (RMSD), and robust contact with the active site residues. Likewise, they all have lower binding free energy with Mpro, demonstrating strong affinity. In free energy landscape profiles, the distances from His41 and Cys145 exhibit a single energy minima basin, implying their preponderance in proximity to Mpro's catalytic dyad. Overall, the computational assessment earmarks promising candidates from the Excelra database, emphasizing on carrying out exhaustive biochemical experiments along with clinical trials. The work lays the foundation for potential therapeutic interventions in treating COVID-19.

Simulation and Experiments on Optimization of Vortex-Induced Vibration Power Generation System Based on Side-by-Side Double Blunt Bodies

In order to improve the utilization efficiency of converting low-flow current energy into electric energy for Reynolds number 10,000 ≤ Re ≤ 40,000, this paper proposes a vortex-induced vibration power generation system based on a side-by-side double blunt body. In this system, the side-by-side double blunt body structure is used in the current energy capture part to enhance the collection of low-flow current energy; the permanent magnet linear motor is used in the electric energy conversion part to improve the efficiency of electric energy conversion; and a laboratory device is constructed for testing. The effects of the blunt body structure parameters and the center spacing ratio on the energy harvesting performance of the system are qualitatively explained by constructing a simulation model. Compared with the single blunt body energy capture structure, the side-by-side double blunt body structure increases the vibration amplitude by 1.04 times and the lift by 1.14 times at the center spacing S/D = 2.4. Meanwhile, energy harvesting can be realized at a lower flow velocity, increasing the vortex-induced vibration’s energy capture range. Finally, the power generation system was experimentally verified in the laboratory, and the results showed that the vibration amplitude of the double blunt body structure was increased by 1.12 times compared to the single blunt body. The maximum output power of the generator is 10.55 W when the water velocity is 0.7 m/s. The energy conversion efficiency of the power generation system can reach a maximum of 52.93%, which is 12.33% higher than that of a single blunt body structure, which proves that the system has a higher power conversion efficiency than that of a conventional single conversion system.

Towards a complete description of the reaction mechanisms between nitrenium ions and water.

Nitrenium ions are intermediates in the metabolic routes producing the highly carcinogenic nitrosamines and binding to DNA molecules. The reaction mechanism of nitrenium molecules with explicit water molecules is sensibly dependent on the number of waters: when a second molecule is involved, it acts as a catalyst for the reaction, lowering intrinsic activation barriers regardless of the substituent. For all cases, the reaction force constants and reaction electron flux indicate highly synchronous reactions for . Conversely, for highly non-synchronous reactions are obtained, involving two separate proton transfers happening early and late in the reaction path. As a test case, for the simplest reactions, orbital interactions within the NBO paradigm, bond orders, and their derivatives indicate that each individual proton transfer is highly synchronous. Molecular geometries were optimized and characterized at the B3LYP/6-311++G(d, p) level. Intrinsic reaction coordinates were calculated. CCSD(T) single point energies with the same basis were computed on all stationary points. The reaction force, reaction force constant, and reaction electron flux are used to study the evolution of the reacting systems. Natural bond orbitals are used to understand the primitive changes driving the reaction.

Effective remediation of toxic metal ions (Cd(II), Pb(II), Hg(II), and Ba(II)) using mesoporous glauconite-based iron silicate nanorods: Experimental and theoretical studies

Glauconite minerals underwent an advanced exfoliation and scrolling process, yielding novel iron silicate nanorods (GRs) with increased surface area, reactivity, and improved physicochemical properties. These structures were introduced as superior adsorbents for the highly efficient adsorption of various toxic metal ions, such as Cd(II), Pb(II), Hg(II), and Ba(II). The GRs exhibited maximum adsorption capacities of 283 mg/g for Cd(II), 247 mg/g for Pb(II), 132.3 mg/g for Hg(II), and 165.2 mg/g for Ba(II). The adsorption properties of the GRs during the adsorption of these four metal ions were elucidated through traditional (Langmuir model) and advanced (monolayer model of single energy site) isotherm analyses. Advanced isotherm modeling revealed that the GR surface was saturated with numerous effective adsorption sites, with densities of 125 mg/g for Cd(II), 68.8 mg/g for Pb(II), 40.9 mg/g for Hg(II), and 57.9 mg/g for Ba(II). Moreover, each site could accommodate approximately four ions of the studied metals, which are vertically oriented and participate in multi-ionic adsorption reactions. Energetic analyses, whether based on classical models (Gaussian energy <8 kJ/mol) or advanced models (adsorption energy <40 kJ/mol), indicated that the adsorption processes are governed by physical mechanisms, including electrostatic attractions, van der Waals forces, and hydrogen bonding. Furthermore, thermodynamic assessments confirmed that the adsorption of these ions occurs through exothermic and spontaneous reactions.

Proton dose calculation with LSTM networks in presence of a magnetic field

Objective.To present a long short-term memory (LSTM) network-based dose calculation method for magnetic resonance (MR)-guided proton therapy.Approach.35 planning computed tomography (CT) images of prostate cancer patients were collected for Monte Carlo (MC) dose calculation under a perpendicular 1.5 T magnetic field. Proton pencil beams (PB) at three energies (150, 175, and 200 MeV) were simulated (7560 PBs at each energy). A 3D relative stopping power cuboid covering the extent of the PB dose was extracted and given as input to the LSTM model, yielding a 3D predicted PB dose. Three single-energy (SE) LSTM models were trained separately on the corresponding 150/175/200 MeV datasets and a multi-energy (ME) LSTM model with an energy embedding layer was trained on either the combined dataset with three energies or a continuous energy (CE) dataset with 1 MeV steps ranging from 125 to 200 MeV. For each model, training and validation involved 25 patients and 10 patients were for testing. Two single field uniform dose prostate treatment plans were optimized and recalculated with MC and the CE model.Results.Test results of all PBs from the three SE models showed a mean gamma passing rate (2%/2 mm, 10% dose cutoff) above 99.9% with an average center-of-mass (COM) discrepancy below 0.4 mm between predicted and simulated trajectories. The ME model showed a mean gamma passing rate exceeding 99.8% and a COM discrepancy of less than 0.5 mm at the three energies. Treatment plan recalculation by the CE model yielded gamma passing rates of 99.6% and 97.9%. The inference time of the models was 9-10 ms per PB.Significance.LSTM models for proton dose calculation in a magnetic field were developed and showed promising accuracy and efficiency for prostate cancer patients.

π-Bond Dissociation Energies: C-C, C-N, and C-O.

Single, double, and triple bond dissociation energies are useful quantities, but students and practicing chemists commonly do not know representative values or how they are obtained. In this paper, select π-bond energies are provided and general methods for their determination are discussed. Relationships between heats of hydrogenation, bond dissociation energies, and π-bond strengths are also addressed.

A Fast 3D Range-Modulator Delivery Approach: Validation of the FLUKA Model on a Varian ProBeam System Including a Robustness Analysis

A 3D range-modulator (RM), optimized for a single energy and a specific target shape, is a promising and viable solution for the ultra-fast dose delivery in particle therapy. The aim of this work was to investigate the impact of potential beam and modulator misalignments on the dose distribution. Moreover, the FLUKA Monte Carlo model, capable of simulating 3D RMs, was adjusted and validated for the 250 MeV single-energy proton irradiation from a Varian ProBeam system. A 3D RM was designed for a cube target shape rotated 45° around two axes using a Varian-internal research version of the Eclipse treatment planning software, and the resulting dose distribution was simulated in a water phantom. Deviations from the ideal alignment were introduced, and the dose distributions from the modified simulations were compared to the original unmodified one. Finally, the FLUKA model and the workflow were validated with base-line data measurements and dose measurements of the manufactured modulator prototype at the HollandPTC facility in Delft. The adjusted FLUKA model, optimized particularly in the scope of a single-energy FLASH irradiation with a PMMA pre-absorber, demonstrated very good agreement with the measured dose distribution resulting from the 3D RM. Dose deviations resulting from modulator-beam axis misalignments depend on the specific 3D RM and its shape, pin aspect ratio, rotation angle, rotation point, etc. A minor modulator shift was found to be more relevant for the distal dose distribution than for the spread-out Bragg Peak (SOBP) homogeneity. On the other hand, a modulator tilt (rotation away from the beam axis) substantially affected not only the depth dose profile, transforming a flat SOBP into a broad, Gaussian-like distribution with increasing rotation angle, but also shifted the lateral dose distribution considerably. This work strives to increase awareness and highlight potential pitfalls as the 3D RM method progresses from a purely research concept to pre-clinical studies and human trials. Ensuring that gantry rotation and the combined weight of RM, PMMA, and aperture do not introduce alignment issues is critical. Given all the other range and positioning uncertainties, etc., not related to the modulator, the RM must be aligned with an accuracy below 1° in order to preserve a clinically acceptable total uncertainty budget. Careful consideration of critical parameters like the pin aspect ratio and possibly a novel robust modulator geometry optimization are potential additional strategies to mitigate the impact of positioning on the resulting dose. Finally, even the rotated cube 3D modulator with high aspect ratio pin structures (~80 mm height to 3 mm pin base width) was found to be relatively robust against a slight misalignment of 0.5° rotation or a 1.5 mm shift in one dimension perpendicular to the beam axis. Given a reliable positioning and QA concept, the additional uncertainties introduced by the 3D RM can be successfully managed adopting the concept into the clinical routine.

Dispersion Parameters of Copper Sulphate Doped PMMA

Films of PMMA and copper sulphate doped PMMA have been prepared by casting method. Absorbance and transmittance spectra were recorded in the wavelength range (300-900) nm in order to calculate, single oscillator energy, dispersion energy, average oscillator strength, the refractive index at infinite wavelength, M-1 and M -3 moments of the optical spectra, it was found that all these parameters were effected by doping.

On the Use of the Multi-Site Langmuir Model for Predicting Methane Adsorption on Shale

Shale gas, mainly consisting of adsorbed gas and free gas, has served a critical role of supplying the growing global natural gas demand in the past decades. Considering that the adsorbed methane has contributed up to 80% of the total gas in place (GIP), understanding the methane adsorption behaviors is imperative to an accurate estimation of total GIP. Historically, the single-site Langmuir model, with the assumption of a homogeneous surface, is commonly applied to estimate the adsorbed gas amount. However, this assumption cannot depict the methane adsorption characteristics due to various compositions and pore sizes of shales. In this work, a multi-site model integrating the energetic heterogeneity in adsorption is derived to predict methane adsorption on shale. Our results show that the multi-site model is capable of addressing the heterogeneity of shales by a wide range of adsorption energy distributions (owing to the complex compositions and different pore sizes), which is different from the single-site model only characterized by single adsorption energy. Consequently, the multi-site model results have better accuracy against the experimental data. Therefore, applying the multi-site Langmuir model for estimating GIP in shales can achieve more accurate results compared with using the traditionally single-site model.

An ANI-2 enabled open-source protocol to estimate ligand strain after docking.

In protein-ligand docking, the score assigned to a protein-ligand complex is approximate. Especially, the internal energy of the ligand is difficult to compute precisely using a molecular mechanics based force-field, introducing significant noise in the rank-ordering of ligands. We propose an open-source protocol (https://github.com/UnixJunkie/MMO), using two quantum mechanics (QM) single point energy calculations, plus a Monte Carlo (Monte Carlo) based ligand minimization procedure in-between, to estimate ligand strain after docking. The MC simulation uses the ANI-2x (QM approximating) force field and is performed in the dihedral space. On some protein targets, using strain filtering after docking allows to significantly improve hit rates. We performed a structure-based virtual screening campaign on nine protein targets from the Laboratoire d'Innovation Thérapeutique-PubChem assays dataset using Cambridge crystallographic data centre genetic optimization for ligand docking. Then, docked ligands were submitted to the strain estimation protocol and the impact on hit rate was analyzed. As for docking, the method does not always work. However, if sufficient active and inactive molecules are known for a given protein target, its efficiency can be evaluated.

Ablation efficiency and laser safety of a novel superpulsed thulium fiber laser: a in vitro study.

To assess the ablation efficiency of the Superpulsed Thulium Fiber Laser (SP-TFL) and investigate the thermal effects of SP-TFL. A SP-TFLwas employed to evaluate ablation efficiency. Fresh ex-vivo pig kidneys and ureters were utilized to evaluate the renal pelvis and ureter temperature changes, different irrigation rates(0, 15, 38mL/min) and a long pulse width were used. The research indicated that as laser output power increased, ablation rates significantly increased. Ablation rates(mg/min) were higher and the energy per ablated mass(J/mg) was lower at lower frequencies(10-50Hz). Under the same frequency and single pulse energy, super short and short pulse widths demonstrated higher ablation rates at higher frequencies (exceeding 100Hz). The temperature of the renal pelvis and ureter decreased with increasing irrigation rates. In the renal pelvis, without irrigation, the temperature quickly reached the critical threshold of 43℃. The irrigation rate was 15ml/min and power was no more than 18W, the renal pelvis temperature did not reach 43℃. When the irrigation rate were 38ml/min, the temperature did not risen to 43℃. In the ureter, without irrigation, the temperature also quickly reached 43℃. The temperature reached 43℃ when the power exceeded to12W with an irrigation rate of 15ml/min. With an irrigation rate of 38ml/min, the temperature reached 43℃ at a laser power of 30W. The SP-TFL demonstrated promising ablation effectiveness especially for lower frequencies and super short and short pulse widths model. Proper irrigation rates, single pulse energy, frequency and pulse width are crucial during lithotripsy.

Dissipative soliton resonance in a normal dispersion all-polarization-maintaining thulium-doped fiber laser

We report a first demonstration of dissipative soliton resonance (DSR) in a normal dispersion thulium-doped fiber laser (TDFL) constructed entirely of polarization-maintaining fibers. A nonlinear fiber loop mirror with a combination of normal and anomalous dispersion fibers inside the loop ensures a self-starting mode-locking operation in the DSR regime. The oscillator generated linearly polarized optical pulses with a fundamental repetition frequency of 5.735 MHz and a nearly constant peak power clamped at ∼ 4.35 W. The pulse duration extends from 217 ps to 6 ns, when the pump power increases from 417 mW to 1.86 W. The maximum average output power of the DSR pulses was 151 mW, corresponding to the single pulse energy of 26 nJ.

High power continuous-wave, Q-switched, and quasi-continuous-wave operation of a diode-pumped Tm,Ho:YAG laser oscillator at 2.09 μm

A high-power diode-pumped 2.09 μm Tm,Ho:YAG laser oscillator with multi-regime operation is demonstrated. In the continuous-wave (CW) regime, it delivers 44.68 W output power with M2 = 3.47, yielding a record-high brightness of 87.57 MW/cm2·sr. In the Q-switched mode, stable laser performance across pulse repetition frequencies (PRFs) from 1 kHz to 4 kHz is achieved. At a PRF of 1 kHz, single pulse energy reaches 22.85 mJ with a pulse width of 443.4 ns and M2 = 3.50. In the quasi-continuous-wave (QCW) operation, a pulse energy of 31.85 mJ is obtained at 200 Hz with a 370 μs pulse width and M2 = 4.49. These pulse energies are the highest reported to date for Tm-Ho co-doped laser oscillators at high PRF. This also marks the first demonstration of a versatile 2 μm laser that can operate at room temperature while delivering high-energy, high-repetition-rate pulses in both nanosecond and microsecond durations within a single device, making it highly appealing for scientific and industrial applications.

Design and performance evaluation of a solar-assisted biogas and air-source heat pump hybrid system for radiant heating applications

Energy and environmental crises have become critical challenges, making the efficient use of clean energy a focus of public attention. Although emerging clean energies have high economic and environmental values, these single energy sources are not well-developed. In colder climates, maintaining a stable fermentation temperature can be difficult, leading to the idling of many biogas digesters. Furthermore, despite extensive research on solar-assisted air-source heat pumps, these systems perform poorly in cold climates with low solar irradiance. To address these issues, this study investigates the integration of a solar-assisted biogas and air-source heat pump hybrid heating system (SBHP-HHS) for capillary radiant heating in rural North China, specifically adapted to local conditions. The system aims to utilize the combined benefits of various clean energy sources to improve system reliability in harsh climates. By combining experimental analysis with TRNSYS simulation, the system performance during the heating season in rural northern China was evaluated. Experimental results show that the hybrid system, utilizing the radiant capillary terminal, is capable of maintaining room temperatures above 20 °C, demonstrating its effectiveness. Simulation results show that the air-source heat pump achieves a coefficient of performance (COP) of 4.6. This efficiency level is generally favorable, ensuring stable operation. Additionally, the system achieves a solar-biogas hybrid utilization rate exceeding 75 % and a primary energy-saving rate of 78 %. Economically, the system offers an 8.23-year dynamic payback period and annual cost savings of 3055 CNY. Environmentally, it reduces CO2 emissions by 4263.6 kg annually. This study provides a solid experimental foundation for further research, offering empirical evidence to optimize system performance.

Inferring ion cluster lifetimes from energy-resolved mass spectra of an EMI-BF4 electrospray plume

The fragmentation of ion clusters within the accelerating fields of ionic liquid ion sources (ILISs) is well documented and degrades ILIS performance and lifetime. Some of the most popular ILIS liquids, such as EMI-BF4 (1-Ethyl-3-methylimidazolium tetrafluoroborate) and EMI-Im (1-Ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide), emit clusters with lifetimes as low as ∼1 ns. Studies of fragmentation within the accelerating field typically rely on measuring the plume energy distribution averaged over all plume species and comparing those measurements with numerical simulations to estimate ion cluster lifetimes. Here, for the first time, we estimate EMI-BF4 cluster lifetimes by analyzing the energy distributions of individual plume species. We use this novel analysis method to estimate mean lifetimes of positive EMI-BF4 ion clusters from previously published experimental data. We find that the mean lifetime ranges from τ=3.7ns to τ=124ns for [EMI+][EMI-BF4] dimers and ranges from τ=1.5ns to τ=23ns for [EMI+][EMI-BF4]2 trimers. Fitting those data to an analytical fragmentation model, we estimate the binding energy and temperature as ΔGS0=0.49eV and T=394K for dimers and ΔGS0=0.40eV and T=365K for trimers. Comparing our results with previous studies supports the conclusion that clusters are emitted with a wide distribution of internal energies, contrary to the common assumption of single internal energy for each species.

Performance Study and Application of Antioxidant Peptides CHCIWM and CHCPWM

AbstractWhen the level of reactive oxygen species (ROS) increases abnormally in the body, it damages tissue cells and is closely related to the occurrence of various chronic diseases. Antioxidant peptides (APs) are a class of active short peptides with antioxidant capacity, which have great research value. In this paper, two APs were designed and synthesized based on the structural features of APs: CHIWM (AP1) and CHCPWM (AP2). The reactive sites and activities of APs were analyzed using the quantum chemical density functional theory (DFT) method. The reactive sites are all located in the indolyl heterocycle of tryptophan residues, among which AP1 has a smaller energy gap and single point energy. Multiple free radical scavenging assays showed that AP1 and AP2 had stronger antioxidant activity than the positive control, glutathione (GSH). For scavenging DPPH free radicals, the IC50s values were 0.613 mg/mL and 0.606 mg/mL, respectively. In cellular assays, both AP1 and AP2 were able to scavenge ROS in the cells and attenuate cellular damage. In conclusion, AP1 and AP2 designed in this paper have good antioxidant activity and have the potential to be applied in the treatment of chronic diseases caused by cellular oxidative damage.

Characterization of Mn3O4/(x)Nb2O5 thin films as a promising material for supercapacitors

The thin films of Mn3O4/(x)Nb2O5 oxides (x = 0, 2, 4, and 6 %) were prepared by the electron beam evaporation method and then annealed at 400 °C. The structural analysis showed a tetragonal phase of Mn3O4, according to X-ray diffraction results. The quantitative elemental analysis of films confirms the stoichiometry of the Mn3O4, while the Nb2O5 couldn’t be recognized using energy dispersive X-ray fluorescence. The chemical speciation of the films was investigated using X-ray photoelectron spectroscopy. Three oxidation states of Mn were recognized at average binding energies of 645.70±0.35, 642.99±0.40, and 641.58±0.48 eV that correlated to Mn4+, Mn3+, and Mn2+, respectively. As the Nb2O5 increases, a slight decrease in the binding energy of Mn4+ is recognized, whereas it is nearly constant for the deconvoluted Mn3+ and Mn2+. Two electronic transitions (Eg1 and Eg2) of the optical band gap showed an increase with increasing Nb2O5 dopant from 2 % to 6 %. The optical parameters such as Eosc, the single oscillator energy, Edis, the dispersion energy, both moments (M−1) and (M−3), the ratio of carrier concentration/effective mass (N/m*), and the lattice dielectric constant (εL) were determined. From electrochemical measurements of pure Mn3O4 films, it is found that the voltammetric current density is directly proportional to the scan rates of cyclic voltammetry CV, with good cyclic stability, and by increasing the scan rates, the values of specific capacitance decrease while the values of specific power increase. In doped thin films, no potential drop is observed in the discharging process, which indicates good interfacial contact between the active material and the substrate during the charge/discharge process. These results suggest that the obtained Mn3O4 nanoparticles are a better candidate for supercapacitor applications.

President’s Column: SPE@Energy Transition

_ Editor’s Note: This is a summary of the October episode of the President’s vodcast. We encourage you to watch the episode to see the full conversation. In this vodcast, Houzé discusses his views on SPE and the energy transition. _ Last month, I presented a to-do list with five main topics: energy transition, finances, governance, membership, and quality. This month, I will discuss the first topic, energy transition. 1. SPE’s polarized views on energy transition As I mentioned last month, there are significant geographic contrasts within our industry and SPE regarding how the energy transition is perceived and developed. In some SPE venues and regions, energy transition is the dominant, if not the only subject, and the P in SPE is indeed a major obstacle for support, even from our traditional stakeholders. In other regions, stakeholders acknowledge the energy transition, but the situation is essentially business as usual. Major investments are underway to sustain production at current levels or even increase hydrocarbon production to record levels, with development plans spanning over several decades. Even in these regions, most, if not all, new projects are designed from the outset to have a low carbon footprint, including mitigating methane emissions. Energy transition is somehow factored in. So, there are very different narratives, which are reflected in a polite but real polarization within SPE—geographic and, occasionally, demographic. As an SPE President visiting many sections and stakeholders, the contrast is often striking. But we are engineers operating on the same planet, and there must be a single energy narrative that accounts for the reality of the energy transition in whatever form or definition this takes, locally or globally. To do this, we can try to assess in the most objective way possible how our license to operate may evolve in the short to medium term. I will start with something you may have seen many times: the three main energy transition scenarios presented by the International Energy Agency (IEA). These scenarios predict the climate impact by the end of the century based on current and predicted variables, including consumption, population, economics, and technologies. These are just models, but as the statistician George Box pointed out: “All models are wrong, but some are useful.” These IEA models are definitely useful (Fig. 1).

1 kHz, 5.3 kW peak power pulse generation from a LiNbO3 acousto-optically Q-switched Er:YAP laser

In this paper, we present, for the first time, a bulk LiNbO3 (LN) acousto-optic (AO) Q-switched Er:YAP laser operating at a pulse repetition frequency (PRF) of 1 kHz. At this frequency, the laser achieves a maximum single pulse energy of 0.34 mJ, with a corresponding pulse duration of 64.5 ns and a peak power of 5.3 kW. Additionally, we characterize the laser beam quality factors, Mx2 = 2.7 and My2 = 2.5, and its central wavelength at 2730.7 nm, all at 1 kHz. These results highlight the significant potential of LN AO Q-switched Er:YAP lasers for enhancing mid-infrared laser applications.

Flexible Thermoelectric and Piezoelectric Hybrid Energy harvester by Adopting a Multifunctional Cellulose-based Composite Film

Thermoelectric (TE) and piezoelectric (PE) energy harvesting are the most advantageous conversion mechanisms for converting energy generated by the human body into electrical energy, respectively. However, the development of a hybrid energy harvester is essential for achieving high energy conversion efficiency since a single energy harvester cannot convert enough energy. Herein, we present a flexible TE-PE hybrid energy harvester (f-TPHEH) based on a self-assembled single film consisting of Bi2Te2.7Sn0.3 particles, BaTiO3 nanoparticles, and cellulose. The fabricated f-TPHEH exhibited a maximum TE output power of 18.47 nW and a PE instantaneous output of 13.5 nW under integrated thermal and bending conditions (ΔT = 40 K, bending displacement = 7 mm). We conducted a theoretical analysis using multiphysics simulation based on finite element analysis to support the experimental measurements of the fabricated f-TPHEH. In addition, to demonstrate the mechanical stability of the cellulose matrix-based f-TPHEH, a durability test was performed up to 5,000 cycles of bending deformation. This study presents a method for developing the flexible hybrid energy harvester using a low-cost and simple fabrication process.