Influence of the interaction potential on simulated sputtering and reflection data

The TRIM.SP program which is based on the binary collision approximation was changed to handle not only repulsive interaction potentials, but also potentials with an attractive part. Sputtering yields, average depth and reflection coefficients calculated with four different potentials are compared. Three purely repulsive potentials (Meliere, Kr-C and ZBL) are used and an ab initio pair potential, which is especially calculated for silicon bombardment by silicon. The general trends in the calculated results are similar for all potentials applied, but differences between the repulsive potentials and the ab initio potential occur for the reflection coefficients and the sputtering yield at large angles of incidence.

Relativistic prediction of the ground state of atomic Lawrencium

In contradiction to the prediction of the Periodic Table but in agreement with earlier suggestions by Brewer and Mann, the ground state configuration of atomic Lawrencium (Z = 103) will not be 7s^2 6d^2 D_3/2 but 7s^2 7p ^2p_1/2. The reason for this deviation from normal trends across the Periodic Table are strong relativistic effects on the outermost 7P_l/2 orbital. Multicontiguration Dirac-Fock calculations are reported for Lawrencium and analogous lighter atoms. These calculations include contributions from magnetic and retardation interactions and an estimation of quantum electrodynamic corrections.

Relativistic atomic level calculations with nuclei carrying additional fractional charges

One-electron energy levels and wavelengths have been calculated for Na-like ions whose nuclei carry quarks with additional charges ±e/3, ±2e/3. The calculations are based on relativistic self-consistent field procedures. The deviations from experimental values exhibit regularities which allow an extrapolation for the wavelengths of 3s - 3p, 3s - 4p, 3p - 3d, and 3p - 4s transitions for the nuclear charge Z = 11± 1/3, ±2/3. A number of transitions are found in the region of visible light which could be used in an optical search for quark atoms.

Finite element method for the accurate solution of diatomic molecules

We present the Finite-Element-Method (FEM) in its application to quantum mechanical problems solving for diatomic molecules. Results for Hartree-Fock calculations of H_2 and Hartree-Fock-Slater calculations of molecules like N_2 and C0 have been obtained. The accuracy achieved with less then 5000 grid points for the total energies of these systems is 10_-8 a.u., which is demonstrated for N_2.

Atomic and ionic radii of superheavy elements

Atomic and ionic radii are presented for the elements E104-E120 and E156-E172. It is shown that a number of effects correlated with the large relativistic contraction of orbitals with low angular momentum leads to smaller atoms for higher atomic numbers. It is expected that Cs is the largest atom in nature.

Electronic configuration in the ground state of atomic lawrencium

Self-consistent relativistic Dirac-Hartree-Fock calculations have been made of some lowlying electronic energies for the atoms of all elements in ground-state ds^2 electron configurations. The results indicate that, contrary to some previous estimates, the ground electronic state of atomic Lr could be in either the 5f^14 6d7s^2 or the 5f^14 7p 7s^2 electron configuration. The separation between the lowest energy level of the 5f^14 6d7s^2 configuration and the lowest energy level of the 5f^14 7p7s^2 configuration is estimated to be (0 ± 3) x 10^3 cm^-1 for atomic Lr.

Observation of two-electron-one-photon transitions in silicon

We report on the observation of K\alpha\alpha X-rays of Si, produced in collisions of 15-28 MeV Si projectiles with various target atoms in the range Z =6 to 29. Energy shifts of X-rays were measured and are compared with theoretical predictions. Cross section ratios for emission of K\alpha\alpha and K\alpha radiation are given.

Molecular effects on the cross section for pair production and their importance in the evaluation of the total photonuclear cross section of {^16}O in the {\delta-resonance} region

The screening correction to the coherent pair-production cross section on the oxygen molecule has been calculated using self-consistent relativistic wave functions for the one-center and two-center Coulomb potentials. It is shown that the modification of the wave function due to molecular binding and the interference between contributions from the two atoms have both sizeable effects on the screening correction. The so-obtained coherent pair-production cross section which makes up the largest part of the total atomic cross section was used to evaluate the total nuclear absorption cross section from photon attenuation measurements on liquid oxygen. The result agrees with cross sections for other nuclei if A-scaling is assumed. The molecular effect on the pair cross section amounts to 15 % of the nuclear cross section in the {\delta-resonance} region.

Time-resolved studies of neutral and ionized Na_n clusters with femtosecond light pulses

The real-time dynamics of Na_n (n=3-21) cluster multiphoton ionization and fragmentation has been studied in beam experiments applying femtosecond pump-probe techniques in combination with ion and electron spectroscopy. Three dimensional wave packet motions in the trimer Na_3 ground state X and excited state B have been observed. We report the first study of cluster properties (energy, bandwidth and lifetime of intermediate resonances Na_n^*) with femtosecond laser pulses. The observation of four absorption resonances for the cluster Na_8 with different energy widths and different decay patterns is more difficult to interpret by surface plasmon like resonances than by molecular structure and dynamics. Timeresolved fragmentation of cluster ions Na_n^+ indicates that direct photo-induced fragmentation processes are more important at short times than the statistical unimolecular decay.

Goldnanoteilchen auf Titandioxid

In dieser Arbeit wurde das Wachstum sowie die ultraschnelle Elektronendynamik des Oberflächenplasmon Polaritons von Goldnanoteilchen auf Titandioxid untersucht. Die Messung der Dephasierungszeit des Oberflächenplasmons von Nanoteilchen mit definierter Form und Größe erfolgte dabei mit der Methode des spektralen Lochbrennens. Die Nanoteilchen wurden durch Deposition von Goldatomen aus einem thermischen Atomstrahl mit anschließender Diffussion und Nukleation, d.h. Volmer-Weber-Wachstum, auf Titandioxidsubstraten hergestellt und mittels einer Kombination aus optischer Spektroskopie und Rasterkraftmikroskopie systematisch untersucht. Dabei lässt sich das Nanoteilchenensemble durch das mittlere Achsverhältnis und den mittleren Äquivalentradius charakterisieren. Die Messungen zeigen, dass die Proben große Größen- und Formverteilungen aufweisen und ein definierter Zusammenhang zwischen Größe und Form der Teilchen existiert. Während kleine Goldnanoteilchen nahezu kugelförmig sind, flachen die Teilchen mit zunehmender Größe immer mehr ab. Des Weiteren wurde in dieser Arbeit die Methode des lasergestützten Wachstums auf das System Gold auf Titandioxid angewendet. Systematische Untersuchungen zeigten, dass sich das Achsverhältnis der Teilchen durch geeignete Wahl von Photonenenergie und Fluenz des eingestrahlten Laserlichts definiert und gezielt vorgeben lässt. Die Methode des lasergestützten Wachstums erschließt damit den Bereich außerhalb der Zugänglichkeit des natürlichen Wachstums. Aufgrund der Formabhängigkeit der spektrale Lage der Plasmonresonanz ist man somit in der Lage, die optischen Eigenschaften der Nanoteilchen gezielt einzustellen und z.B. für technische Anwendungen zu optimieren. Die Untersuchung der ultraschnellen Elektronendynamik von Goldnanoteilchen auf Titandioxid mit äquivalenten Radien zwischen 8 bis 15 nm erfolgte in dieser Arbeit mit der Methode des spektralen Lochbrennes. Hierzu wurde die Dephasierungszeit des Oberflächenplasmons systematisch als Funktion der Photonenenergie in einem Bereich von 1,45 bis 1,85 eV gemessen. Es zeigte sich, dass die gemessenen Dephasierungszeiten von 8,5 bis 16,2 fs deutlich unter den in der dielektrischen Funktion von Gold enthaltenen Werten lagen, was den erwarteten Einfluss der reduzierten Dimension der Teilchen demonstriert. Um die Messwerte trotz verschiedener Teilchengrößen untereinander vergleichen und den Einfluss der intrinsischen Dämpfung quantifizieren zu können, wurde zusätzlich der Dämpfungsparameter A bestimmt. Die ermittelten A-Faktoren zeigten dabei eine starke Abhängigkeit von der Plasmonenergie. Für Teilchen mit Plasmonenergien von 1,45 bis 1,55 eV wurde ein Dämpfungsfaktor von A ~ 0,2 nm/fs ermittelt, der lediglich Oberflächenstreuung als dominierenden Dämpfungsmechanismus widerspiegelt. Hingegen wurde für Teilchen mit Plasmonenergien oberhalb von 1,55 eV ein drastischer Anstieg der Dämpfung auf A ~ 0,4 nm/fs beobachtet. Die erhöhte Dämpfung wurde dabei dem zusätzlichen Vorliegen einer chemischen Dämpfung durch das Titandioxidsubstrat zugeschrieben. Zusammenfassend zeigen die Ergebnisse somit, dass eine starke Abhängigkeit der chemischen Dämpfung von der Photonenenergie vorliegt. Es konnte erstmals nachgewiesen werden, dass die chemische Dämpfung erst ab einer bestimmten unteren Schwelle der Photonenenergie einsetzt, die für Goldnanoteilchen auf Titandioxid bei etwa 1,6 eV liegt.

Relativistic LCAO with Minimax Principle and New Balanced Basis Sets

Relativistic density functional theory is widely applied in molecular calculations with heavy atoms, where relativistic and correlation effects are on the same footing. Variational stability of the Dirac Hamiltonian is a very important field of research from the beginning of relativistic molecular calculations on, among efforts for accuracy, efficiency, and density functional formulation, etc. Approximations of one- or two-component methods and searching for suitable basis sets are two major means for good projection power against the negative continuum. The minimax two-component spinor linear combination of atomic orbitals (LCAO) is applied in the present work for both light and super-heavy one-electron systems, providing good approximations in the whole energy spectrum, being close to the benchmark minimax finite element method (FEM) values and without spurious and contaminated states, in contrast to the presence of these artifacts in the traditional four-component spinor LCAO. The variational stability assures that minimax LCAO is bounded from below. New balanced basis sets, kinetic and potential defect balanced (TVDB), following the minimax idea, are applied with the Dirac Hamiltonian. Its performance in the same super-heavy one-electron quasi-molecules shows also very good projection capability against variational collapse, as the minimax LCAO is taken as the best projection to compare with. The TVDB method has twice as many basis coefficients as four-component spinor LCAO, which becomes now linear and overcomes the disadvantage of great time-consumption in the minimax method. The calculation with both the TVDB method and the traditional LCAO method for the dimers with elements in group 11 of the periodic table investigates their difference. New bigger basis sets are constructed than in previous research, achieving high accuracy within the functionals involved. Their difference in total energy is much smaller than the basis incompleteness error, showing that the traditional four-spinor LCAO keeps enough projection power from the numerical atomic orbitals and is suitable in research on relativistic quantum chemistry. In scattering investigations for the same comparison purpose, the failure of the traditional LCAO method of providing a stable spectrum with increasing size of basis sets is contrasted to the TVDB method, which contains no spurious states already without pre-orthogonalization of basis sets. Keeping the same conditions including the accuracy of matrix elements shows that the variational instability prevails over the linear dependence of the basis sets. The success of the TVDB method manifests its capability not only in relativistic quantum chemistry but also for scattering and under the influence of strong external electronic and magnetic fields. The good accuracy in total energy with large basis sets and the good projection property encourage wider research on different molecules, with better functionals, and on small effects.

Ferrocene-Based Pyridylphosphine Ligands – Coordination Chemistry of Group 10, 11 and 12 Metals

English: The present thesis describes the synthesis of 1,1’-ferrocendiyl-based pyridylphosphine ligands, the exploration of their fundamental coordination chemistry and preliminary experiments with selected complexes aimed at potential applications. One main aspect is the synthesis of the bidentate ferrocene-based pyridylphosphine ligands 1-(Pyrid-2-yl)-1’-diphenylphosphinoferrocene, 1-(Pyrid-3-yl)-1’-diphenylphosphinoferrocene and 1-[(Pyrid-2-yl)methyl]-1’-diphenylphosphinoferrocene. A specific feature of these ligands is the ball-bearing like flexibility of the ferrocenebased backbone. An additional flexibility element is the rotation around the C–C single bonds. Consequently, the donor atoms can realise a wide range of positions with respect to each other and are therefore able to adapt to the coordination requirements of different metal centres. The flexibility of the ligand also plays a role in another key aspect of this work, which concerns the coordination mode, i. e. bridging vs. chelating. In addition to the flexibility, also the position of the donor atoms to each other is important. This is largely affected by the position of the pyridyl nitrogen (pyrid-2-yl vs. pyrid-3-yl) and the methylen group in 1-[(Pyrid-2-yl)methyl]-1’-diphenylphosphinoferrocene. Another interesting point is the combination of a soft phosphorus donor atom with a harder nitrogen donor atom, according to the HSAB principle. This combination generates a unique binding profile, since the pi-acceptor character of the P site is able to stabilise a metal centre in a low oxidation state, while the nitrogen sigma-donor ability can make the metal more susceptible to oxidative addition reactions. A P,N-donor combination can afford hemilabile binding profiles, which would be ideal for catalysis. Beyond 1,2-substituted ferrocene derivatives, which are quite successful in catalytic applications, 1,1’-derivatives are rather underrepresented. While a low-yield synthetic pathway to 1-(Pyrid-2-yl)-1’-diphenylphosphinoferrocene was already described in the literature [I. R. Butler, Organometallics 1992, 11, 74.], it was possible to find a new, improved and simplified synthetic pathway. Both other ligands were unknown prior to this work. Satisfactory results in the synthesis of 1-(Pyrid-3-yl)-1’-diphenylphosphinoferrocene could be achieved by working in analogy to the new synthetic procedure for 1-(Pyrid-2-yl)-1’-diphenylphosphinoferrocene. The synthesis of 1-[(Pyrid-2-yl)methyl]-1’-diphenylphosphinoferrocene has been handled by the group of Prof. Petr Stepnicka from Charles University, Prague, Czech Republic. The synthesis of tridentate ligands with an analogous heterodentate arrangement, was investigated briefly as a sideline of this study. The major part of this thesis deals with the fundamental coordination chemistry towards transition metals of the groups 10, 11 and 12. Due to the well-established catalytic properties of analogous palladium complexes, the coordination chemistry towards palladium (group 10) is of particular interest. The metals zinc and cadmium (group 12) are also of substantial importance because they are redox-inert in their divalent state. This is relevant in view of electrochemical investigations concerning the utilisation of the ligands as molecular redox sensors. Also mercury and the monovalent metals silver and gold (group 11) are included because of their rich coordination chemistry. It is essential to answer questions concerning aspects of the ligands’ coordination mode bearing in mind the HSAB principle.

Selbstassemblierte Monolagen und magnetisch strukturierte Submonolagen mit dia- und paramagnetischen Molekülen sowie Einzelmolekülmagneten auf Phthalocyanin- und Subphthalocyanin-Basis

Mit dieser Arbeit wurde die Selbstassemblierung von dia- und paramagnetischen Molekülen sowie Einzelmolekülmagneten auf Goldsubstraten und magnetisch strukturierten Substraten untersucht. Dazu wurden drei verschiedene Klassen an Phthalocyaninderivaten verwendet: Diamagnetische Subphthalocyanine, paramagnetische Phthalocyaninatometalle und Diphthalocyaninatolanthanidkomplexe. Alle synthetisierten Verbindungen sind peripher thioethersubstituiert. Die Alkylketten (a: n-C8H17, b: n-C12H25) vermitteln die Löslichkeit in vielen organischen Solventien und sorgen für eine geordnete Assemblierung auf einer Oberfläche, wobei die Bindung auf Gold hauptsächlich über die Schwefelatome stattfindet. Die aus Lösung abgeschiedenen selbstassemblierten Monolagen wurden mit XPS, NEXAFS-Spektroskopie und ToF-SIMS untersucht. Bei der Selbstassemblierung auf magnetisch strukturierten Substraten stehen die Moleküle unter dem Einfluss magnetischer Streufelder und binden bevorzugt nur in bestimmten Bereichen. Die gebildeten Submonolagen wurden zusätzlich mit X-PEEM untersucht. Die erstmals dargestellten Manganphthalocyanine [MnClPc(SR)8] 1 wurden ausgehend von MnCl2 erhalten. Hier fand bei der Aufarbeitung an Luft eine Oxidation zu Mangan(III) statt; +III ist die stabilste Oxidationsstufe von Mangan in Phthalocyaninen. Der Nachweis des axialen Chloridoliganden erfolgte mit Massenspektrometrie und FIR- sowie Raman-Spektroskopie. SQUID-Messungen haben gezeigt, dass die Komplexe 1 vier ungepaarte Elektronen haben. Bei den Subphthalocyaninen [BClSubpc(SR)6] 2 wurde der axiale Chloridoligand mit dem stäbchenförmigen Phenolderivat 29-H substituiert und die erfolgreiche Ligandensubstitution durch NMR- und IR-Spektroskopie sowie Massenspektrometrie an den Produkten [BSubpc(SR)6(29)] 30 belegt. Der Radikalcharakter der synthetisierten Terbiumkomplexe [Tb{Pc(SR)8}2] 3 wurde spektroskopisch nachgewiesen; SQUID-Messungen ergaben, dass es sich um Einzelmolekülmagnete mit einer Energiebarriere U des Doppelpotentialtopfs von 880 K oder 610 cm-1 bei 3a handelt. Zunächst wurden die SAMs der Komplexverbindungen 1, 2, 30 und 3 auf nicht magnetisch strukturierten Goldsubstraten untersucht. Die Manganphthalocyanine 1 bilden geordnete SAMs mit größtenteils flach liegenden Molekülen, wie die XPS-, NEXAFS- und ToF-SIMS-Analyse zeigte. Die Mehrzahl der Thioether-Einheiten ist auf Gold koordiniert und die Alkylketten zeigen ungeordnet von der Oberfläche weg. Bei der Adsorption findet eine Reduktion zu Mangan(II) statt und der axiale Chloridoligand wird abgespalten. Das beruht auf dem sog. Oberflächen-trans-Effekt. Im vorliegenden Fall übt die Metalloberfläche einen stärkeren trans-Effekt als der axiale Ligand aus, was bisher experimentell noch nicht beobachtet wurde. Die thioethersubstituierten Subphthalocyanine 2 und 30 sowie die Diphthalocyaninatoterbium-Komplexe 3 sind ebenfalls für SAMs geeignet. Ihre Monolagen wurden mit XPS und NEXAFS-Spektroskopie untersucht, und trotz einer gewissen Unordnung in den Filmen liegen die Moleküle jeweils im Wesentlichen flach auf der Goldoberfläche. Vermutlich sind bei diesen Systemen auch die Alkylketten größtenteils parallel zur Oberfläche orientiert. Im Gegensatz zu den Manganphthalocyaninen 1 tritt bei 2b, 30a, 30b und 3b neben der koordinativen Bindung der Schwefelatome auf Gold auch eine für Thioether nicht erwartete kovalente Au–S-Bindung auf, die durch C–S-Bindungsbruch unter Abspaltung der Alkylketten ermöglicht wird. Der Anteil, zu dem dieser Prozess stattfindet, scheint nicht mit der Molekülstruktur zu korrelieren. Selbstassemblierte Submonolagen auf magnetisch strukturierten Substraten wurden mit dem diamagnetischen Subphthalocyanin 2b hergestellt. Der Nachweis der Submonolagen war schwierig und gelang schließlich durch eine Kombination von ToF-SIMS, NEXAFS Imaging und X-PEEM. Die Analyse der ToF-SIMS-Daten zeigte, dass tatsächlich eine Modulation der Verteilung der Moleküle auf einem unterwärts magnetisch strukturierten Substrat eintritt. Mit X-PEEM konnte die magnetische Struktur der ferromagnetischen Schicht des Substrats direkt der Verteilung der adsorbierten Moleküle zugeordnet werden. Die Subphthalocyanine 2b adsorbieren nicht an den Domänengrenzen, sondern vermehrt dazwischen. Auf Substraten mit abwechselnd 6.5 und 3.5 µm breiten magnetischen Domänen binden die Moleküle bevorzugt in den Bereichen geringster magnetischer Streufeldgradienten, also den größeren Domänen. Solche Substrate wurden für die ToF-SIMS- und X-PEEM-Messungen verwendet. Bei größeren magnetischen Strukturen mit ca. 400 µm breiten Domänen, wie sie aufgrund der geringeren Ortsauflösung dieser Methode für NEXAFS Imaging eingesetzt wurden, binden die Moleküle dann in allen Domänen. Die diamagnetischen Moleküle werden nach dieser Interpretation aus dem inhomogenen Magnetfeld über der Probenoberfläche heraus gedrängt und verhalten sich analog makroskopischer Diamagnete. Die eindeutige Detektion der Moleküle auf den magnetisch strukturierten Substraten konnte bisher nur für die diamagnetischen Subphthalocyanine 2b erfolgen. Um die Interpretation ihres Verhaltens bei der Selbstassemblierung in einem inhomogenen Magnetfeld weiter voranzutreiben, wurde das Subphthalocyanin 37b dargestellt, welches ein stabiles organisches TEMPO-Radikal in seinem axialen Liganden enthält. Das paramagnetische Subphthalocyanin 37b sollte auf den magnetisch strukturierten Substraten in Regionen starker magnetischer Streufelder binden und damit das entgegengesetzte Verhalten zu den diamagnetischen Subphthalocyaninen 2b zeigen. Aus Zeitgründen konnte dieser Nachweis im Rahmen dieser Arbeit noch nicht erbracht werden.

Electron-Nuclear Correlation in Laser-Driven Diatomic Molecules

The interaction of short intense laser pulses with atoms/molecules produces a multitude of highly nonlinear processes requiring a non-perturbative treatment. Detailed study of these highly nonlinear processes by numerically solving the time-dependent Schrodinger equation becomes a daunting task when the number of degrees of freedom is large. Also the coupling between the electronic and nuclear degrees of freedom further aggravates the computational problems. In the present work we show that the time-dependent Hartree (TDH) approximation, which neglects the correlation effects, gives unreliable description of the system dynamics both in the absence and presence of an external field. A theoretical framework is required that treats the electrons and nuclei on equal footing and fully quantum mechanically. To address this issue we discuss two approaches, namely the multicomponent density functional theory (MCDFT) and the multiconfiguration time-dependent Hartree (MCTDH) method, that go beyond the TDH approximation and describe the correlated electron-nuclear dynamics accurately. In the MCDFT framework, where the time-dependent electronic and nuclear densities are the basic variables, we discuss an algorithm to calculate the exact Kohn-Sham (KS) potentials for small model systems. By simulating the photodissociation process in a model hydrogen molecular ion, we show that the exact KS potentials contain all the many-body effects and give an insight into the system dynamics. In the MCTDH approach, the wave function is expanded as a sum of products of single-particle functions (SPFs). The MCTDH method is able to describe the electron-nuclear correlation effects as the SPFs and the expansion coefficients evolve in time and give an accurate description of the system dynamics. We show that the MCTDH method is suitable to study a variety of processes such as the fragmentation of molecules, high-order harmonic generation, the two-center interference effect, and the lochfrass effect. We discuss these phenomena in a model hydrogen molecular ion and a model hydrogen molecule. Inclusion of absorbing boundaries in the mean-field approximation and its consequences are discussed using the model hydrogen molecular ion. To this end, two types of calculations are considered: (i) a variational approach with a complex absorbing potential included in the full many-particle Hamiltonian and (ii) an approach in the spirit of time-dependent density functional theory (TDDFT), including complex absorbing potentials in the single-particle equations. It is elucidated that for small grids the TDDFT approach is superior to the variational approach.

Theoretical study of magnetism, structure and chemical order in transition-metal alloy clusters

Research on transition-metal nanoalloy clusters composed of a few atoms is fascinating by their unusual properties due to the interplay among the structure, chemical order and magnetism. Such nanoalloy clusters, can be used to construct nanometer devices for technological applications by manipulating their remarkable magnetic, chemical and optical properties. Determining the nanoscopic features exhibited by the magnetic alloy clusters signifies the need for a systematic global and local exploration of their potential-energy surface in order to identify all the relevant energetically low-lying magnetic isomers. In this thesis the sampling of the potential-energy surface has been performed by employing the state-of-the-art spin-polarized density-functional theory in combination with graph theory and the basin-hopping global optimization techniques. This combination is vital for a quantitative analysis of the quantum mechanical energetics. The first approach, i.e., spin-polarized density-functional theory together with the graph theory method, is applied to study the Fe$_m$Rh$_n$ and Co$_m$Pd$_n$ clusters having $N = m+n \leq 8$ atoms. We carried out a thorough and systematic sampling of the potential-energy surface by taking into account all possible initial cluster topologies, all different distributions of the two kinds of atoms within the cluster, the entire concentration range between the pure limits, and different initial magnetic configurations such as ferro- and anti-ferromagnetic coupling. The remarkable magnetic properties shown by FeRh and CoPd nanoclusters are attributed to the extremely reduced coordination number together with the charge transfer from 3$d$ to 4$d$ elements. The second approach, i.e., spin-polarized density-functional theory together with the basin-hopping method is applied to study the small Fe$_6$, Fe$_3$Rh$_3$ and Rh$_6$ and the larger Fe$_{13}$, Fe$_6$Rh$_7$ and Rh$_{13}$ clusters as illustrative benchmark systems. This method is able to identify the true ground-state structures of Fe$_6$ and Fe$_3$Rh$_3$ which were not obtained by using the first approach. However, both approaches predict a similar cluster for the ground-state of Rh$_6$. Moreover, the computational time taken by this approach is found to be significantly lower than the first approach. The ground-state structure of Fe$_{13}$ cluster is found to be an icosahedral structure, whereas Rh$_{13}$ and Fe$_6$Rh$_7$ isomers relax into cage-like and layered-like structures, respectively. All the clusters display a remarkable variety of structural and magnetic behaviors. It is observed that the isomers having similar shape with small distortion with respect to each other can exhibit quite different magnetic moments. This has been interpreted as a probable artifact of spin-rotational symmetry breaking introduced by the spin-polarized GGA. The possibility of combining the spin-polarized density-functional theory with some other global optimization techniques such as minima-hopping method could be the next step in this direction. This combination is expected to be an ideal sampling approach having the advantage of avoiding efficiently the search over irrelevant regions of the potential energy surface.