The advancement in the eld of cold-atoms has generated a lot of interest in the condensed matter community. Cold-atom experiments can simulate clean, disor-der/impurity free systems very easily. In these systems, we have a control over various parameters like tuning the interaction between particles by the Feshbach resonance, tuning the hopping between lattice sites by laser intensity and so on. As a result, these systems can be used to mimic various theoretical models, which was hindered because of various experimental limitations. Thus, we have an ex-perimental tool in which we can start with a simple theoretical model and later tune the model experimentally and theoretically to simulate the real materials. This will be helpful in studying the physics of the real materials as we can control interactions as well as the impurities can also be taken care of. But the advance-ment in the eld of cold atoms has seen a roadblock for the fermions in optical lattices. The super uid and anti-ferromagnetic phases has not been achieved for fermions in optical lattices due to the \cooling problem" (entropy issues). In this thesis, we have addressed the issue of the \cooling problem" for fermions in optical lattice systems and studied the system with determinant quantum Monte Carlo technique. We start by giving a general idea of cold-atoms and optical lat-tice potentials, and a brief review of the experimental work going on in the cold-atomic systems. Experimental limitations like \fermionic cooling problem" have been discussed in some detail. Then we proposed a bilayer band-insulator model to circumvent the \entropy problem" and simultaneously increasing the transi-tion temperature for fermions in optical lattices. We have studied the attractive Hubbard model, which is the minimal model for fermions in optical lattices. The techniques that we have used to study the model are mean- eld theory, Gaussian uctuation theory and determinant quantum Monte Carlo numerical technique. . Chapter-1 : provides a general introduction to the ultra-cold atoms, optical lattice and Feshbach resonance. In this chapter we have discussed about cold-atom experiments in optical lattice systems. Here, we have brie y discussed the control over various parameters in the experiments. The goal of these experiments is to realize or mimic many many-body Hamiltonians in experiments, which until now was just a theoretical tool to describe various many-body physics. In the end we give a brief idea for introducing disorder in the cold-atom experiments discuss the limitations of these experiments in realizing the \interesting" super uid and anti-ferromagnetic phases of fermionic Hubbard model in optical lattices. Chapter-2 : gives a brief idea of \Determinant Quantum Monte-Carlo" (DQM C) technique that has been used to study these systems. In this chapter we will discuss the DQM C algorithm and the observables that can be calculated. We will discuss certain limitation of the DQM C algorithm like numerical instability and sign problem. We will brie y discuss how sign problem doesn't occur in the model we studied. Chapter-3 : discusses the way by which we can bypass the \cooling problem" (high entropy state) to get a fermionic super uid state in the cold atom experi-ments. In this chapter we propose a model whose idea hinges on a low-entropy band-insulator state, which can be tuned to super uid state by tuning the on-site attractive interaction by Feshbach resonance. We show through Gaussian uctua-tion theory that the critical temperature achieved is much higher in our model as compared to the single-band Hubbard model. Through detailed variational Monte Carlo calculations, we have shown that the super uid state is indeed the most stable ground state and there is no other competing order. In the end we give a proposal for its realization in the ultra-cold atom optical lattice systems. Chapter-4 : discusses the DQM C study of the model proposed in chapter- 3. Here we have studied the various single-particle properties like momentum distribution, double occupancies which can be easily measured in cold-atom ex-periments. We also studied the pair-pair and the density-density correlations in detail through DQM C algorithm and mapped out the full T U phase diagram. We show that the proposed model doesn't favor the charge density wave for the interaction strengths we are interested in. Chapter-5 : gives a brief idea of the e ect of adding an on-site random disorder in our proposed bilayer-Hubbard model. We study the e ect of random disorder on various single-particle properties which can be easily veri ed in cold-atom ex-periments. We studied the suppression of the pair-pair correlations as we increase the disorder strength in our proposed model. We nd that the critical value of the interaction doesn't change in the weak-disorder limit. We estimated the critical disorder strength needed to destroy the super uid state and argued that the tran-sition from the super uid to Bose-glass phase in presence of disorder lies in the universality class of (d + 1) XY model. In the end, we give a schematic U V phase diagram for our system. Chapter-6 : We studied the bilayer attractive Hubbard model in different lattice geometry, the bilayer honeycomb lattice, both in presence and in absence of the on-site random disorder. We discussed how the pair-pair and density-density cor-relations behave in the presence and absence of disorder. Through the finite-size scaling analysis we see the co-existence of the super fluid and the charge density wave order at half- lling. An in nitesimal disorder destroys the CDW order com-pletely while the super uid phase found to be robust against weak-disorder. We estimated the critical interaction strength, the critical temperature and the critical disorder strength through nite-size scaling, and provide a putative phase diagram for the system considered.