Synthetic, EPR spectroscopic, magnetic and X-ray crystallographic structural studies on copper(II) complexes of the tridentate N2S donor ligand formed from 6-methyl-2-formylpyridine and S-methyldithiocarbazate (Hmpsme)

New copper(II) complexes of general empirical formula, Cu(mpsme)X center dot xCH(3)COCH(3) (mpsme = anionic form of the 6-methyl-2-formylpyridine Schiff base of S-methyldithiocarbazate; X = Cl, N-3, NCS, NO3; x = 0, 0.5) have been synthesized and characterized by IR, electronic, EPR and susceptibility measurements. Room temperature mu(eff) values for the complexes are in the range 1.75-2.1 mu(beta) typical of uncoupled or weakly coupled Cu(II) centres. The EPR spectra of the [Cu(mpsme)X] (X = Cl, N-3, NO3, NCS) complexes reveal a tetragonally distorted coordination sphere around the mononuclear Cu(II) centre. We have exploited second derivative EPR spectra in conjunction with Fourier filtering (sine bell and Hamming functions) to extract all of the nitrogen hyperfine coupling matrices. While the X-ray crystallography of [Cu(mpsme)NCS] reveals a linear polymer in which the thiocyanate anion bridges the two copper(II) ions, the EPR spectra in solution are typical of a magnetically isolated monomeric Cu(II) centres indicating dissociation of the polymeric chain in solution. The structures of the free ligand, Hmpsme and the {[Cu(mpsme)NO3] center dot 0.5CH(3)COCH(3)}(2) and [Cu(mpsme)NCS](n) complexes have been determined by X-ray diffraction. The {[Cu(mpsme)NO3]0.5CH(3)COCH(3)}(2) complex is a centrosymmetric dimer in which each copper atom adopts a five-coordinate distorted square-pyramidal geometry with an N2OS2 coordination environment, the Schiff base coordinating as a uninegatively charged tridentate ligand chelating through the pyridine and azomethine nitrogen atoms and the thiolate, an oxygen atom of a unidentate nitrato ligand and a bridging sulfur atom from the second ligand completing the coordination sphere. The [Cu(mpsme)(NCS)](n) complex has a novel staircase-like one dimensional polymeric structure in which the NCS- ligands bridge two adjacent copper(II) ions asymmetrically in an end-to-end fashion providing its nitrogen atom to one copper and the sulfur atom to the other. (c) 2005 Elsevier B.V. All rights reserved.

Preparation of high purity ZnO nanobelts by thermal evaporation of ZnS

High purity one-dimensional ZnO nanobelts were synthesized by thermally evaporating commercial ZnS powders in a hydrogen-oxygen mixture gas at 1050 degrees C. It was found that these ZnO nanobelts had a single crystal hexagonal wurtzite structure growing along the [0001] direction. They had a rectangle-shaped cross-section with typical widths of 20 to 100 nanometers and lengths of up to hundreds of micrometers with lattice constants of a = 0.325 nm and c = 0.520 nm. The self-catalytic hydrogen-oxygen assisted growth of ZnO nanobelt is discussed. The photoluminescence (PL) characterization of the ZnO nanobelts shows strong near-band UV emission (about 383 nm) and one broad peak at 501 nm, which indicates that the ZnO nanobelts have good potential application in optoelectronic devices.

Molecular architecture influences on material properties of pharmaceutical compounds

Salt formation has extensively been studied as a strategy to improve drug solubility but it has not been explored as a strategy to improve mechanical properties. A better understanding of which factors of the solid state can have an influence in the mechanical properties of pharmaceutical powders can help to optimise and reduce cost of tablet manufacturing. The aim of this study was to form different series of amine salts of flurbiprofen, gemfibrozil and diclofenac and to establish predictive relationships between architectural characteristics and physicochemical and mechanical properties of the salts. For this purpose, three different carboxylic acid drugs were selected: flurbiprofen, gemfibrozil and diclofenac, similar in size but varying in flexibility and shape and three different series of counterions were also chosen: one with increasing bulk and no hydroxyl groups to limit the hydrogen bonding potential; a second one with increasing number of hydroxyl groups and finally a third series, related to the latter in number of hydroxyl groups but with different molecular shape and flexibility. Physico-chemical characterization was performed (DSC, TGA, solubility, intrinsic dissolution rate, particle size, true density) and mechanical properties measured using a compaction replicator. Strained molecular conformations produce weaker compacts as they have higher energy than preferred conformations that usually lie close to energy minimums and oppose plastic deformation. It was observed that slip planes, which correspond to regions of weakest interaction between the planes, were associated with improved plasticity and stronger compacts. Apart from hydrogen bonds, profuse van der Waals forces can result in ineffective slip planes. Salts displaying two-dimensional densely hydrogen bonded layers produced stronger compacts than salts showing one-dimensional networks of non-bonded columns, probably by reducing the attachment energy between layers. When hydrogen bonds are created intramolecularly, it is possible that the mechanical properties are compromised as they do not contribute so much to create twodimensional densely bonded layers and they can force molecules into strained conformations. Some types of hydrogen bonding network may be associated with improved mechanical properties, such as type II, or R (10) 3 4 using graph-set notation, versus type III, or R (12) 4 8 , columns. This work clearly demonstrates the potential of investigating crystal structure-mechanical property relationship in pharmaceutical materials.

Study of neutral Fe(III) complexes of pyridoxal-N-substituted thiosemicarbazone with desolvation induced spin-state transformation above room temperature

The preparation and characterization of two new neutral ferric complexes with desolvation-induced discontinuous spin-state transformation above room temperature are reported. The compounds, Fe(Hthpy)(thpy).CH3OH.3H2O (1) and Fe(Hmthpy)(mthpy).2H2O (2), are low-spin (LS) at room temperature and below, whereas their nonsolvated forms are high-spin (HS), exhibiting zero-field splitting. In these complexes, Hthpy, Hmthpy, and thpy, mthpy are the deprotonated forms of pyridoxal thiosemicarbazone and pyridoxal methylthiosemicarbazone, respectively; each is an O,N,S-tridentate ligand. The molecular structures have been determined at 100(1) K using single-crystal X-ray diffraction techniques and resulted in a triclinic system (space group P1) and monoclinic unit cell (space group P21/c) for 1 and 2, respectively. Structures were refined to the final error indices, where RF = 0.0560 for 1 and RF = 0.0522 for 2. The chemical inequivalence of the ligands was clearly established, for the "extra" hydrogen atom on the monodeprotonated ligands (Hthpy, Hmthpy) was found to be bound to the nitrogen of the pyridine ring. The ligands are all of the thiol form; the doubly deprotonated chelates (thpy, mthpy) have C-S bond lengths slightly longer than those of the singly deprotonated forms. There is a three-dimensional network of hydrogen bonds in both compounds. The discontinuous spin-state transformation is accompanied with liberation of solvate molecules. This is evidenced also from DSC analysis. Heat capacity data for the LS and HS phases are tabulated at selected temperatures, the values of the enthalpy and entropy changes connected with the change of spin state were reckoned at DeltaH = 12.5 0.3 kJ mol-1 and DeltaS = 33.3 0.8 J mol-1 K-1, respectively, for 1 and DeltaH = 6.5 0.3 kJ mol-1 and DeltaS = 17.6 0.8 J mol-1 K-1, respectively, for 2

A method of fundamental solutions for two-dimensional heat conduction

We investigate an application of the method of fundamental solutions (MFS) to heat conduction in two-dimensional bodies, where the thermal diffusivity is piecewise constant. We extend the MFS proposed in Johansson and Lesnic [A method of fundamental solutions for transient heat conduction, Eng. Anal. Bound. Elem. 32 (2008), pp. 697–703] for one-dimensional heat conduction with the sources placed outside the space domain of interest, to the two-dimensional setting. Theoretical properties of the method, as well as numerical investigations, are included, showing that accurate results can be obtained efficiently with small computational cost.

Two-dimensional viscoelastic discrete triangular system with negative-stiffness components

Potential applications of high-damping and high-stiffness composites have motivated extensive research on the effects of negative-stiffness inclusions on the overall properties of composites. Recent theoretical advances have been based on the Hashin-Shtrikman composite models, one-dimensional discrete viscoelastic systems and a two-dimensional nested triangular viscoelastic network. In this paper, we further analyze the two-dimensional triangular structure containing pre-selected negative-stiffness components to study its underlying deformation mechanisms and stability. Major new findings are structure-deformation evolution with respect to the magnitude of negative stiffness under shear loading and the phenomena related to dissipation-induced destabilization and inertia-induced stabilization, according to Lyapunov stability analysis. The evolution shows strong correlations between stiffness anomalies and deformation modes. Our stability results reveal that stable damping peaks, i.e. stably extreme effective damping properties, are achievable under hydrostatic loading when the inertia is greater than a critical value. Moreover, destabilization induced by elemental damping is observed with the critical inertia. Regardless of elemental damping, when the inertia is less than the critical value, a weaker system instability is identified.

Three-Dimensional Operational Calculi for Nonlocal Evolution Boundary Value Problems

Иван Христов Димовски, Юлиан Цанков Цанков - Построени са директни операционни смятания за функции u(x, y, t), непрекъснати в област от вида D = [0, a] × [0, b] × [0, ∞). Наред с класическата дюамелова конволюция, построението използва и две некласически конволюции за операторите ∂2x и ∂2y. Тези три едномерни конволюции се комбинират в една тримерна конволюция u ∗ v в C(D). Вместо подхода на Я. Микусински, основаващ се на конволюционни частни, се развива алтернативен подход с използване на мултипликаторните частни на конволюционната алгебра (C(D), ∗).

Ionized cluster beam deposition: Modeling, cluster size measurement and applications

Clusters are aggregations of atoms or molecules, generally intermediate in size between individual atoms and aggregates that are large enough to be called bulk matter. Clusters can also be called nanoparticles, because their size is on the order of nanometers or tens of nanometers. A new field has begun to take shape called nanostructured materials which takes advantage of these atom clusters. The ultra-small size of building blocks leads to dramatically different properties and it is anticipated that such atomically engineered materials will be able to be tailored to perform as no previous material could.^ The idea of ionized cluster beam (ICB) thin film deposition technique was first proposed by Takagi in 1972. It was based upon using a supersonic jet source to produce, ionize and accelerate beams of atomic clusters onto substrates in a vacuum environment. Conditions for formation of cluster beams suitable for thin film deposition have only recently been established following twenty years of effort. Zinc clusters over 1,000 atoms in average size have been synthesized both in our lab and that of Gspann. More recently, other methods of synthesizing clusters and nanoparticles, using different types of cluster sources, have come under development.^ In this work, we studied different aspects of nanoparticle beams. The work includes refinement of a model of the cluster formation mechanism, development of a new real-time, in situ cluster size measurement method, and study of the use of ICB in the fabrication of semiconductor devices.^ The formation process of the vaporized-metal cluster beam was simulated and investigated using classical nucleation theory and one dimensional gas flow equations. Zinc cluster sizes predicted at the nozzle exit are in good quantitative agreement with experimental results in our laboratory.^ A novel in situ real-time mass, energy and velocity measurement apparatus has been designed, built and tested. This small size time-of-flight mass spectrometer is suitable to be used in our cluster deposition systems and does not suffer from problems related to other methods of cluster size measurement like: requirement for specialized ionizing lasers, inductive electrical or electromagnetic coupling, dependency on the assumption of homogeneous nucleation, limits on the size measurement and non real-time capability. Measured ion energies using the electrostatic energy analyzer are in good accordance with values obtained from computer simulation. The velocity (v) is measured by pulsing the cluster beam and measuring the time of delay between the pulse and analyzer output current. The mass of a particle is calculated from m = (2E/v$\sp2).$ The error in the measured value of background gas mass is on the order of 28% of the mass of one N$\sb2$ molecule which is negligible for the measurement of large size clusters. This resolution in cluster size measurement is very acceptable for our purposes.^ Selective area deposition onto conducting patterns overlying insulating substrates was demonstrated using intense, fully-ionized cluster beams. Parameters influencing the selectivity are ion energy, repelling voltage, the ratio of the conductor to insulator dimension, and substrate thickness. ^

Nonlinear quantum transport in low-dimensional electronic devices

The study of transport processes in low-dimensional semiconductors requires a rigorous quantum mechanical treatment. However, a full-fledged quantum transport theory of electrons (or holes) in semiconductors of small scale, applicable in the presence of external fields of arbitrary strength, is still not available. In the literature, different approaches have been proposed, including: (a) the semiclassical Boltzmann equation, (b) perturbation theory based on Keldysh's Green functions, and (c) the Quantum Boltzmann Equation (QBE), previously derived by Van Vliet and coworkers, applicable in the realm of Kubo's Linear Response Theory (LRT). ^ In the present work, we follow the method originally proposed by Van Wet in LRT. The Hamiltonian in this approach is of the form: H = H 0(E, B) + λV, where H0 contains the externally applied fields, and λV includes many-body interactions. This Hamiltonian differs from the LRT Hamiltonian, H = H0 - AF(t) + λV, which contains the external field in the field-response part, -AF(t). For the nonlinear problem, the eigenfunctions of the system Hamiltonian, H0(E, B), include the external fields without any limitation on strength. ^ In Part A of this dissertation, both the diagonal and nondiagonal Master equations are obtained after applying projection operators to the von Neumann equation for the density operator in the interaction picture, and taking the Van Hove limit, (λ → 0, t → ∞, so that (λ2 t)n remains finite). Similarly, the many-body current operator J is obtained from the Heisenberg equation of motion. ^ In Part B, the Quantum Boltzmann Equation is obtained in the occupation-number representation for an electron gas, interacting with phonons or impurities. On the one-body level, the current operator obtained in Part A leads to the Generalized Calecki current for electric and magnetic fields of arbitrary strength. Furthermore, in this part, the LRT results for the current and conductance are recovered in the limit of small electric fields. ^ In Part C, we apply the above results to the study of both linear and nonlinear longitudinal magneto-conductance in quasi one-dimensional quantum wires (1D QW). We have thus been able to quantitatively explain the experimental results, recently published by C. Brick, et al., on these novel frontier-type devices. ^

Formation of Polycyclic Aromatic Hydrocarbons and Nitrogen Containing Polycyclic Aromatic Compounds in Titan's Atmosphere, the Interstellar Medium and Combustion

Several different mechanisms leading to the formation of (substituted) naphthalene and azanaphthalenes were examined using theoretical quantum chemical calculations. As a result, a series of novel synthetic routes to Polycyclic Aromatic Hydrocarbons (PAHs) and Nitrogen Containing Polycyclic Aromatic Compounds (N-PACs) have been proposed. On Earth, these aromatic compounds originate from incomplete combustion and are released into our environment, where they are known to be major pollutants, often with carcinogenic properties. In the atmosphere of a Saturn's moon Titan, these PAH and N-PACs are believed to play a critical role in organic haze formation, as well as acting as chemical precursors to biologically relevant molecules. The theoretical calculations were performed by employing the ab initio G3(MP2,CC)/B3LYP/6-311G** method to effectively probe the Potential Energy Surfaces (PES) relevant to the PAH and N-PAC formation. Following the construction of the PES, Rice-Ramsperger-Kassel-Markus (RRKM) theory was used to evaluate all unimolecular rate constants as a function of collision energy under single-collision conditions. Branching ratios were then evaluated by solving phenomenological rate expressions for the various product concentrations. The most viable pathways to PAH and N-PAC formation were found to be those where the initial attack by the ethynyl (C2H) or cyano (CN) radical toward a unsaturated hydrocarbon molecule led to the formation of an intermediate which could not effectively lose a hydrogen atom. It is not until ring cyclization has occurred, that hydrogen elimination leads to a closed shell product. By quenching the possibility of the initial hydrogen atom elimination, one of the most competitive processes preventing the PAH or N-PAC formation was avoided, and the PAH or N-PAC formation was allowed to proceed. It is concluded that these considerations should be taken into account when attempting to explore any other potential routes towards aromatic compounds in cold environments, such as on Titan or in the interstellar medium.

Nonlinear quantum transport in low-dimensional electronic devices

The study of transport processes in low-dimensional semiconductors requires a rigorous quantum mechanical treatment. However, a full-fledged quantum transport theory of electrons (or holes) in semiconductors of small scale, applicable in the presence of external fields of arbitrary strength, is still not available. In the literature, different approaches have been proposed, including: (a) the semiclassical Boltzmann equation, (b) perturbation theory based on Keldysh's Green functions, and (c) the Quantum Boltzmann Equation (QBE), previously derived by Van Vliet and coworkers, applicable in the realm of Kubo's Linear Response Theory (LRT). In the present work, we follow the method originally proposed by Van Vliet in LRT. The Hamiltonian in this approach is of the form: H = H°(E, B) + λV, where H0 contains the externally applied fields, and λV includes many-body interactions. This Hamiltonian differs from the LRT Hamiltonian, H = H° - AF(t) + λV, which contains the external field in the field-response part, -AF(t). For the nonlinear problem, the eigenfunctions of the system Hamiltonian, H°(E, B) , include the external fields without any limitation on strength. In Part A of this dissertation, both the diagonal and nondiagonal Master equations are obtained after applying projection operators to the von Neumann equation for the density operator in the interaction picture, and taking the Van Hove limit, (λ → 0 , t → ∞ , so that (λ2 t)n remains finite). Similarly, the many-body current operator J is obtained from the Heisenberg equation of motion. In Part B, the Quantum Boltzmann Equation is obtained in the occupation-number representation for an electron gas, interacting with phonons or impurities. On the one-body level, the current operator obtained in Part A leads to the Generalized Calecki current for electric and magnetic fields of arbitrary strength. Furthermore, in this part, the LRT results for the current and conductance are recovered in the limit of small electric fields. In Part C, we apply the above results to the study of both linear and nonlinear longitudinal magneto-conductance in quasi one-dimensional quantum wires (1D QW). We have thus been able to quantitatively explain the experimental results, recently published by C. Brick, et al., on these novel frontier-type devices.

Síntese e caracterização de ligantes nitrogenados derivados de isatina, 4-benziltiossemicarbazida e seus complexos com metais d¹º

Neste trabalho são apresentadas e discutidas as estruturas cristalinas e moleculares do ligante (1), isatina-3-(toluilsulfono-hidrazona), dos complexos [bis(2-acetilpiridina-N4 - benziltiossemicarbazona-N,N,S)Cd(II)], (1), [bis (isatina-3-N4 -benziltiossemicarbazonaN,S)Hg(II)].Etanol, (2) e [bis (isatina-3-N4 -benziltiossemicarbazona-N,S,O)Zn(II)].DMF, (3). Cristais amarelos vítreos do ligante (1) foram obtidos a partir da evaporação lenta de etanol do ensaio de cristalização. Seus dados cristalográficos indicam que duas moléculas interagem através de ligações de hidrogênio do tipo N1-H···O1, formando unidades dímeras. A reação entre 2-acetilpiridina-N4 -benziltiossemicarbazona e Cd(CH3COO)2.2H2O, em presença de etanol, KOH, e após evaporação lenta da mistura de acetona e DMF(2:1), resultou em cristais amarelos do complexo (1). As interações do tipo C(10)-H(10)···S(1)···H(1)-N(4), e N(8)- H(29)⋅⋅⋅S(2) permitem a dimerização do complexo, e a formação de uma cadeia unidimensional. Os cristais laranja do complexo (2) foram obtidos da reação entre o ligante isatina-3-N4 -benziltiossemicarbazona e Hg(NO3)2.H2O, na presença de metanol, KOH, e após evaporação lenta de uma mistura de tolueno e acetona (2:1). As moléculas do complexo (2) estão associadas por ligações de hidrogênio do tipo N(63)-H(4)···O(21), essas interações centrossimétricas conduzem a formação de dímeros. A reação entre o ligante isatina-3-N4 - benziltiossemicarbazona e Zn(CH3COO)2.2H2O, em presença de etanol e KOH resultou em cristais de coloração laranja do complexo (3). A estrutura do complexo apresenta múltiplas ligações de hidrogênio, com formação de dímeros através das interações N1-H1···O1 e C3- H3···N(7). Os dímeros associam-se por interações N4-H4···O2 numa cadeia unidimensional ao longo da direção cristalográfica [100]. A polimerização bidimensional é observada 7 considerando-se as interações do tipo C20-H20···S1, N8-H8···S2 ao longo da direção cristalográfica [010], bem como das interações, N5-H5···O3DMF e C31-H31B···Car, que ocorrem através da molécula de solvente DMF.

Novel laser-induced luminescence resulting from benzophenone/O-propylated p-tert-butylcalix[4]arene complexes. A diffuse reflectance study

Laser-induced room temperature luminescence of air-equilibrated benzophenone/O-propylated p-tert-butylcalix[ 4] arene solid powdered samples revealed the existence of a novel emission, in contrast with benzophenone/p-tertbutylcalix[ 4] arene complexes, where only benzophenone emits. This novel emission was identified as phosphorescence of 1-phenyl-1,2-propanedione, which is formed as the result of an hydrogen atom abstraction reaction of the triplet excited benzophenone from the propoxy substituents of the calixarene. Room temperature phosphorescence was obtained in air-equilibrated samples in all propylated hosts. The decay times of the benzophenone emission vary greatly with the degree of propylation, the shortest lifetimes being obtained in the tri- and tetrapropylated calixarenes. Triplet - triplet absorption of benzophenone was detected in all cases, and is the predominant absorption in the p-tert-butylcalix[ 4] arene case, where an endo-calix complex is formed. Benzophenone ketyl radical formation occurs with the O-propylated p-tert-butylcalix[ 4] arenes hosts, suggesting a different type of host/guest molecular arrangement. Diffuse reflectance laser. ash photolysis and gas chromatography - mass spectrometry techniques provided complementary information, the former about transient species and the latter regarding the final products formed after light absorption. Product analysis and identification clearly show that the two main degradation photoproducts following laser excitation in the propylated substrates are 1-phenyl-1,2- propanedione and 2- hydroxybenzophenone, although several other minor photodegradation products were identified. A detailed mechanistic analysis is proposed. While the solution photochemistry of benzophenone is dominated by the hydrogen abstraction reaction from suitable hydrogen donors, in these solid powdered samples, the alpha-cleavage reaction also plays an important role. This finding occurs even with one single laser pulse which lasts only a few nanoseconds, and is apparently related to the fact that scattered radiation exists, due to multiple internal reflections possibly trapping light within non-absorbing microcrystals in the sample, and is detected until at least 20 mus after the laser pulse. This could explain how photoproducts thus formed could also be excited with only one laser pulse.

Characterization of the volatile profile of Brazilian Merlot wines through comprehensive two dimensional gas chromatography time-of-flight mass spectrometric detection.

Wine aroma is an important characteristic and may be related to certain specific parameters, such as raw material and production process. The complexity of Merlot wine aroma was considered suitable for comprehensive two-dimensional gas chromatography (GCGC), as this technique offers superior performance when compared to one-dimensional gas chromatography (1D-GC). The profile of volatile compounds of Merlot wine was, for the first time, qualitatively analyzed by HS-SPME-GCxGC with a time-of-flight mass spectrometric detector (TOFMS), resulting in 179 compounds tentatively identified by comparison of experimental GCxGC retention indices and mass spectra with literature 1D-GC data and 155 compounds tentatively identified only by mass spectra comparison. A set of GCGC experimental retention indices was also, for the first time, presented for a specific inverse set of columns. Esters were present in higher number (94), followed by alcohols (80), ketones (29), acids (29), aldehydes (23), terpenes (23), lactones (16), furans (14), sulfur compounds (9), phenols (7), pyrroles (5), C13-norisoprenoids (3), and pyrans (2). GCxGC/TOFMS parameters were improved and optimal conditions were: a polar (polyethylene glycol)/medium polar (50% phenyl 50% dimethyl arylene siloxane) column set, oven temperature offset of 10ºC, 7 s as modulation period and 1.4 s of hot pulse duration. Co-elutions came up to 138 compounds in 1D and some of them were resolved in 2D. Among the coeluted compounds, thirty-three volatiles co-eluted in both 1D and 2D and their tentative identification was possible only due to spectral deconvolution. Some compounds that might have important contribution to aroma notes were included in these superimposed peaks. Structurally organized distribution of compounds in the 2D space was observed for esters, aldehydes and ketones, alcohols, thiols, lactones, acids and also inside subgroups, as occurred with esters and alcohols. The Fischer Ratio was useful for establishing the analytes responsible for the main differences between Merlot and non-Merlot wines. Differentiation among Merlot wines and wines of other grape varieties were mainly perceived through the following components: ethyl dodecanoate, 1-hexanol, ethyl nonanoate, ethyl hexanoate, ethyl decanoate, dehydro-2-methyl-3(2H)thiophenone, 3-methyl butanoic acid, ethyl tetradecanoate, methyl octanoate, 1,4 butanediol, and 6-methyloctan-1-ol.

Development of multicomponent models for the prediction of lubricant oil-fuel dilution in internal combustion engines

The current environmental crisis is forcing the automotive industry to face tough challenges for the Internal Combustion Engines development in order to reduce the emissions of pollutants and Greenhouse gases. In this context, in the last decades, the main technological solutions adopted by the manufacturers have been the direct injection and the engine downsizing, which led to the rising of new concerns related to the fuel-cylinder walls physical interaction. The fuel spray possibly impacts the cylinder liner wall, which is wetted by the lubricant oil thus causing the derating of the lubricant properties, increasing the oil consumption, and contaminating the lubricant oil in the crankcase. Also, concerning hydrogen fuelled internal combustion engines, it is likely that the high near-wall temperature, which is typical of the hydrogen flame, results in the evaporation of a portion of the lubricant oil, increasing its consumption. With regards on the innovative combustion systems and their control strategies, optical accessible engines are fundamental tools for experimental investigations on such combustion systems. Though, due to the optical measurement line, optical engines suffer from a high level of blow-by, which must be accounted for. In light of the above, this thesis work aims to develop numerical methodologies with the aim to build useful tools for supporting the design of modern engines. In particular, a one-dimensional modelling of the lubricant oil-fuel dilution and oil evaporation has been performed and coupled with an optimization algorithm to achieve a lubricant oil surrogate. Then, a quasi-dimensional blow-by model has been developed and validated against experimental data. Such model, has been coupled with CFD 3D simulations and directly implemented in CFD 3D. Finally, CFD 3D simulations coupled with the VOF method have been performed in order to validate a methodology for studying the impact of a liquid droplet on a solid surface.