CONTEXTThis research is based on the theoretical study of seven push-pull molecules composed of conjugated bridges based on two different organometallic rings, these bridges are linked at their ends by acceptor groups (-NO2) and donor groups (-N(CH3)2) on the α position of the rings mentioned above. The location of the donor and acceptor groups revealed that the addition of the acceptor groups near the rings (Titanol, Ferrol and Nickelol) improves the NLO response in comparison with the grafting of these groups on the Zinkol ring and also influences the positioning of the π electrons at the level of the chromophores studied. The molecule 2B gave the highest values of static first hyperpolarisabilitiy (βtot) and static second hyperpolarisabilitiy (γav), knowing that: βtot (2B) = 135.79 * 10-30 esu and γav (2B) = 135.79 * 10-35 esu. The highest values of dynamic first __-mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"-__ β ‖ λ - 2 ω ; ω , ω and second __-mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"-__ γ ‖ λ - 2 ω ; ω , ω , 0 hyperpolarisabilities are assigned to the molecule 1C with the following values: __-mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"-__ γ ‖ λ - 2 ω ; ω , ω , 0 =1,218,310.00 * 10-30 esu and __-mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"-__ γ ‖ λ - 2 ω ; ω , ω , 0 =1,324,520,000 * 10-35 esu. The metal Zn is considered as an acceptor group and the remaining metals (Ti, Fe and Ni) are considered as donor groups. The specific solvents for the seven molecules are water, ethanol and acetonitrile. The maximum wavelengths recorded for all molecules in combination with all solvents are in the range of 421.39 to 765.28 nm. λ METHOD: The calculations were performed using Gaussian 16 software to perform DFT calculations with B3LYP functional. The LanL2DZ basis-set was used for transition metals, while the 6-31 + + G(d,p) basis-set was used for nonmetal atoms. The functionals used are: CAM-B3LYP, LC-wPBE, LC-BLYP, M11, wB97X, M08-HX, M06-2X, MN12SX, MN15, and M06HF. The basis-sets used are: 6-31G(d,p), 6-31 + + G(d,p), cc-pVDZ, aug-cc-pVDZ, 6-311G(d,p), 6-311 + + G(d,p), cc-pVTZ, and aug-cc-pVTZ. The Natural Bond Orbital (NBO) calculations are performed by the NBO program incorporated by default in the Gaussian 16 program. The solvation models studied are the CPCM model (conductor polarizable continuum model) and the SMD model (Solvation Model Density). Excited states calculations are calculated by the time-dependent DFT method (TD-DFT).